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SUMMARY 


A study is made of a new theoretical approach for -the description of 
gas flows having strong deviations from local translational equilibrium 
(sufficient "rarefaction interaction "), when the Navier-Stokes equations 
become invalid. The usual approach in dealing with the limitations in the 
macroscopic description is to consider the kinetic level of description., 
governed by the Boltzmann equation, and perhaps to relate a kinetic analysis 
to higher approximations in the macroscopic level, by the Chapman-Enskog 
procedure. Because of (a) the general intractability of the Boltzmann 
equation or models of it, and (15 ) the limitation of the Chapman -Enskog 
theory to near-local-equilibrium flows, the present approach considers the 
development and possible utility of a directional level of description. 
Properties and equations on the directional level are obtained by integrat- 
ing the molecular properties and the Boltzmann equation over all magnitudes 
of molecular velocity. The resulting dependent variables (properties of a 
class of molecules) then depend on space, time, and molecular-velocity 
direction , but not magnitude. 

The development on the directional level, between the kinetic and 
macroscopic levels, is intended to exploit the significant directional 
aspects of translational nonequilibrium in describing the gas flow in 
regions of rapid variations of the flow variables. It is useful to con- 
sider a directional mean free path , defined as the average distance 
(measured relative to the observer) travelled by a molecule, moving in a 
designated direction, between collisions with other molecules. This free 
path is useful in formulating a physically intuitive model for approximating 
the "gain-term" collision integrals in the equations of change on the 
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directional level. The gain -term collision integrals on the directional 
level are replaced by terms that are related to appropriate corresponding 
"loss terms" evaluated at a point a directional mean free path away and at 
a mean collision time in the past. Several intuitive direotional-average 
collision models are proposed for use in the directional -mean-free -path 
approximat ion . 

Use of the directional -mean-free -path approximation leaves the 
equations on the directional level in a form where use of Lagrange f s 
expansion for further approximation appears natural. For more than one 
independent variable, generalizations of Lagrange r s expansion are needed, 
and simple forms that are directly applicable are derived. A perturbation- 
expansion scheme based on these generalizations is also developed. This 
scheme is then applied to the equations in the directional -mean-free -path 
method, to simplify the theory. 

The use of the directional -mean-free -path method is illustrated by 
an outlined application to the problem of steady-flow shock -wave structure 
in a monatomic gas. The main value of the method would be expected to be 
realized only after an extension of the theory to include the effects of 
boundaries, for application to problems that may not be tractable by more 
detailed kinetic -theory methods. 
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CHAPTER I 


INTRODUCTION 


This study introduces a new theoretical approach to the description 
of gas flow. This approach may eventually lead to the solution of gas- 
dynamic prohlems not accurately governed by the Navier-Stokes or higher- 
order continuum descriptions, and not sufficiently tractable by kinetic- 
theory methods . 

The limitations of the Navier-Stokes equations, and of the higher- 
order approximations of the macroscopic-moment equations on the continuum 
level (obtained, e.g. , by the Chapman-Enskog procedure), are well-known 
in the calculation of gas flows that are too far out of translational 
equilibrium (see, e.g., Grad 1958, i960, 1963, 1966). The usual approach 
is then to consider the kinetic level of description, either by synthetic 
representations of the velocity distribution function, or by the theory 
of the Boltzmann equation or models of it. An analysis on the kinetic 
level may be related to higher approximations on the macroscopic level, 
as by the Chapman-Enskog procedure. 

Synthetic methods of approximating the velocity distribution function 
are represented by: the bimodal-distribution method of Mott-Smith (1951); 

generalizations of this, including a trimodal distribution function, by 
Krook (1959) (see Anderson and Macomber 1965); the "two-stream" and 
"two- fluid" methods , which are somewhat analogous to the bimodal-distribu- 
tion method, by Lees (1959) (see also Liu and Lees 1961) and Glansdorff 
(1961, 1962); and the ellipsoidal-distribution- function method of Holway 
(1965, 1967). These approximate methods are of considerable value. 

Some limitations have been discussed by Liepmann, Narasimha, and Chahine 

( 1966 ). 

Models of the Boltzmann equation, in particular the BGK, or Krook- 
Welander, model (Bhatnagar, Gross, and Krook 195*+; Welander 195*0 have 
been useful. Recently a more sophisticated conceptual extension and 
generalization (of the BGK model), the ellipsoidal-statistical model, has 
been presented by Holway (1966). The validity of the BGK model may be 
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limited in some problems by: (a) its equivalence to replacing the "gain- 

term” collision integral for a given point in the flow by an expression 
that assumes molecules come out of collisions in equilibrium with the gas 
at that point (cf. Liepmann, Narasimha, and Chahine 1962), or (b) its 
equivalence to a "first iteration" of the Boltzmann equation about a 
local-equilibrium velocity distribution (Gross 195*+; also Rott 1964). 

In spite of the possible limitations, the BGK model does retain many 
important aspects of the Boltzmann equation for arbitrary nonequilibrium 
and has been, and will undoubtedly continue to be, extremely useful for 
kinetic-theory investigations. The noted limitations may be important 
in some problems, however. In addition, the model equations are still 
highly nonlinear integro-partial differential equations and are generally 
difficult to deal with. 

The main purpose of the present study, then, is to investigate a 
new approach on what may be called the "directional level of description." 
The directional level deals with molecular quantities, or properties of 
a class of molecules, that have been integrated over, or averaged over, 
all magnitudes of molecular velocity; the dependent variables then all 
may depend on molecular-velocity direction , as well as on configuration- 
space coordinates and time. For each molecular property of interest, 
the appropriate equations are developed by one integration of the 
Boltzmann equation (over all molecular -velocity magnitudes). The treat- 
ment on the directional level is therefore between the kinetic level 
and the macroscopic level. 

Theoretical analysis on the directional level may be more tractable 
than on the kinetic level, and may at the same time be capable of treating 
substantial translational nonequilibrium in gas flows, or substantial 
"rarefaction interaction." Detailed discussions motivating this approach 
are given in Chapters II and III. The basis may be noted here: that in 

many problems the translational nonequilibrium is essentially a directional 
phenomenon. Classic examples are shock waves and boundary layers, which 
are very thin nonequilibrium regions of rapid transition to a local- 
translational-equilibrium state. In those cases, the flow properties 
change most rapidly in a direction normal to the thin region. In 
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accounting for the significant directional aspects of translational 
nonequi librium , the concept of a mean free path that varies with direction 
is useful. 

In the most general form, the equations on the directional level are 
expected to he appropriate for approximate calculation of flows in 
arbitrary translational nonequilibrium because of the allowance for the 
variation of the flow properties with molecular-velocity direction. 

This study of the new approach on the directional level may be 
considered to be preliminary. At various steps in the development, 
intuitive approximations are used; and the formulation is limited to 
apply only to sufficiently rarefied ("perfect") gases of similar, neutral, 
spherically symmetric, monatomic molecules. Also the effects of "surfaces," 
or condensed regions, adjacent to the gas flow have not been adequately 
considered. (The theory could in the future be much more highly developed, 
if such further developments appear justified. ) 

Chapters II and III include appropriate definitions and interpreta- 
tions that are intended to both motivate and lay a formal rigorous 
foundation for the treatment on the directional level in terms of the 
directional mean free path. 

Then in Chapter IV the directional-mean- free-rath approximation for the 
gain-term collision integrals is introduced, along with several proposed intui- 
tive "directional-average collision models"; and a scheme is outlined for 
making the system of equations on the directional level determined. The 
postulated intuitive models involved in the method are to be partially 
justified physically, but final justification for use of any model will 
lie in the tractability and utility it provides, in the qualitative or 
quantitative realism and consistency of the results it yields, and in its 
possible advantages over other methods . 

Use of the directional-mean-free-path approximation will be seen to 
leave the equations on the directional level in a form amenable to treat- 
ment by vector generalizations of Lagrange 1 s expansion, certain useful 
forms of which are developed in Chapter V. A perturbation-expansion 
scheme based on those generalizations, also developed in Chapter V, is 
then applied in Chapter VI to the equations for the directional-mean- free- 
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path method. As an illustration of the method, a description of the 
solution procedure for steady- flow shock-wave structure, according to 
a simplified first-order system of equations, is given in Chapter VII. 
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CHAPTER II 


DISCUSSION OF MEAN-FREE-PATH CONCEPT 
2.1 Definitio n of Mean Free Path 

For proper understanding of the developments to follow, it is impor- 
tant to define precisely the mean free path. For a specified class of 
molecules. Maxwell’s mean free path is defined as the average distance 
travelled by one molecule of that class between successive collisions 
with all other molecules. (Refer, e.g. , to Jeans, 195^, 1952; Chapman 
and Cowling, 1961; or Vincenti and Kruger, 1965 .) According to this 
definition, if £ rel is the magnitude of molecular velocity measured 
relative to some specified coordinate system, Maxwell’s mean free path 
is formulated for the class of molecules under consideration as 


"class 


^r el/* average for class 
0 


class 


(2.1a) 


where 


average distance travelled per free path 
terminated by a molecule of the class 


(2.1b) 


<^rel5 


av for class 


average distance travelled per unit time 

by a molecule of the class (measured in 

some coordinate system) (2.2a) 


and where 


average velocity magnitude for molecules of 

the class , measured relative to some 

specified coordinate system (2.2b) 


class 


average number of free paths terminated per 

unit time by each molecule of the class 

within a very small volume (say dV^.) 

about the point under consideration (2.3a) 


l 
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= average number of collisions per unit time 

experienced by each molecule of the class (2.3b) 

"“probable number of molecules of the class" 
undergoing collisions per unit time in a 
Lsmall volume about the po int J (2 3c) 

[ probable number of molecules of the class] 
in the small volume about the point J 

which is denoted as the "average collision frequency per molecule of the 
class . If * 

In Chapter III, it will be convenient to consider various classes of 

molecules, to be defined. The Collision frequencies, 0 _ , are 

class 

calculated by considering velocities of molecules relative to another 
molecule, so concern over a specified reference frame does not enter 


*Note that for the particular class consisting of all molecules within a 
small differential volume d^^ about a point, the number of molecules 
undergoing collisions per unit time is twice the number of collisions per 
unit time, since each collision involves two molecules; thus, for this 
class , 

number of collisions per unit time in 

Q — 2 ^ — . ■ ■ " i — — — ■ — ... - ■■ ■■ 

number of molecules in d^ 

However, the collision frequency of interest in eq. (2.1) is that denoted 
by 0 and calculated by (2.3a), (2.3b ) 3 or (2.3c), and not just the 
second factor in the above equation. For alternate, but similar, 
discussion see Yincenti and Kruger (1965), pp. ^8-5 ^- The appropriate 
collision frequencies will be calculated as needed, for each class 
considered, in a manner that does not involve concern over the factor 
of 2, namely by use of eqs . (2.3). 
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into the calculation of ® c q ass » that is, the collision frequency per 
molecule of the class is not affected by the mean velocity of mass 
motion (e.g., see Jeans, 195*1 » p. 35 ). However, as noted in the defini- 
tions (2.2), each mean free path to be calculated will need to have the 
reference frame specified in which ? re q is measured. The distance 
travelled in unit time, or the velocity, measured in one reference frame 
is different from that measured relative to a reference frame moving with 
a different velocity. This fact is important to the discussions in the 
remainder of this chapter; and these discussions are important to 
understanding the method developed in Chapter IV. 

2.2 Definiti_ons_for Calculation of Mean Free Paths 

2.2.1 Configuration space, velocity space, phase space , 
and velocity distribution function 

Consider the physical configuration space in which molecules of a 
gas are moving and interacting with each other, and let the reference 

frame of the observer be a system of Cartesian coordinates denoted by 

. -y -y -y 

x , x , and x , with the orthogonal unit vectors e , e , and e , 

respectively in the directions of the three space coordinates. A point 

in this space is defined by the position vector 

r = ex+ex+ex (2.4) 

1 1 2 2 3 3 

Denote by 

d °j/ = dx dx dx (2.5) 

¥ r 12 3 

a differential element of volume in the physical configuration space, and 

denote the range of all values of r within a certain dY r by 

(r, &y ) * which can also be thought of as the volume av about a 
—y 

certain point r. 
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Denote the velocity of a molecule, measured in this frame of 


reference, by 


->■ ->■ -y ~y 

E, = eE,+eE,+eZ 

1 1 2 2 3 3 


( 2 . 6 ) 


where the E,. are the Cartesian components of the molecular velocity. 
Consider molecular- velocity space (in which t is a "position vector") 
and denote a "differential volume of velocity space" by 


i¥ r 


s d{ d; as 

1 2 3 


( 2 . 1 ) 


Then denote by (£, ) the range of all values of E, within a certain 

d °(/ or, equivalently, the volume d^ about a certain point Z in 
velocity space. 

It is customary to speak of the "phase space" of a molecule as the 

six-dimensional hyper-space having coordinates x,x,x, £,£,£. 

12 3 12 3 

A point in the molecular phase space can be defined by the six-dimensional 
"position vector" x with components x , x^ , x^ , E, ^ , Z^ • A differen- 
tial element of "volume" in this space about the point x can be denoted 
by (x, d^ ), where 


d y - dx dx dx d£ d£ d£ 

X 1-23123 


= d 




( 2 . 8 ) 


If the total probable number of molecules in (x, d^^.), (i.e., 

located in (r, d ^) and having molecular velocity in the range (5,d^)), 
at time t is denoted by dN|, then the probable number density at point 
x(r ,1) in phase space at time t, called the velocity distribution 
function and denoted by f(r,£,t) , is defined by 



If the total number of molecules in (r, d %) at time t is 
denoted by 

dN = / d N (2.10a) 

a 111 | 

then the number density at point r in configuration space is 

n(r,t) = ^ f(r,f,t) d'V (2.10b) 

r r all? 9 


2.2.2 Collision frequency of a molecule 


Denote by 


a( v? 


the probable total number of molecules in 


(r, d r r ). (I, d^ ) that experience collisions with other molecules in 

unit time. Denote by | the velocity of each other molecule undergoing 

a -* 

a collision with a molecule of velocity £. Let g be the relative 
velocity before collision: 



(2.11a) 


and denote its magnitude by g: 

g = |g| = |f a 4| (2.11b) 

-> -> 

For each £ * consider the plane that contains the center of the £ 
a 

molecule and that is perpendicular to the relative-velocity vector g. 

Let b be the distance from the center of the £ molecule in this 
plane, and let e be an angle measured in this plane, centered at the 
| molecule. (Refer to Chapman and Cowling, 1961, P- 6l, fig. 6. To 
aid understanding of the impact parameters b and e , see also Chapman 
and Cowling, 1961 , p. 57 » fig- 3; and Vincent i and Kruger, 1965 : 
fig. 5 in Chap. II (p. 37) and fig. 3 in Chap. IX (p. 35l)-) Then, as 
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is shown in many books on kinetic theory of gases. 


a % d % i* r i zn J ff a gbdbdiT| (2.12a) 

? all? ^aLo o J 

a 

where a is an effective "molecular diameter" (which in reality is 

infinite but which may be "truncated" at an appropriate value beyond 

which the effects of molecular interaction are negligible, depending on 

the nature of the intermolecular force potential; cf. Chapman and 

Cowling, 196l), and where f = f(r,^,t) and f = f(r,"g ,t). Then, 

a a 

since f(r,^,t) is independent of , equation (2.12a) can be written 

a 

as 


d(N ch 


0 + dN^ 

5 ? 


(2.12b) 


where dN is given by equation (2.9) and where 

I 


0^ = 0^(r,?,t) 

K ? 


d ( N ) ->■ 
c 5 
dN^ 


/ * 
^ic a 



gbdbde 


(2.12c) 


is the probable fraction of all the dN^. molecules in (r.d^.), (t,d^) 

undergoing collisions per unit time. Thus 0^ is the collision frequency 

per molecule of the class of cLDL^, as define! by equation (2.3c) 

£ 

The total number of molecules in (r, d ^r ) undergoing collisions 
per unit time, denoted by dN , is found by integrating equation (2.12a) 
over all £ : 


dN = / 0 r dN+ = df j 

c aLl? 4 ? all? 


% 


I * 

all ? 


r r 2 n r a 

/ / f 

L o o 


^gbdbdej 


(2.13a) 
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This may also he written as 


ail = 0 dW (2.13b) 

C 

where dN is the number in (r, dV) (eq. 2.10)) and where, for the class 
of dN (all molecules in (r,d^)), 0 is, according to equation (2.3c), 

and also from (2.13a) and (2.13h), 


0 = 0(r,t) = -r^ 2 - = — / 0^(r,t,t) f(r,|,t) aV (2.13c) 

n ail? 5 K 

This is the probable fraction of all the dN molecules undergoing 
collisions per unit time. The collision frequencies 0^. and 0 depend 
on the nature of the interaction potential and on the particular appropri- 
ate corresponding specification of a. 


2.2.3 Average_ values of mo lecul ar properties; mean mass 
v elocity and pe c ul ia r (random) velocity 


Each molecule may be said to have certain properties that may depend 
on its position, velocity, and time, and are denoted generally by 

— y 

<f> = <f>(r,5,t). Examples are: mass/unit mass = 1; momentum /unit mass = 

1 > > 

energy/unit mass = — ; ££ ; 0^. ; etc. (These examples for <J) can 

all be considered to be "specific" values of extensive molecular proper- 
tiesO For any such molecular property, <|>(r,g,t), one may define its 
average value, for the class of all molecules in (r,d^ ) (regardless 
of molecular velocity) , as 


<f> av (r,t) = <4>(r,|,t)]> 


J ^ +(r,t,t) dN^ 

allg K 

/ dW 
all| t 


= r 1 • " / + <Kr,|,t) f(r,t,t) &V 

n(r,t) a 11C 


(2.14) 
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Particular examples of average values that are of concern in this 
chapter are found as follows: First, the mean mass velocity of the gas 

is defined by 

v =<yy (2.15) 

Then the molecular velocity relative to the mean mass motion, denoted as 
the "random molecular velocity" or as the "peculiar velocity" is defined 
hy 


C 



V 


(2.l6a) 


= eC+eC+eC 
1 1 2 2 3 3 

where the Ch are the Cartesian components of (5. The magnitudes of 
the vectors % and $ are, respectively, 


5 = 




C 




(2.17) 


Quantities that will be of interest below are the average values of £ 
and C, found by taking = £ and = C in equation (2.lU) to obtain, 
respectively, and . Note also that taking <f>(r", £,t) = 0+ 

(the probability of a molecule in the class 

in unit time) one finds (cf. eq. (2.13c)): 


dN->- undergoing a collision 



(2.19) 
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2 . 3 Mean Free Path for a_Gas at R est ; and Mean Free Path 
Relative_ t_o Me an M o tion of the Gas , X 


Maxwell’s mean free path was originally calculated for a gas at 
rest, that is, with no overall mass motion, and in a uniform state (see 
derivations outlined by Jeans, 1952 » and by Chapman and Cowling, 1961). 
In that case, C = £ , and the average velocity magnitude measured 
relative to the reference frame of the observer (i.e., relative to the 
coordinates of r space; see § 2 . 2 . 1 ) is then <^C^> . Then Maxwell’s 
mean free path for the class of all molecules in a small volume about a 
point is, from (2.1a), 


^(V= 0 ) 


0 


(2.20a) 


For a gas composed of rigid elastic spheres, in equilibrium, and with no 
mean mass motion (see above references). 


0 = /2 n n a 2 <^C^ 

so in that case (2.20a) becomes 


i ^ 

($■0) 


1 

/2 n n a 2 


(2.20b) 


(2.20c) 


The notion of mean free path has been extremely useful in the exten- 
sion to gases with nonvanishing mean mass motion (V ^ 0 ) for the 
calculation of transport phenomena (e.g. , see especially Chapman and 
Cowling, 1961 ; Present, 1958 ; Lighthill, 1956 ; Hirschf elder , Curtiss, 
and Bird, 1964 ; Guggenheim, i960; and Vincenti and Kruger, 1965). 
However, the mean free path used in that extension has been calculated 
from equations (2.20), or 
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That is, the relative velocity magnitude in equations (2.1) and (2.2), or 
the average distance travelled hy a molecule in unit time, has "been 
measured in a reference frame that moves with the mean motion of the gas. 
One may therefore refer to A, defined by equation (2.2l), as the mean 
free path relative to the mean mass motion . As noted above (at the end 
of § 2.1), although 0 is not affected by the mean velocity of mass motion 
V, the calculation of the average velocity magnitude, 

depend on the reference frame in which it is measured. This is further 
discussed in § 2.5* 


X 


rel/av 


does 


2.b Absolute Mean Free Path, A* 

It will be found useful to consider what may be denoted as the absolute 
mean free path , defined by 

X* E ^ (2.22) 

From the definitions (2.1) and (2.2), for the class of all molecules in 
(?,dV r ), since £ is the distance per unit time travelled by a molecule 
relative to the observer’s reference frame (i.e., relative to the coordin- 
ates of the r space defined in § 2.2.1), Af^ is the average distance 
travelled , relative to the observer T s reference frame , by a. molecule 
between successive collisions . The significance of A*, and its comparison 
with A, are discussed in § 2.5* 

2 . 5 Comparison of A* with A; and Significance of A* 

In comparing A*, defined by equation (2.22), with A, defined by 
equation (2.21), it will first be noted that if V = 0, A* is the same as 
A. However, at least for an equilibrium distribution with V ^ 0, A* > A. 
For example, consider a uniform flow of a gas in translational equilibrium 
with very large V. Then within the mean collision time of a molecule. 



t = 1/0 , nearly all the molecules move a very large distance in the 
observer's reference frame. In fact A* 00 as V « . However, A 
does not depend on V. This illustration is made more precise by the 
calculation in the following paragraph. 

Consider calculation of A* and A for a gas having a Maxwellian 
distribution function 


f = f 

e 


n(e e /n ) 3/2 



(2.23a) 


and having mean mass velocity 



with magnitude V, where 


_ _m_ _ _JL_ _ Y/2 

e " 2kT 2RT a 2 


(2.23b) 


and where m is the mass of a molecule, k is Boltzmann's constant, T 
is the gas temperature, y is the ratio of specific heats of the gas 
(c /c ), and a = (YRT) 1 / 2 is the speed of sound in the gas. Without 
loss of generality, one can take V = e^u for this calculation, since the 
x axis can arbitrarily be aligned with the velocity vector. Then 

C 2 = (5 -u) 2 + s 2 + K 2 (2.23c) 

1 2 3 


It is readily found by substituting equations (2.23a-c) into equation 
(2.lU) with <j> = C that 



(2.23d) 


(cf. Chapman and Cowling, 1961 , p. 7*0 * and with 


(j> = £ that 


«> 



(1 + 2f3 u 2 ) 
e 


erf (/g - u) 
e 


(2.23e) 
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where, by definition, the "error function" is 


erf n 


(2//ff) / ' e" ? d? 


With definition of the Mach number 


(2.24) 


M = u/a , or /fT u = /y/s* M = M (2.25) 

and with use of the appropriate asymptotic expansions in equation (2,23e) 
for M -> 0 or M -* » , one finds: as /jT[~ u = M -* 0: 

= <C>(1 + |m 2 -^+^oM 6 -+ ) (2.23f) 

and as /jT~ u = M °° : 

e 


<*> 


~ u{l + i M 2 + n 1/>2 M 1 e M 


[- J M 4 + 0(M 6 )]} 


(2.23g) 


Therefore, from the definitions (2.21) and (2.22), for this case of the 
Maxwellian distribution. 


A 

A 


* 



1 + M 2 + 0(M 4 ) 


as M 0 


and 


\ 


n 


l A - 

■ M 


h 


-M 2 . 


(2.23h) 


(2.23i) 


Hence, for small M, A* % A ; and for large M, A*/A = 0(M). 

One indication of the significance of A* is the following: 
Suppose we wish to consider the appropriate dimensionless parameter 
characterizing a molecular flow. Let L be a characteristic reference 
length of the problem, 0^ be a reference collision frequency per 


l6 



molecule, and be an appropriate reference velocity. Since X* = <^£^/0 

is the actual mean free path relative to the observer, it is suggested 
that 


Kn* ^ 





L 0/ 


ref 


( 2 . 26 a) 


(a "Knudsen number based on the absolute mean free path") is the most 
appropriate characteristic dimensionless parameter, or 


Kn* = 



where 


c 


r 


<«> 


( 2 . 26 b) 


(2.27) 


The parameter c_^/L0^ Vincent i and Kruger T s ( 1965 , p. 380) parameter 

g or Liepmann, Narasimha, and Chahine’s ( 1962 , p. 1319) parameter l/a Q 
(the appropriate expansion parameter for the Chapman-Enskog procedure) . 

Those writers made appropriate specifications of c^ for low speed flows 
and for high Mach numbers. Since, from (2.23f,g), <^> ^ as M 0 

and ^ u as M -* °° , we see from (2.26) and (2.27) that use of A*, 

rather than A, has led to a result for c^ (eq. (2.27)) that a u t omat i c al ly 
includes the separate specifications of c^ for low-speed and high Mach 
number flows that were made by both Liepmann et al. ( 1962 ) and Vincent i 
and Kruger ( 1965 ). 

The simplified calculations of the transport coefficients, from the 
elementary kinetic theory, treat molecules as carrying with them values 
of various properties (e.g., momentum) that are characteristic of the 
location of their last collision, a mean free path away (see references 
mentioned above, prior to equation (2.21)). Although A* is the actual 
mean free path relative to the observer, use of A has been convenient, 
and its use is justified especially by the fact that the actual average 
free path of the molecules moving normal to the flow (used, e.g. , in 
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calculating shear rate) is close to X, However , for flows in which V 
is large and/or changes rapidly, use of X* (in a certain way) is suggested 
to he perhaps more appropriate than X, and, in fact, use of a directional 
absolute mean free path (to he introduced in Chapter III) may he most 
appropriate in cases of significant translational nonequilibrium, where 
the directional mean free path of the molecules varies significantly with 
the direction of the molecular velocities. 
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CHAPTER III 


THE DIRECTIONAL LEVEL OF DESCRIPTION 
AND ITS RELATION TO THE KINETIC AND MACROSCOPIC LEVELS 


3-1 Introductory Remarks 


The purpose of this chapter is to introduce the basic concepts , 
definitions, and formal equations that constitute the directional level 
of description. This formal development is to form the basis of a method 
for treating gas flows in translational nonequilibrium (introduced in 
Chapter IV). Because of the suggested importance of the directional 
aspects of significant translational nonequilibrium, it is proposed to 
treat each molecular- velocity direction separately, which is possible on 
the directional level. The role of the directional level of description 
in its relation to the kinetic and macroscopic levels is emphasized in 
the formulation. 

It is convenient first to discuss the Boltzmann equation (§3.2), 
which involves both the collision frequency per molecule of the class 
dN-+ (i.e., all molecules in (r,d^£) with velocity in (|\d%.)* cf. 
eq. (2.9)) and the inverse collision frequency per molecule corresponding 
to that class. Then introduction of spherical coordinates in molecular- 
velocity space (§3.3) facilitates the subsequent formulation. 

In Chapter II, different classes of molecules have already been 
mentioned. In particular were discussed: the class of all molecules in 

in 

(r,dV r ) regardless of molecular velocity. In this chapter (in § 3.^-) it 
is convenient to distinguish, and define precisely, four classes of mole- 
cules. For each of these classes, one can define certain properties and 
write equations of change for those properties. The Boltzmann equation 
applies on the kinetic level of description. One then can define the 
intermedi ate levels of description, involving marginal distributions of 


(r,d Y r ) with velocity in (£,d*^), and the class of all molecules 


19 



molecular properties. Finally, one finds the full moment level, for which 
Enskog's general equation of change governs the various macroscopic proper- 
ties, which are molecular properties integrated over (or averaged over) 
all molecular velocities. The directional level of description comprises 
an intermediate class, and its consideration will be found to be useful. 

In this connection, directional properties of the molecular motion are 
defined, including a mean free path that can be different for each 
molecular-velocity direction. Molecules with velocity in a certain direc- 
tion can be thought of as constituting a "species" in a mixture of species 
of all directions, with each direction having its own mass density, momentum 
density, energy density, collision frequency per molecule, mean free path, 
etc. Equations of change can then be written for each of the "directional 
properties" (§ 3-12). 


3.2 The Boltzmann Equation 


The Boltzmann equation is an equation of change for the probable 
number of molecules in (r,d^),(f,d^) at time t (eq. (2.9)), 


dN+ = f(r,|,t) dT d Y 

K * 


(3.1) 


or for the probable number density, f, in phase space. Many derivations 
of the Boltzmann equation have been given (cf. Grad, 1958; Chapman and 
Cowling, 1961, pp. 46, 47, 63-66; or Hirschf elder , Curtiss, and Bird, 
1964, pp. 444-452). The equation of change for dN^. in phase space 
involves a "net source term" owing to collisions. In writing the net 
source term one needs an expression for the number of molecules in 
(r ,d?£) ,(|*,d^") undergoing collisions per unit time, and hence the number 
lost from that class per unit time owing to collisions (cf. eqs . (2.12)): 


d(N )+ = 0+ d!4 

c £ £ ? 


(3.2) 


where 0.* 
6 


defined by (2.12c), is the fraction of the number dE* 

K 
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undergoing collisions per unit time, and hence the fraction lost per unit 
time, and is denoted as the average collision frequency for a molecule of 
the class of dN^. . One also needs an expression for the number of mole- 
cules undergoing "inverse collisions" per unit time that end up in (r,d*>K) 5 
Cc,d^) after the collisions, or the probable number gained by that class 
per unit time owing to collisions: 


— 2 H o' ~ 

d(N')-*. = d % d %> f » M / f 'f ' g'b' db' de' (3.3a) 

C ? r * all? ' 5 a L o o a 


- ^ *%[ I** i'iw*) 

a 


f g b db de (3.3b) 
a 


] 


©4. dN-y 

5 5 


(3.3c) 


where 


E / d % 


a 


all? a L o o \ a 


/ 2 vp) 

o ova 7 


f g b db de 
a 


(3.3d) 


is defined by equation (3.3c) as the fraction of the number dN+ that are 
put into (r,d?£) ,(!,d^ ) per unit time by inverse collisions. This 


fraction , 


Q k 


may be denoted as the "inverse collision frequency per 


molecule" for the class dN-* . In equations (3-3), 


f f = f(?,f',t) , f' = t) 

a a 

and 

s’ = Ib'I = if; - V i 


\ (3.4) 


where E ! and E 1 are the velocities of two molecules before an "inverse 
a ^ 

collision" or after a" direct collision." The velocities and are 
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not independent, but are related through the impact parameters and inter- 

4 4- , . 

molecular potential to £ and £ . In equation (3- 3b) the well-known 

a 

relations have been used: 

d y, ; g' = g ; b' = b, e' = e (3.5) 


In terms of dN_^ , the Boltzmann equation is written as 


F 

(dttj + ? r * (f dN^) + ^ dNj = d(r - d(N c ) 


e (0* - 0 ) dN 

-4 -4 y -4 

5 5 5 


(3.6a) 

(3.6b) 


or, upon dividing by a^ 


|£ + v 

8t r 


(If) + V - ( -£■ f) 

5 m 


(04- - 0+) f 

6 K 


(3.7) 


where 


8 

1 

8 

4 

8 



1 9x 1 

+ S 2 

9x 2 

+ 6 3 

8x 3 



8 

i ^ 

4 

+ e 

2 

3 

9? 2 

4 

+ e 

3 

8 

9? 3 


> (3.8) 


4 

and F is the external body force acting on a molecule. 
Jd 


3. 3 Distribution Function in Spherical Coordinates in 
Molecular-Velocity Space 

For subsequent developments it is convenient to transform the 
coordinates of velocity space to spherical coordinates. 

With the molecular velocity £ defined by equation (2.6) and its 
magnitude by (2.17) » the transformation is 


22 



s 

5 


1 

2 
3 


£ COS <p 

5 sin <p cos 0 
£ sin <jp sin 0 


l (3.9) 


where, for a molecule located at point r and having velocity £ (see 

figure 3.1):, <P is the angle between the axis (which is parallel to 

the x^ axis) and the molecular-velocity vector X ; and 0 is the 

angle of the projection of the vector % on the plane of £ , £ (which 

2 3 

is parallel to the x^ , x 3 plane). The angle 0 is measured from the 
axis, as shown on figure 3.1. 

The unit vector in the direction of % is then defined by this 
transformation (3*9) and equation (2.6) as 


e+ = = e cos <p + e sin <p cos 0 + e„ sin cp sin 0 (3.10a) 

£ t 1 2 3 

Noting that e*. is a function only of 0 and <p , we will also find it 
convenient to use the notation 


0<p 




-> 


( 3 . 10t ) 


The transformation of the velocity-space volume element, 
given by 



is 


ay { = a 5i a « 2 a { 3 


9 ( *? 2 » C 3 ) 


9(S, 6, <p) 


d? d0 dcj£> 


(3.11a) 


= 5 2 sin <p d£ d0 dc p 


(3.11b) 


8(5 1S 5 2 ,S 3 ) 

where ttt — s — is the Jacobian determinant of the transformation. 

0, <p) 

(From Fig. 3.1, the three edges of the volume element are £d(p , £ sin<pd0, 
and d£ which when multiplied together also give as eq . ( 3 . lib ) . ) 
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It will be found convenient to also consider the solid angles in 
molecular-velocity space (see fig. 3 . 2 ): 

dft 0 ^ = sin q? dtp de (3.12) 

e=n 

aa = / da = 211 sin 9 d<p (3.13) 

T e=-n v 

in terms of which equations ( 3 . 1 l) may also be written as 



S 2 


d£ da rt ^ 

e<P 


(3.lUa) 



(3.lVb) 


(The significance of the subscripts E, , Q$ , and <p is discussed in 
§ 3.b.) 

The velocity distribution function defined by equation ( 2 . 9 ) may be 
written with 0 , (p as arguments. Thus define 

f(r,|,t) = f*(r,C, 0 ,<p,t) ( 3 . 15 ) 


From this point on, however, the superscript * in equation (3.15) will 
be dropped, with the understanding that f (r, g ,0 ,<^,t ) is actually f* 
defined in equation (3.15). Thus from equations (2.9)* (3.1l), (3.12), 
(3.13)* and (3.1*0 ve now have 



f(r,S,e,<p,t) dj/ 

f S 2 dfi 0<p 


( 3 . 16 a) 


- f « 2 (§) ia <f 


( 3 . 16 b) 


26 


f E, 2 dC d0 sin <p d<p 


( 3 . 16 c) 



3 . 4 Classes of. Molec ules^ Considered; Total Number 
and Notation for Each Class 


The remainder of this chapter involves consideration of four classes 
of molecules: 

(a) The class of all molecules in (r, d % ) and having velocity in 
(t,d^) contains the number 

dN+ = f (r (3.1Ta) 

where the notation ) has been defined in § 2.2.1. This 

class is denoted as the class dN|- . The subscript % is to be used 
generally to denote quantities in this class, which can be functions of 
X (or of 5 , 0 , *). as well as of r and t . 

(b) The class of all molecules in (r,d^) and having velocity 

direction in ( 0 , dfi__) (regardless of the magnitude 5 ) contains 

T 0*P 

the total probable number 

£*=oo oo 

dN e<p E / = d % dn ea> / f 5 2 (3.171) 

r 5=0 ^ *0 

where ( 0 , <D , dft_) indicates the range of 0 and a> in the solid angle 

Q( r + -> 

d^ 0 ^ about the vector = e ~£ • (See fig. 3-2.) This class is 

denoted as the class • The subscript 0<p is to be used generally 

to denote quantities corresponding to this class which may be functions 
of 0 and <p as well as of r and t (but not functions of the velocity 
magnitude 5 ) . 

(c) The class of all molecules in (r,d^) and having velocity 
angles in (<p,d^) (regardless of the values of 5 and 0) contains 
the total probable number 

0=n 5 =°° . n 

dN E / / dN t = d Y an — / d0 / f 5 2 d5 (3.17c) 

e=-n 5= O 5 v -n O 
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where (<p , dQ ^ ) indicates a direction in the range d<p about the angle 
C p anywhere in the solid angle dQ^ (see fig. 3.2). This class is 
denoted as the class dN^. The subscript (p is to be used generally to 

denote quantities corresponding to this class, which may be functions of 
<p, as well as of r and t (but not functions of £ and e). 

(d) The class of all molecules in (r,d^) (regardless of g - ) 
contains the total probable number (cf. eqs . (2.10a,b)) 

cp=Tl 0=11 £=°o n n 00 

dN = / J / dlL*. = dfY J d$?(sin<p) / d0 / f d? (3.17d) 

<p=o e=-n 5=0 5 r o -n o 

This class is denoted as the class dN . Quantities corresponding to this 
class may be functions of r and t (but not of £ , 0, or <p ) . 


3.5 Total Number of Molecules Undergoing Collisions per Unit Time in 
Each Class of Molecules 


For the class dN.^ , the total number of molecules undergoing 
collisions per unit time is given by equation (2.12b) as 


d(N U 
c 5 


= 0_> dN. 


&V r djS 0+ f 


( 3 . 18 a) 


where 0_*. , defined by equation (2.12c), is the fraction of the number 
dN.* that undergo collisions per unit time (cf. eq. (2.3c)). 


For the class dN _ , 

e<p ’ 

collisions per unit time is 


the total number of molecules undergoing 


d(N 


c 6(p 


/ 

5=o 


d(N 


= dV 


dfi 


0<p 


| 0> f ( 2 d? 


0 ecp dN e<p 


(3.18b) 


which therefore defines 


0<p 


as the fraction of the number 


undergo collisions in unit time (cf. eq. (2.3c)). 


diL that 

0q> 
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For the class*' dN^p, the total number of molecules undergoing 
collisions per unit time is 

0=n n °° 

d(H ) = / / d(N U = ^ / a 0 / 0 , f dc 

c <P e=_n 5=0 ? ^ -no? 

= 0^ dN^ ( 3 . 18 c) 

which therefore defines 0 as the fraction of the numher dlU> that 

<P T 

undergo collisions in unit time (cf. eq. (2.3c)). 

For the class dN, the total number of molecules undergoing 
collisions per unit time is 

<p=n 6 =n e ,=°° ' n n °° 

dN = J / J d(N )-► = / d<p( sin ) / d0 / 0 4 f d£ 

c <p= o e=-n 5=0 G ? r o _n o ? 


= 0 dN 


(3.l8d) 


which defines 0 as the fraction of the number dN that undergo 
collisions in unit time (cf. eqs . (2.3c) and (2.13b)). 


3 • 6 Total Number of Molec ules _Undergo ing_ "Inverse Collisions " Per 
Unit Time Co rresponding to Each Cilass of Molecules 

For the class dN-> , the total number of molecules undergoing 
"inverse collisions" (molecules that end up in the class dN|. after the 
collision) per unit time is given by equation (3.3c) as 

d(N')-*- = 04- dm = &¥ 04- f (3.19a) 

c 5 5 5 r ? 5 

where 04 . , defined by equation (3.3d) 3 is the fraction of the number 
&KU. that are put into the class dm per unit time by collisions. 

5 ? 
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For the class dN , the total member of molecules undergoing 
inverse collisions that end up in (r,d^) with molecular velocity in 
the direction range ( 0 , <p, dJ2 Q ^) per unit time is 


d <”c>e4> = I d( "c ) i " a % d V / 6 f f 5 2 d « 

T E=o t o 9 


= 0 ' dN 

6 <p 0<p 


(3.19b) 


which defines 0 ' as the fraction of the number dN that are put into 

— Of 0f> 

the class dN_ per unit time owing to collisions. 

0 (p 

For the class dN^ , the total number of molecules undergoing inverse 
collisions that end up in (r,dY) with molecular velocity in the direc- 
tion range ((pi dfyp) per unit time is 


d(w ) 

C <f 


0=n 5=°° . n °° 

/ / d(N* U = d Y an / d e / 04 f ? 2 d? 

e=-n 5=0 ? -n o K 


0 ' dN_ 

<p <p 


:3.19c) 


which defines 0 ! as the fraction of the number dN that are put into 

<p <p 

the class dN^ per unit time owing to collisions. 

For the class dN, the total number of molecules undergoing inverse 
collisions in (r,d^) per unit time is 


<p=n o=n g=~ 

dN T E / / / d(r)_ 

C J J J C i 

<p=o 0 =-n 5=0 


n n 

d^f J d<p (sin <p> / do / 04 . f c 2 

r o -no? 


0' dN 


(3.19a) 


which defines 0 ' as the fraction of the number dN that undergo collisions 
per unit time, and is therefore equal to 0, a well-known result (see 
§ 3.12 below) . 
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3*7 Molecular-Property-Density Distributions for Each Class 
of Molecules 


3*7*1 Definition of property-density distributions 
fo r each class 

For any molecular property <f>(r 9 ^,t) (cf. § 2.2.3 above and § 3.7*2 
below) one can define expressions for the density distributions of the 
property <j> for each class. 

For the class dN^ , the density of the molecular property <j> per 
unit volume of configuration space and per unit volume of velocity space 
is defined simply as 


- - ♦ 

V r,5>t) s = * f 


(3.20a) 


For the class dN^ , the density of the molecular property <j) 
per unit volume of configuration space at r and per unit solid angle 
of velocity space at 0 9 cp is defined as 

/_ * ^ 

* Q<p (?,e,<P,t) s k W ~^ l — = / <P f c 2 d5 (3.20b) 

r T 0<p O 


For the class dN<p , the density of the molecular property <p per 

~y 

unit volume of configuration space at r and per unit solid angle of 
velocity space at angle <p (integrated over all 0) is defined as 


0=11 £=~ 

/ / <j> dN+ n 

5 — 1 = iff I de 1 * f E 2 d 5 

r r ? -II o 


(3.20c) 


31 



For the class dU, the density of the molecular property <f> per 
unit volume of configuration space at r is defined as 


$>(r,t) 


=n e=n £=<*> 

I I I <p 

< P=Q Q=— II g=Q _ ^ 

a y 

'r 


n n °° 

= / d<p (sin <p ) / d 0 / ()> f ? 2 dC (3.20d) 

o -no 


3.7*2 Specific examples of molecular properties of interest, 
and their dens ity_ distributions in each class 


Examples of molecular properties that may be of particular interest 


are : 


4) 


(£) _ 


, , *>■ — ->■ p. 

= 555...E = e e . e . . . e . 5 , £=0,1,2,, 


£ times 


5 ? C 


£ times 

(which is a tensor of order £) , 


(3.21a) 


( 0 ) 

(j) = m (()_>_ , £=0,1,2,... 

£ 


( 3 • 21h ) 


1 f £ ) 

<t> = 4 I : m *+ £=0,1,2,... t 

2 K 


(3.21c) 


For any two Cartesian tensors of second order or higher, say 


3 3 


3 3 


A - I 1 A. . e.e. and B = 7 £ e,e B 

M A U 1 J ,1, ofln k £ 


i=l j=l 


k=l £=1 


k£ 


the "double-dot" product is defined by 


3 3 

A * B = l l A,, B„ 


i=i j=i 


1J 


where the components A— 311(1 may themselves be tensors of zeroth 

1 J KX, 
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(where 


I i Vj 


+ + SO 6 -, 

2 2 3 3 


is the unit tensor, or "idemtensor" ) 


<f> = , £=0,1,2,... 


<p = f(r,l,t) 


<j> = Xi£ (to be defined below in i 3.11) 


<f> = ♦ 5 


4> - <P 0->. 

K 

4> = 4> ei 

5 


(3.21d) 

(3.21e) 

(3.21f) 

( 3 . 21g ) 

(3.21h) 

(3.21i) 


where <f> may* he any other molecular property. Some of these examples are 

discussed and the property densities for each class indicated below. 

The number density in each class is found by taking <j> from equation 


(3.21a) with £ = 0 and using equations (3.20). Thus, with (j> = c^° = 1* 


/■v ->■ . 

n+(r,£,t) 


dN^ 


= f(r,|,t) = 




(3.22a) 


n e(p ( r ,0,<p, t) 


/ f € 2 d? 
o 


dN 


0<J> 


d V 

r r 0Cf> 


(3.22b) 


order or higher, i.e., scalars, vectors, or higher-order tensors. The 


orders of the component tensors 


and B 

orders lower than the orders of the tensors 
if A and 


A 


are scalars. 


B are second-order tensors, then the components 


are by definition two 
and B respectively; 

and 
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(3.22c) 


y?,<p,t) = ^ / ae / f 5 2 a? = 


m. 


2n 


-n 


&Y r an 9 


n(r,t) = / df(sin<p) / d0 / f E, 2 d£ = ^7- 


n 

f 

-n o 


(3.22a) 


In particular, may "be denoted as the "directional number density" 

(which is the same as n^ if f does not vary with 0). 

The mass density in each class is found by taking $ from equation 
(3.21b) with £ = 0 and using equations (3. 20), or by simply multiplying 
each of equations (3.22) by the molecular mass, m, to obtain, respectively. 


= mn^ , p 0<p = mn 0<p , p f = mn^ , P = mn 


(3.23) 


The quantity p may be denoted as the "directional mass density." 

0<P 

If we define = ^class vit]rl ^ g iven "b y equation (3.21b), 

where the property densities ^ c q ass are ^ e ^i ne(3 - by equations (3-20) , then 

*(0) 

P class 'class 


The momentum density (or mass flux) in each class is found by taking 
$ from equation (3- 21b) with £=1 (cj)=m5) and using equations (3.20) to 
obtain 


^ = mff = e_^ mgf 

y x OO OO 

$2 = / mffg 2 d£ = m e / f 5 3 d£ 

0(0 J s s s 0<P J ^ * 

T o o 


/ \ n n <» n °° 

= 2? / de J = If I d0 (O / f C 3 d ? 

-no -no 


<p 


pV = = / d<p(sin <p) f dd f mff£ 2 d E, = m J d<?(sin q>) j d0(e )/ f£ 3 d£ 


n 


n 


n 


(3.24a) 

(3.24b) 

(3.24c) 


-n o 


-n 


e<p 


(3.24d) 
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The quantity could he denoted as the "directional momentum density" 

(or as the "directional flux of mass"). 

The energy density in each class is found hy taking <J> from equation 
(3.21c) with Z=2 ($ = ^ m^ 2 ) and using equations (3.20) to obtain: 


i|: *<*>. i 

n ■ • — o 


m f 


(3.25a) 


— | • = 

2 *• 6<p 


CO oo 

= / | m f ^ dg = \ m / f 5 4 d? 


(3.25b) 


| I: 


/ \ n 00 n °° 

i = / d0 / | mfS 4 d£ = | m / de / f £ 4 d? (3.25c) 


<P 


211 2 
-n o 


-n o 


1 ( 2 ) 1L 11 w 1 1 il n °° 

o’ |: = J df(sincp) / d0 / — mfC 4 d? = — m / d<p(sin q?) / d0 / f£ 4 d£ 


-II o 


-n o 


(3.25d) 


The quantity calculated in equation (3.25b) may he denoted as the "direc- 
tional energy density." 

The quantity 


( 2 ) 

V f m ft f dS = e e m / f C 4 d 5 

o T ^ o 


is the "directional flux of momentum." The quantity 


. q . OO OO 

\ 1 : V E ^ I m 5 2 f f 5 2 dC = e \ m / f S 5 d? 
o T o 


(3.26) 


(3.27) 


is the directional flux of energy." 

By taking <p from equation (3.21g) and using equations (3. 20) one 
obtains for each class the dens ity of the molecular property <p or the 
flux of the molecular property <f> . For example, the quantity 
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(3.28) 


= [ m 4 &) f 52 d5 

( £-l) 

is the directional flux of the molecular property m <f>_^ as defined 

in equation (3.21b) with (3.21a) (the directional density of the latter 

property being ^ ). 

By taking <j> from equations (3.21h) or (3.21i) and using equations 
(3.20) one obtains for each class the density of the property cj) respec- 
tively lost or gained per unit time by each class owing to collisions. 

For example, with (J>= 1 (c^=0|-), equations (3.20) give the number density 
of molecules undergoing collisions per unit time for each class 
(respectively: <1(^/1^ d(N c ) 0<p /d% dfi 0<p ; dtNj^/d^ d£^; and 

dN c /d^/; cf. eqs . ( 3 . 18 )). 


3.8 Average Values of Molecular_Properties for Each Class 
of Molecules 


For any molecular property (f>(r,^,t), an average value of the 
property possessed by the molecules of each class is defined by 


^^^lass 


/ * » 

class ^ _ class 

f dN-> n class 

class 5 


(3.29) 


where the $ are defined by equations (3.20) for each <j> and the 

C _L_ cLS S 

n are defined by equations (3.22). In § 2.23, only the averages for 

eras s 

the class dN (all molecules in (r, d Yu thus averaged over all t) 
were considered. However, the definition (3.29) can apply to any of the 
four classes defined in § 3.^. 

The subscripts g, 0tp, Cf are used to indicate what part of the 
molecular velocity the resulting quantity can be a function of. 
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The average value of 
therefore 


for the molecules in the class is 

K 




5 


(3.30a) 


The average value of $ for the molecules in the class dN e< p is 


/ <f> f s 2 a? 

<*>e<P 5 *4 

I t £ 2 d{ 


n 


9jp 

09 


(3.30b) 


The average value of <J> for the molecules in the class 

n 00 

/ ae / <f> f £ 2 

-n o _ 

n n- 

/ de / f £ 2 d? 9 

-n o 



is 


(3.30c) 


The average value of cf> for the molecules in the class dN (cf. eq. 
(2.14)) is 

n n °° 

/ dtp (sin <p) J d 6 / <|> f £ 2 d£ 

<?> 5 ° n T 4 ' f (3. 30d) 

/ dcp ( sin 9 ) / d 0 / f 5 2 d£ 

o _no 

Note that in equations (3.30) the quantity $ c q ass corresponding to any 

molecular property is defined by equations ( 3 - 20 ). 

Particular average values of interest are discussed in §§ 3-9, 3.10, 
and 3-11- One particular average value that is convenient to use is 
defined by taking <f> = £ (eq. ( 3 . 21 d)) and using equations ( 3 * 30 ); 
thus define 

-ills 5 <«^lass (3 ' 31a) 
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or 





I 

o 


f K 


2+Z 





1 1_ 
n<p 2 n 


f U de f t e +l d S 

-n o 


(3.31b) 


(3.31c) 


V 


U) 



1 

n 


n n 

J &<p (sin <p) / d0 / f^ 2 d£ 
o -no 


(3.31d) 


The quantity 1 =4> has already found use in Chapter II. For the 
case £=1, it is often convenient to omit the superscript; thus 


v 


class 


v 


(i) 

class 


<?>c 


olass 


(3.31e) 


3.9 Collision Frequency Per Molecule of Each Class 


For each of the classes defined in § 3-^, the total number in the 

class that undergo collisions in unit time was found in § 3.5- Also 

defined in § 3-5 were the quantities 0 , which represent the fraction 

eras s 

of the total number of the class that undergo collisions in unit time. 
Starting with 0^. , defined in equation (2.12c) as the molecular property 
<f> a one finds from equations ( 3 . 18 ) and ( 3 . 30 ) that, for each class. 


a(l, c> class 


dh 


= 0 


class 


class 


- < 0 t>c 


class 


Thus : 


d(N )+/dH* = 

c 5 5 


0+ 

5 




E 0 


09 


= /0*\ = / ©-> f ? 2 d£ 

\ £'69 _ r 


e<p o 5 


(3.32) 

(3.33a) 

(3.33b) 
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d(N L/cUSL 

C <p <P 



1_ 1_ 

n <j> 211 


n °° 

/ ae / 0+ f 52 ae 
-n o ^ 


(3.33c) 


n n 00 

/ d<p( sin <p) / d6 / 0+ f ? 2 d? (3.33d) 
o -II o € 


(cf. eqs . (2.13c) and ( 2 . 19 )). 


cLN /dN = 
c 


0 = 


- < e ?> ■ 


1 

n 


3 . 10 " Inver s e Coll ision Frequency;" JPer Molecule of Each Class 

For each of the classes defined in § the total number of mole- 

cules undergoing inverse collisions in unit time that end up in the class 
was found in § 3.6. Also defined in § 3.6 were the ’’inverse collision 
frequencies," 0 cq ass for ea °h class, which represent the fraction of the 
total number in the class that undergo inverse collisions in unit time. 
Starting with 0_|. , defined by equation (3.3d), as the molecular property 
4> , one finds from equations ( 3 . 19 ) and (3. 30) that, for each class. 


a(H ;»clas s 


dN 


class 


^class 


- < 9 ?>c 


class 


(3.34) 


Thus : 


d(N' )-*/ dN.^ = 
c K % 

0 - 5 - 

5 

(3.35a) 

d(N’). /dN_ 

c &p 6<p 

G e*p " 

(3.35b) 

d( ^cV /dN< P 5 

e » = < e i>» 

(3.35c) 

cLN'/ dN = 0' 

c 

= < 0 ?> 

(3.35d) 
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3*11 Absolute Mean Free Path for Each Class, Including Directional 
Mean Free Path 


According to the definition (2.1a) (and with use of eq. (3.32)), 
the mean free path measured relative to the observer (or the "absolute 
mean free path"; cf. eq. (2.22)) for each class is 


^class 


_ <«>c 


class 


«>c 


class 


class 


lass 


(3.36) • 


Thus for molecules in the class dN-^ , the absolute mean free path is 


= 5/0+ 

5 5 


(3.37a) 


(For a gas composed of rigid elastic spheres in equilibrium* this free 
path is given by eq. (8) on page 95 in Chapman and Cowling (1961). Note 
that c in their notation is 5 in present notation.) 

For molecules in the class cLN 5 the absolute mean free path is 






0<p 


0<p 


< e ?>< 


5 /0cp 


/ f S 3 d? 

o 

00 

/ e + f ? 2 as 

o € 


(3.3Tb) 


For molecules in the class dN^ * the absolute mean free path is 


A* 

<P 


/E\ /E\ / d0 / f 5 3 d ? 

\ /<p _ V/fr -n o 


< G t^ 


II 00 

/ a 0 / e->- f s 2 as 

-no? 


(3.37c) 
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For molecules in the class dN (cf. eq,. (2.22)), the absolute mean 
free path is 


- 6 <S p 


n n 

/ d<p( sin 4>) f dd J f £ 3 d£ 
o -II o 

n n 

/ d$)( sin <p) / d0 / 0+f £ 2 dg 
o -It o 5 


(3.37d) 


Note that for the classes dN+ and dfh_ one can define a vector 

5 G<p 

mean free path by: 


A* 


->• 

5 


e+ A_v 
? 5 


(3.38a) 


% 


<0k 


a 5 


0 0<p ^ '' e< p 


( 3 . 38 b) 


Either the vector A*^ or its magnitude, A*^ , may be denoted as the 
directional mean free path . Note also that if f does not depend on 

0, then A* = A* . 

0<p <p 

It may also be convenient to define an average absolute free path 
calculated in a different vay, defined by: 


“class 5 < A |>class < 3 ' 3 9> 



' t> 


c 

1 

(3.40a) 

A eq> 

<“?>e» 

(3.40b) 

A 9 = 

<i?> 
X g'<P 

(3.40c) 

A = 

<1> 

(3.40d) 
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3.12 Equations of Change for Molecul ar-Pr operty Distributions 
for Each Class of Molecules 


For any molecular property <f)(r,£,t) one can write an equation of 
change for the property density 4> class = n class < <j>^ lass for each of 

the four classes of molecules defined in § 3.4, equations (3.17). These 


will be of use, particularly for the class dN 


e<> : 


in the directional 


level of description, in § 4.1. For convenience in the following develop- 
ments it is assumed there are no external body forces acting on a molecule 
other than those involved in the isolated molecular encounters, that is, 


f b = 0 . 


3.12.1 Equation of change _for number density in each class 


For molecules in the class dN. 


l ^r a 


on the kinetic level, the equation 


of change for dN-> is given by (3.6b). That equation is divided by 

d to obtain an equation of change for the number density in phase 
r 5 

space } n-> = f(r,g,t) (eq. (3.7) » a form of the Boltzmann equation): 


3f 

9t 


+ V 


(If) 


(04. - a*) f 

5 K 


(3.41a) 


In equation (3.4la) 3 0_>. f is the probable number density of molecules 

lost from (r, ctfp,(t,cL ^ ) per unit time because of collisions, and 
©4. f is the probable number density of molecules gained by (r, aV), 

j, 

( 5 , d9^ ) because of collisions (see §§ 3.5 and 3.6 ). 

For the class dN , equation (3.6b) is integrated over all velocity 

Q<P qy 

magnitudes, £ , from 0 to infinity, then divided by d/£ dQ to 

obtain an equation of change for the number density n (r,e,m,t), the 
probable number of molecules per unit volume of configuration space at 
r and per unit solid angle of molecular-velocity space, in the direction 

0,<p : 


3n 


6<p 


at 


+ V r‘ n e<p n eo/ e e<p 


6(jp 0<P 


e<P' 


( 3 . Ulb ) 
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In equation (3.^1b), n Qyri 0. is the probable number lost, owing to 

u(p U<P 

collisions, per unit time, per unit volume at r, and per unit solid angle 
of velocity space in the direction 0 , qf> ; n e<p ®0<p ^ e V ro ^ >a ^^- e 

number gained, owing to collisions, per unit time, per unit volume at 
r, and per unit solid angle of velocity space in the direction 0, Cp 
(see § § 3.5 and 3.6). 

For the class cLN^ , equation (3.6b) is integrated over all velocity 
magnitudes £ from 0 to infinity and over all 0 from -II to n 5 
then divided by d^£ dfyp to obtain an equation of change for the number 
density n^(r,<p,t), the probable number of molecules per unit volume at 
r and per unit solid angle of velocity space in the direction <p (averaged 
over all 0 ) : 


3n 


<P 


3t 




(3.Ulc) 


In equation (3.Hlc), n^ 0^ is the probable number lost, owing to colli- 
sions, per unit time, per unit volume at r, and per unit solid angle in 
velocity space in the direction of <p (averaged over all 0 ); n<p 0^ is the 
probable number gained, owing to collisions, per unit time, per unit 
volume at r, and per unit solid angle of velocity space in the direction 
(p (averaged over 0) (see §§ 3.5 and 3.6). 

For the class dN, equation (3.6b) is integrated over all that 

is, over all £ from 0 to infinity, over all 0 from -II to n, and 
over all (p from 0 to II; then divided by dY to obtain an equation 
of change for the number density n(r,t), the probable number of molecules 
per unit volume at r: 


3n 

at 


+ V 


[n <|>] = n(0' - 0) = 


(3. Hid) 


In equation (3- Hid), n 0 is the total number density of molecules under- 
going "direct collisions” per unit time, which is the same as the total 
number density of molecules undergoing "inverse collisions" per unit time, 
n 0 1 (see §§ 3.5 and 3.6). 
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3.12.2 Equations of cha nge fo r other molecular- 
property distributions in each class 


For any other molecular property, the equations of change for each 
class of molecules defined in § 3.4 are found as follows: 

For the class dN^- , the density distribution of the property 
(J)(r,^,t) in phase space is (eq. (3.20a)): 


To obtain an equation of change for the distribution <f>f , multiply 
equation (3.4la) by <j> (where <j> may be a tensor of any order): 

|f + <j> v r • (If) = f 4>(e^ - G|) 

and rearrange in the form 

|t- (f+) + V • (ff<j>) - f[ + f-V <J>] = f [ 4>( ©4- - a*)] (3.^2a) 

o u r do r £* ^ 

To obtain equations of change for the property-density distributions 
in each of the other three classes, multiply equation (3.42a) by 

d£ dG sin (}? dtp and integrate first over d£ , then over d0 , 
then over &<p to obtain respectively the results : 

Tt ^ r: S9 "©Xcp ' + 7 r’ - n ecp'4''W' -n 6f»^3t"^W 

= -<* 0 fV (3 - 1,2b) 

It ["<p<*V * v'vCs+V - n *KftV<H+V 

= y <+ e i>,, - ^■ k2 °'l 





I 


If [»<*>] * V r - [=<?♦>] - nt<H> +<«'V> ] 


(3.42d) 


In each of equations (3.42a, b, c, d) the first term on the right 
side is the density, for the class, of the molecular property <J> gained 
by the respective class of molecules owing to collisions, and the second 
term on the right side is the density, for the class, of the molecular 
property tj> lost from the respective class owing to collisions. 

Note that equations (3.42) all have a very similar form and that 
equation (3.42d) is Enskog T s general equation of change (with = 0). 
(See, e.g., Hirschf elder , Curtiss, and Bird, 1964, p. 460 . ) 

1 9 / 

For the particular cases of c(>=m, cf>=mg , and <f>= t? 1 ? (the 
summational invariants of a collision) one obtains from (3.42d) a form 
of the conservation equations of gasdynamics (for which ( 04 - - 0->)^ = 0): 


+ v r ‘ [p<(t^>] = P<0| “ 0 |^> = 0 (3.43a) 

ft + V r - = p<J(0i - 0|)> = o (3.43b) 

ft [ \ P<5 2 >] + v r ‘ t | P<C5 2 I)>] = P<f 5 2 (0| - 0|)> = 0 (3.43c) 
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CHAPTER IV 


PRELIMINARY FORMULATION OF THE 
DIRECTIONAL-MEAN-FREE-PATH METHOD 


4.1 Directional-Mean-Free-Path Approximation for the "Gain- 
Term" Collision Integrals in Flows Far From Boundaries 


In this section a model for the "gain-term" collision integrals is 
proposed that simplifies the treatment on the directional level of 
description; that is, in consideration of the class dN (cf. § 3.4), 
or the class of molecules in (r,d^) with molecular velocity in the 
direction range (0, (P, dfi^) . 

Up to this point, all developments and definitions made have been 
completely rigorous and formal, with no approximations made other than 
those involved in assuming validity of Boltzmann's equation in the form 
of equation (3.7), with Eg subsequently taken to be zero. 

The model to be proposed here involves a physically intuitive approxi- 
mation of the "gain term" of the collision integrals in equations (3.4lb) 
and (3.42b). It makes use of the concept of the directional mean free 
path, introduced in § 3.11. The form of the approximations will be seen 
to be particularly amenable to treatment by certain generalizations of 
Lagrange's expansion, to be developed in Chapter V and to be applied in 
Chapter VI. 

The equation of change for dN , or for n^ , can be written 
as 


d f dl 

¥ r 




+ 



a * N o4<p - d( Ve<p 


0 e<p) 


= d % d %<p V ( w 


(4.1) 
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where d(N c ) is the total probable number of molecules undergoing 

collisions per unit time in (r,d^),( 0 , <P , 9 and hence lost from 

(r,d^ ) and the solid-angle range of velocity ( 0 , (p , dft ) per unit 
r . 04 * 

time because of collisions; and probable number of 

molecules undergoing "inverse collisions" that occur elsewhere in phase 
space (i.e., other than in any n% within (r ,d^) , ( 0 , d^ Q ^) ) 

that are put into (r,d^),(0,^ 5 dfi Q< p) i n time at time t because 

of collisions. 


4,1.1 Definitions and_as su mptions of th e di r ect ional-mean- 
free-path a pproximation for the gain term 

The aim of the direct ional-mean-free-path approximation is to replace 
the difficult gain term on the directional level by an appropriate corres- 
ponding loss term. 

The inverse collisions all occur in elements of phase space, 

back in the direction -e nir . at locations within the solid angle opposite 

0*r 

to ( 0 , dp, dft ) from the point r (see fig. 4.1). 

One may think of the gain term in equation (4.1), therefore, as 
resulting, on the average, from collisions occurring at a directional mean 
free path , A*^ , back from the point r in the solid angle equal to 

d^ 0 (p a " lDOU ' {: ~^Q(p anc * a me an collision time in the past (corresponding 

to this direction), 1/0 . Let the representative midpoint and time of 

the inverse collisions be denoted by r and t, defined implicitly by 


r + A* = r 

0<P 


( 4 .- 2 a) 


t + 1 / 0 e<p = t 


( 4 . Zb ) 


" e e<p ^Ae4>^ e ea> 


(4.3a) 



(4.3b) 


© 


<* 5 V r e ’*’ t) = < e J>e<p 


<^00 = [<£>e<p evaluated at Lt] = f 


(4.3c) 



n 


6<p 


(r,e,q>,t) 


(4.3d) 


in which A* and 1/0 are respectively the directional absolute mean 

0<p % 0<p 

free path and mean collision time that are characteristic of the direction 
0,<j> at the point r and at time t. 

The number of molecules, d(N^ ) q^ > put into (r,dV r ) with velocity 

direction in (e,<?>> dQ a/n ) at time t .per unit time owing to collisions 

0CP 

is accordingly assumed to "be emitted per unit time from collisions in 
(r,d^) with average velocity v^ E v tth direction in (e^jd!^^) 

at time t. The appropriateness of this assumption can be seen from 
figure U.2 and the following considerations: The volume (r,d^) is 

equal in size to the volume (r,d^). Let B denote the area, in the 
plane perpendicular to e at point f, that subtends the solid angle 

dft ^ rom the point r, with point r in the center of area B. Thus 

B = (A* ) 2 dn 

0<P 6 <p 


Let duo be the solid angle, centered also at point r, and subtended by 
the cross-sectional area of the volume (rjd^); denote the cross section 


of 


(r,d9^) 


by dB, so that 


dB = (1*^)2 da) 


Let the ratio of dn^ to doo be j 


thus 


d % = J dw 


B = j dB 


48 





Now divide the area B into j equal parts, dB, so that each subtends 

a solid angle dco from the point r. Consider the small volume dJY 

r 

centered on each of the j small areas, dB. Each of these j elements 

has a cone of solid angle dfi in the direction e . Each of these 

&p dtp 

cones is subtended by a surface A in the plane that is normal to e rt 

0<p 

and passes through the point r. Each of these areas. A, is equal to 
B, and each A may be divided into j equal parts, dA = dB, one of 
which contains the cross-section of (r, a %)■ The number of molecules 
emitted per unit time in ( 0 ) with average velocity v from each 
d f r is assumed to be ^ fraction, 1/j 5 of this number per 

unit time is received by (r,dY r ) at time t = t + l/Q from - each of 

the j volume elements d r on area B. Hence (r, d y ) receives 
j(l/j)d(N ! ) = d(N f ) molecules per unit time at time t with velocity 

C yjCf C 

direction in (e,(P,da ), where d(N T ) Q _ is also the number originally 
assumed to be emitted from (r, d%) by collisions per unit time at time 
t with velocity direction in (6,<p,dQ_ ). 

It is significant that, on the average, the d( ) ^molecules , assumed 
to be emitted from (r , d^ ) in (0,<p,dQ ) per unit time at time t, 

travel the distance A* to point r in the time interval 1/cL =A* /v rt , 

dip 0cp Qcp Qcp 

where 1/0 q^ is the average frequency of direct collisions at t 

in direction range ( 0 , (p , &Q ) . Note that d(N f ) is yet to be 

^ utp c 0<p 

determined. 


4.1.2 Directional-average collis ion models to b e used with the 

direct ional-mean-free-path approximation for t he gain term 
For the direct ional-me an- free-path approximation partially introduced 
above, all the m °l ecu l es have been assumed, on the average , to 

be emitted from collisions taking place within (r, dV) at time t, 
per unit time. This number d(N’)^^ must now be estimated in some way. 
For this estimate, just as (?,d^) was taken to be a small volume at a 
representative average location in which the inverse collisions occur, so 
there is also now assumed to be some representative average di recti on 
that the molecules had before the "inverse collisions" in 
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at time t that resulted in the final average velocity 

in the direction range (0,<p,dQ_ ). The number d(N T )^ is estimated 

d(p c 0cp 

in terms of a representative average direction (0 T ,q>' ) (see fig. 4. 3 ) 

for the d(N')_ molecules before collision in (r ,dy ) at time t , 
c 6<p * r 5 

per unit time, according to one of several "directional-average collision 
models" for the gain term, as described below in §§ 4. 1.2.1, 4. 1.2. 2, and 
4.1.2. 3. These models for the gain term simply represent approximate 
ways of evaluating the representative direction 0',<p f and of relating 
the gain term collision integral, in the directional-mean-free-path approxi- 
mation, to a corresponding loss term , which is easier to deal with. The 
resulting collision integrals must then still be evaluated in terms of 
some process of molecular interaction, with some assumed representation of 
a molecular-interaction potential, to be considered in §§ 4.1.3 and 4.3.3. 
Thus, the number per unit time, d . is to be related to d(N ) t . , 
the number of molecules undergoing direct collisions per unit time in 
(?,d^) at time t and having initial velocities in the solid angle 

range dft , f about e f f , where (cf. eqs . (3.10)) 

0 cp 0 


e'f 


5 e e <? (e ’’ < P’) 


(4.4) 


is the unit vector in the direction 0 T ,<p T . 

For any directional-average collision model to be discussed, one 
will first need to consider the number of molecules in (?,d^) at time 

t with velocity direction in ( 0 f , 4 / ,dft , f ) (see fig. 4.3): 

• 0 


diL,^ f = rT tvvt d Y dO . . 
0 T (p f 0<^ T r r 0 ! cp 


= n e<p (? ’ e ’’f d % d V<p' 


(4.5a) 

(4.5b) 


(where the direction 0',<p' is yet to be defined for the particular 
directional-average collision models). During the time interval dt 
after time t, a fraction. 


0 e'<?' dt 


0 6(p( r ’ e ' ,< * > ' ’O dt 


(4.6) 
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ill i ill i linn i hi i 


nil I III I II I 


of this number dN Q , , undergo collisions with other molecules that have 
average velocity ^ = V(?,t), where V(r,t) = )> . Thus, the number, 

out of the original cLN , , , that undergo collisions in time dt is 

D 4? 


d( VeV dt = e e V ®eV dt 


n e'<p' ®e'<p’ d % dn e'<p' dt 


(4.7a) 

( 4 . Tb ) 


This n umb er is scattered in all directions and with all velocities by the 
collisions. (in the center-of-mass frame of reference of each pair of 
colliding molecules, at least for hard elastic spheres, they are scattered 
isotropically; cf. Jeans, 1952. This has been a useful concept in the 
formulation and interpretation of the BGK collision model; see Liepmann, 
Narasimha, and Chahine, 1962, p.l3l8.) 

For any directional-average collision model to he considered, it is 
convenient to define certain "inverse-collision factors," Z(r , 0 ,<f>,t ) , 

by 


■V <* e |>e«> 


Z n 


eV <? 0 !/ > e , <p' 


(4.8a) 


and in particular 


<+ u>0 A 


ecp 


£/0<p 


U) 7 ~ 


(4.8b) 


where ^ , % =0,1,2 For various collision models (I, II, 

III) one can denote Z = Zj , Z = , or Z = Z^. 

In terms of these defined inverse collision factors, equations (3.4lb) 
and (3-42h) become 



+ e 6V ' V r 

= (0)z n 0y^'^' " n e<p < \ 0 |/ > e<p 


(4.9a) 
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St ^<* u ^ + V'r'“*<* < " 1, V 

- U>z Vr <* U>e !>ev - n e«.<* U) ®|>6* ih - 9b) 

and for any other cj>(r,f\t), 

St [n e<><»Xcp ] + ®e<p ' 7 p [n 8»< 5 *^» ] ‘ n e*[ < clt>et. + 

* z Vf 5®|>e’<p' ' n e ¥ <+ e t)fe«' (l, ' 9c) 

where, for any function F(r,0,<p,t), the notation F indicates: 

F = P(r,0,«p,t) 


where 


r e r - t* (?,6,qp s t) = r - 


<«>e 


0 £ 


e<p e<p ^ 


< e t> fl 


d«p 


t H t - l/<§ > ( 


|/9<P 


(4.10) 


and where the particular collision model considered determines the direc- 
tion 0 T 3 <p f corresponding to each 0 ,<£ and determines the factors Z, 
to he considered helow. 

A particular constraint to he satisfied hy the collision models should 
he noted: Since in local translational equilibrium (with a local Maxwellian 

molecular-velocity distribution function) the sum of the collision terms 
on the right side of the Boltzmann equation (e.g., in the form (3.4la)) 
vanishes, the right side of each of the integrals of that equation, 
equations (4.9) , must also vanish in equilibrium. It is evident then 
that the specification of the collision models should allow, in the 

— y- 

limiting condition of equilibrium, e fl f , to become coincident with e 


^0 


'0<p 
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and each Z to become unity. The task in the remainder of this section 
is to consider various direct ional~average collision models for determining 
and the corresponding factors Z that satisfy these conditions. 


0 f Cp f 


h. 1.2.1 The simplest directional-average collision model, I 


The molecules whose number is given by equations (^.7 ) 5 scattered out 

of dQ , , in all directions during dt by collisions in (r, dV), 

0 <p r 

tend to favor the original direction before collision (even though the 

colliding molecules may be scattered isotropically in the frame of 

reference fixed relative to the center of mass of each pair of colliding 

molecules) because of the effect of persistence of velocities (cf. Chapman 

and Cowling, 1961 ; Jeans, 1952; and Present, 1958). We therefore take, 

for the simplest model to consider, that the representative direction 

before the collisions (0 T ,<p’) is approximately the same as (0,<p). 

Then take as a reasonable intuitive assumption that, on the average, this 

number d(N ) , scattered out of (r, dY), (e,q>, &Q ) per unit time, 

c 6 <P r 

can represent approximately the number put into a volume dV of the same 
size a directional mean free path away at r = ? + A* (?,9,<p,t) per unit 

utp 

time at time t = t+l /0 (r, 0 ,<p,t), with final velocity in the solid angle 

0<p 


( 0 ,<p,dQ e<p ). at r, 


, owing to collisions at point r 
assume for the model that 


at time t . Thus we 


0 ? <p T 0<p 


(model I) 


(4.lia) 


and that 


(0) 


z i = 1 


(4.11b) 


in equation (4.9a) (so that the right side of (4.9a) would become 

n e<p ®e<p ” n e<p 0 e<3> ^ * + + ~ 

Furthermore, assume that certain other molecular properties <J>(r,£,t), 
which are characteristic properties of the molecules at point r at time 
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t, are carried with this number of molecules 
assumed to be put into (r, d$p, (e,<p,dfl ) 

Thus it -will be assumed for this model that, 
equations (4.8) and (4.9) : 


" d( Ve<>- 

in unit time 


for certain 


that has been 
at time t. 

<f>( r,f,t) in 


Z 


I 


= 1 


(4.11c) 


The number and nature of the properties assumed to satisfy (4.11c) 
would be specified as required to make a determined system of equations 
and to provide most realistic results from the model in problems where 
it can be used. 

It is assumed that the main justification for use of this model in 
any given problem would be determined by the realism and consistency of 
the quantitative or qualitative results it yields and by the tractability 
it provides in solving the equations of gas flow where significant trans- 
lational nonequilibrium is present. 


4. 1.2. 2 D irectional-average co llision model II 

The next simplest directional-average collision model to be considered, 
denoted here as II, corresponds in a sense to the BGK model of the 
Boltzmann equation (cf. Liepmann, Karas Imha, and Chahine, 1962 ) (also 
called the Krook-Welander model; cf. Bhatnagar, Gross, and Krook, 1954, 
and Welander, 1954). 

An important aspect of the BGK kinetic model is its interpretation of 
molecules being emitted from collisions in local equilibrium with the 
state of the gas at the point of collision (see Liepmann, Narasimha, and 
Chahine, 1962 ). In the BGK model, the gain-term collision integral of the 
Boltzmann equation is replaced by the local Maxwellian distribution func- 
tion times the local collision frequency; that is, the term 0 ^. f in 
equation (3.4l) is replaced by 0 f^ where f e is the local Maxwellian 
(eq. ( 2 . 23 a)). 
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This same notion is used here to define collision model II, except 
that here the molecules are assumed to he emitted from collisions in 
equilibrium with the local state at point r when considering the gain 

— ^ — y- 

term for the point r. For this model, the initial direction e^ r ^ t at 
point r is immaterial, so we take: 


e 


0’<p’ 



(model II ) 


and assume that the right side of each of equations (4.9) is 


(4.12a) 


^ n e<f^e ^ 0 |^p*U " n 09 < C ct) 0 |) > 0 ^ 


where subscript e indicates a local-equilibrium value, evaluated using 
a local-Maxwellian velocity distribution function (to be described in 
§ 4.3.2 and to be used in Chapter VII). This is equivalent to defining, 
for each in equations (4.9): 


U) 


J II 


n e<p 

n e’<p’ ^ 


^ 9<P h^ U) e^] e 

n e'4>' 


(4.12b) 


It will be seen that the first appr oximat i on in a scheme described 
In Chapter VI gives a result on the directional level for this model, II, 
that corresponds to direct integration of the BGK kinetic equation. 


4. 1.2. 3 Directional-average collision model III 

The following-described model, denoted as III, estimates a direction 
e_ !/r)f that is different from e _ in regions of spatially varying flow. 
The corresponding evaluation of the inverse-collision factor, j ? 

follows directly. 
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For estimating e^, , > the number per unit time coming out of 

collisions, d(N ) , , , given by (U.7), is assumed to have average 

c 0 cp 

velocity e rtl ^. v rt . . within before collision and average velocity 

47 0 T (p B (f 0 <p 

e v _ within dfL,*, after collision (see fig. 4.3). (The relationship 

0q> ocp 0<p 

between dft^ and dS'T will be given below.) Assume that the average 

0<J> 0<p 

velocity of all the other molecules with which this number collides is 
V T before collision and V after collision. An equation that could be 
interpreted as a momentum- conservation equation between "pairs of average 
molecules", one with original velocity e 0f ^ T an< ^ Tinal velocity 

eL,„v„ , the other with respective initial and final velocities and 

e<p e<p 

V, is 


© v + V * 

0 T q> ! e’9 T 


e v + V 

0? 0<p 


(4.13a) 


The quantity 


^ c (r,t,6,<p) 


1' - 1 


e 6<p V 0q> ” e 0 ! <p ! V 0'cp ! 


(4.13b) 


could be called a "relative overall collision velocity change for the 
direction (0,cp) at (r,t)." It is a measure of the degree of spatial 
nonequilibrium, since it vanishes when e^,^, = e g c ^ * One could represent 
V by any appropriate quantity that is known to vanish in equilibrium. As 
a special case that appears particularly appropriate, which we may denote 
as model Ilia , it is taken as a reasonable assumption that 


v =V' - V = V - $ (model Ilia) (U.l4) 

c 


that is , that the average change in the overall mass velocity of the 
"other molecules" colliding with the d(N c )Q f ^ f molecules can be repre- 
sented approximately by V-V, the change in overall mass velocity in one 
directional mean free path in the direction e g^> • The approximate 
validity of this model. Ilia, or of any other model for in (4.13b) y 
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could "be tested by solving the resulting flow equations for a given problem. 

Refer now to figure 4.4, which: (a) represents the vector addition 

of equations (4.13), and (b) illustrates the construction of the solid 
angles for estimating d(N f ) (and consequently Z TTT ) a to be discussed 
in the following. 

Note first that the vector equation (4.13b) 5 along with the trigono- 
metric law of cosines : 


v c = 44 + 4)^- - 2 We-?- e e<p' e e'<r 


(4.15) 


can be solved for the direction e^ , , and the magnitude v f f i n terms 

of e_ _ , e , v A ^ , and V , where 

5 c 5 QCp c 5 


V 


it 


and 


H /V 
c c 


(4.16) 


and in particular. 


itt 


(model Ilia) 


(4.17) 


Thus, taking the dot product of e with equation (4.13b) gives an 

0<P 

expression which, when substituted into equation (4.15), yields 


v 


e'4>- 


[V 2 + v 2 - 2 v V (e *e )] 1 ^ 2 
c e? e? c e<p c 


(4.18a) 


Then substitution of (4.l8a) back into (4.13b) gives 


e -<p- 


e e? v e? - e c v c 


[V 2 + v 2 - 2 v 


V. (S,_4 )] l/2 


( 4 . 18b ) 


c 0<p Q(p c 0<p c 


The direction (0 ? ,<p T ) is therefore defined explicitly in terms of 

V c = e^V^ . (Note that model Ilia could be regarded as a first approximation 


6o 




and then higher order corrections made to ensure that 


for V o , 




-y ~ -y ~ 

e e<jp v e<p _ e e *<p * v e *<p* 


) 


The factor Z = Zjjj for this model is now found by the following 
procedure. If each possible vector of magnitude Vq,-, , and having 

any direction within ( 6 ' ,<f' ,dfig t< p, ) is added to (see fig. 4.4), 

all the resultant vectors would fill a solid angle, denoted above as 

dfL • In general dfi Q/yi is not the same as dfl . The magnitude of 

6<Z) o<P t) 

dfi is found as indicated by the construction in figure 4.4: 

0<p 

The ratio of the area vjL d9 0<p on figure 4.4 to the area v^ t , x 
is l cos “I = l®0<j>‘ e 0'<p' I ♦ Therefore 


d9 


0<jp 


/ V 9 , q>' V 

v V 0(f> ' 


; 0<P* e 0'<p'' d %'<p' 


(4.19) 


The number of the molecules emitted from collisions per unit time in 

(r,d^) at time t in the angle range (e,<p,dfi e ^) is 

times the number emitted in (0,0,dft._ ). Thus, from above assumptions 

(with use of (4.7b) and (4.19)) > 


d(1 4V 


d U 




^ d(N ) , , = 

c 0 '<p ' 

exp 




dft 


0<P 


dfi 


0CJ> 


dfi 0’<y' 


d( "c>e' T ' 


; 'e<p V V«' a 9’a»' da e » 


0'4»' 


6<p 6 <p 


(4.20) 


Now substitution of (4.l8a) and of the scalar product of e Q( p with equation 
(4.l8b) into equation (4.20) gives for collision model III: 
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d(N') 0 -. = n ©' d °/ dQ Q ^ 

c 0^> Q<P 0<p r 0<p 


(4.21a) 


o'q , , 0 Q .^. d*j/ dft Q ._ 
0 t cp l 0 Cp r r r 0< jP 


(X) 

lA v e<p' 


2 v i 1 / 2 

+ 1 - 2(t.t) ° 


(4.21b) 


0<J> c v 


e?J 


v 

l-(e n *e) ~— 


'0<f> c v 


ecp 


Comparison of equations (4.21) and use of the definitions in (4.8) give 
for collision model III: 


<o) r 


in 


17 v c \ 2 

v - 

-1/2 

V 

[(*<*) 

+ 1 - 2( e. • e ) 

e <P c v e<p _ 


1_< ^ 


-1 


(4.22) 


where model Ilia can he used to estimate: 


V 


1 1 

v-v 


and 


e = 


V-V 

1 ?-?i 


(model Ilia) 


(4.23) 


It may he assumed that certain molecular properties cj)(r,^ T 5 t), which 
are characteristic properties of the molecules at point r and time t, 
with velocity having direction (G’^cp*)* are carried with the number 

per unit time, ^^<Pe<p 9 on avera S e » as estimated hy equation (4.21b) 
and therefore that 


U), 


- (0) r 


j iii 


hi 


(4.24) 


for use in equation (4.9b). 

Further reduction of this model, and its use, are illustrated below. 
It is assumed that the main justification for use of this intuitive 
collision model would be determined by the results it yields. 
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4.1.3 Assumption on the collision frequency per molecule, and 


results 


Since A h , defined by equation (3.40b), could reasonably play 
0sP 

the same role as A* (defined by (3.3Tb)) in the development of the 

0 V 

direct ional-mean-free-path approximation above, it is regarded as consistent 
with the rest of the development to consider 0_>. (cf. eqs . (2.12c) and 
(3.33a)) to be independent of the molecular-velocity magnitude , £. 

The results are that 

©| = (b.25) 

and therefore that (cf. eqs. (3.37) to (3.40)) 


A 


6<p 



(4.2 6 ) 


In addition to the above argument, another indication that the assumption 
of 0-* independent of £ is realistic is the fact that for Maxwell 
molecules 0^ does not vary with any part of the molecular velocity. 
Equation (4.25) still allows variation with 9 and <p , as well as with 
r, and t. 

The assumption leading to (4.25) and (4.26) leads to further very 
convenient results, mainly embodied in the form taken by the directional 
equations of change (eqs. (3.4lb) and (3.42b) or eqs. (4.9)) which, with 
the directional-mean-free-path approximation described above in §4.1.1 
and with any of the proposed directional-average collision models (§4.1.2), 
become 


at ^ n e<p < \‘* ) Xp-' + e e<p' v r^ n eqK^X<p-l ~ n e<p "^at’^ > e<p +e e<p' ^r^^X<p 


= z 


n e'<j>’ ^X'qp' ® e ’cp ' n eq>^4^i: 

where for any function F( r , 9 ,<p,t-) : 


"e<$> 0 9<p 


(4.27) 
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F 


F(r,0,<p,t) 



r e e<p 


t 


t - 1/0 


e<p 


^ (U.28) 


where e is defined by equations (3.10) and where the molecular proper- 
ties <j> = <f>(r,£,t) are to be appropriately chosen (see § 4.2 below). 


4.-2 Directional Equations and Some "Boundary Conditions" in Terms of 

Dimensionless Integrals of Directional Property Distributions 

If one chooses to work on the d irectional level , as described above, 
one should still ensure satisfaction of the conservation equations on the 
macroscopic (full moment) level. A proposed method for accomplishing this 
is described here and is used in the following sections (§§ 4.3 and 4.4), 
as well as in Chapter VII. The procedure amounts to defining new dimension- 
less dependent variables, on the directional level, that are certain 
integrals of the previously defined directional property distributions. 

These new variables are defined in a form such that: (a) the directional 

equations of change can be written in terms of derivatives of the new 
variables, and (b) the macroscopic conservation principles can be pre- 
scribed as certain boundary conditions on the new variables on the direc- 
tional level. The procedure also provides a convenient means for 
evaluating all the macroscopic variables of interest. 

Since the macroscopic conservation equations, (3.43) 9 are integrals 
of the directional equations of change (eq. (3.42d), represented approxi- 
mately by (4.27) and (4.28)), with <j> respectively m, m£, and — - m£ 2 , 
it is appropriate to use these choices for (f> in equation (4.27). In 
writing the resulting equations it is convenient to use the following 
definitions (cf. §§ 3.7 and 3.8): 
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Mass density: 


$ e® } = P e<p (r,S,<p,t) = m I f 5 2 d£ 

o 


, » n n 

$ = p(r,t) = / d<P (sirup) / d0 [p 0<(> (r ,0 ,<p,t ) ] 


-n 


Mass flux, or moment-urn density : 


> (1) 




0<P 

e«p p e<p v e<p _ ^e<p 


~-y ( 1 ) -¥ r _ o . 

= e Qm P oyn v A _ = e Q/K m J f 5 dC 

o 


= J(r,t) = P<(t)> = f d<p (sin <p) / d6[j (r,0,(p 5 t)] 


n 


n 


_n 


Moment-uni flux: 


> (2) _ . 

e<p 0cp 


= p 9 T <II> 6? 


-+ -* ( 2 ) -» -> ,°° u 9 
= e e p v = e e„ e f P dr 

0(p 0q> M 0<P 0<p 9 <p 0<j? J o s ^ 


/ \ °° n 

> = P(r,t) = p J d<p (sin <f>) / cL© [ |» (r ,0 ,<p,t ) ] 

o -n V 


Energy density 


.t 


I • = ** h I 0ip ( ^’ e, ' t ’’ t) 


/ 1 p2 \ _ 1 (2) _ 1 r . , 

P 0(p\ 2 ^ /0cp 2 P 0f V 0<P 2 m ' f ■ 


tSee footnote on pp. 32 and 33* 


(4.29a) 

(4.29b) 

(4.29c) 

(4.29d) 

(4.29e) 

(4.29f) 

4 d£ ( 4 . 29g) 
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1 / « \ n n 

2 I : *'■ = E(r,t) = p 2 < 5 2 > = / dtp ( sin <P) / de[E (r,0,<p,t)] (U.29h) 

o -II “ 

Energy Flux : 

2< : *99 5 V ?,9 ’ , ’’ t) s P e<p H 2 «>e<p 

/ \ °° 

e e<p "2 P e<p v e<p _ e e<p 2" m £ f ? 5 d? (4.291) 

I 1 ^ ^ E Q(?,t) = p(|? 2 t)= / d<jp ( sin (p ) / de[Q (?,0,9,t)] (4.29J) 

/->, TT * 


Then the macroscopic conservation equations, (3.43), are 


3p/3t 

+ V 

r 

-> 

• J = 

0 

(4.30a) 

3J/3t 

+ V 

r 

p = 

0 

(4.30b) 

3E/3t 

+ V 

r 

-> 

* Q = 

0 

(4.30c) 


and for <f> - m, m£, and —■ m£ 2 the corresponding forms of equation 
(4.27) are: 


3p Q /n/3t + V 

0<P r 

. t - ( 0 ) 

J 0<P 

3J fliri /3t + 7 ' 

0<P r 

r r-t 

II 

0 ,® 


P e’cp’ 0 0'<p' P 0<p 0 0<p 
J 0’<p’ 0 0 '<p’ “ J 0(p 9 ecp 


(it. 31 a) 
( 4 . 31 b ) 


+See footnote on pp. 32 and 33. 
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3 V at 


+ V 




_ (2) 


Z E 


e'<p' 0 e'<p' 


E 09 0 0<p 


(4.31c) 


where 


-> 

r 



r J 09^00 0 e<p 


(4.31d) 


t 


t - 1/0 


0<p 


(U.31e) 


4.2.1 Dimens io nless variables and equations of change 


For the purpose of defining dimensionless variables that will facili- 
tate solution of the equations, consider the constant reference values: 
arbitrary length L, density p , velocity u , and characteristic 

Ot Oi 

collision frequency 0^ . Define the dimensionless constant: 


e = u /L 0 
a a 


(4.32) 


It is convenient (here and in succeeding chapters) to consider time as a 

fourth dimension, with unit vector e^ that is orthogonal to , e^ , 

and e : that is 

3 


e . 

i 



<5. • (i,j = 1 »2 ,3,4) 


= 1 if i=j 

= 0 if i^J 


(4.33) 


(All the three-dimensional vectors and tensors discussed above can be 
considered as four-dimensional -with zero component in the e^ direction.) 
Then define the dimensionless four-component (space-time) position vectors 
(e.g., see Karamcheti, 1967 ) : 
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->■ 

c 


i=l 


+ _ i - A - V 

e i C i = L r + S 4 ^T 


+ - ^2 + ^3 h. V 

S 1 L ®2 L ®3 L + e 4 L 


t - y - _ 1 * A - V 

= i=l S± Z± = L r + e 4 — 


x, x 


= e 


1 -> 9 ->-3 

1 L + e 2 ~ + e 3 L“ + 


u t 

a 

e 4 L 


and the dimensionless vector partial derivative (holding 0 
constant in a function of t, 0, and <P): 


I e. 

i 




i=i - 


L V + t f- 

r 4 u 9t 

a. 


Also define (for — H < 0 <_ II , 0 <. <p <. II , with v and 

"dummy variahles" for 6 and (p ) : 


F (0) (t,6,(p,e) 


^ l} (t ,0,<P»e) = 


n 0 

f dn (sin n) / dv 

Vn * 

-n 


'<f> 


n e 

jf dn (sin n) / dv 

<P in 



F (2) (t,e,^,s) 


n e 

f dn (sin n) / dv 

* in 


lgg(r»v,n,t) ' 

p u 2 
a a 


G (0) (!,e,»,c) 


n 


e 


L d n (sin n) / dv 

v -n 


E e<p (?>v,n>t) 

p 

a a 


(4.34a) 

(4.34b) 

and <p 

(4.35) 

n, respectively } 

(4.36a) 

(4.36b) 

(4.36c) 

(4.36d) 
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6 


JI 0 Q ( X* ^ 

4t,6,<f>,e) = f dn (sin n) / dv 

9 -n L p »; 


3 

a a 


H (?,e,<J>,e) = 0 Qm (r,0,<p,t)/0 


so that: 




l a 2 F (0 *(e.e.<i>.s) 
sin <p d(f> 3 0 




P u 
a a 


■* (1) 

e 6(p p 9 d> Qcp 

P u 
a a 


i 3 2 F (l) (t,e,q>,t) 


sin <p 3q>3 0 


p e<0 ( ?, e ,(P, t) 


p u2 

a a 


-*■ -»■ ( 2 ) 

e e<p e 9<p p eq? e<p 

p u 2 

a a 


-1 3 2 F (2) (t,e,<p,e) 


sin <p 3 <p 30 


0 -*■ 6<p 

e 9<P e 0(P u2 

a a 


E 64> (r > e ^’ t) 

P u 2 
H a a 


1 _ r (2) 

is 

P u 2 
a a 


2 P e<p v e<p 


-1 3 2 G (0) (t,6,<p,e) 

sin <p 3<p 30 


Qq^( r 9 Q ) 

p u 3 

a a 


-* 1 (3) 

e 9<P 2 p 9<P V 64> 

p u 3 

a a 


-i 3 2 G (1) (c,e,<p,g) 


sin <p 3<p 30 


(where e ajn is given hy (3.10)). 

For any function f ( £ ,0 5 <p, e ) 9 denote 


f = f(z,0,<P, e) 


(4.36e) 


(4.37) 


(4.38a) 


(4.38b) 


\ (4.38c) 


(4.38d) 


(4.38e) 


(4.39) 
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Then define the four- component dimensionless vector 


y (£»8 ,<P,e) 


-^L(r,e,(p,t) 

/ d 

u /0 
a a 


i=l 




^ /e 


8 < p\ u „ 


•(l/H) 



- e 


^ V 0 


e<p. 


/ s 2 ? (l) 

/ 3 2 F^ \ 

+ t 

y 3<p3 e i 

f 909 0 J 

+ e 4 


> (U.uo) 




with which equations (4.31d) and (4.31e) become 


->■ -> 

C + e y 


(4.ia) 


where 


->■ 

U 


_ ->/->• \ .L 

= y (. z , 6 ,<p , e j = - — 

® TT 



(3 2 ? (1)N ] 


1 

\ 9<P90 / 

~^y 

H 

/ a 2p(°) N 

+ % 


_ V 3<p3 0 ) 

f 


(4.42) 


The directional equations of change, (4.31a,b,c), are now: 


? 3 f(0 ’ + , • 

3 2 F (1) 

1 

(0) ZH' 

sin <P 


3 2 F (0) 
— 11 „ 

(4.43a) 

3? 4 3<p3 e c 

3<P3 9 

e 


sin 

\ 3<ps e / 

9090 

3 3 1 ^ 

3 2 F (2) 

1 

(1, ZH' 

sin <P 

/Vf^V 

TT 3 2 F (1) 

(4.43b) 
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3t 9<P30 4 

4 

3<p3 6 

£ 

sin <p f 

y 3<p3 e J 

3<p3 e 


• 3 2 G (1) 

3036 


1 

e 


(4.43c) 


33 G (0) 

3 3<P3 0 


+ V 




sm <p 


where 


3 2 p ( 2 } - - 3 2 G (0) 

3 <ps e e e<P e eq> 3 «> 3 e 


dcpae j 


(4.43d) 


and where ( )' indicates a function evaluated at 0', , according to 

the collision model. We note from the definitions (4.36) that: 


when either 6 = -II or <p = II : 

F (°) = j^ 1 ) = = 0 (4.44) 

and (cf. eqs . (4.29)): when 6 = II an d <p - 0 : 

Cc,n ,0,e ) = p(r,t)/p a (4.45a) 

F (l) = F (l ) (t,n,0,e) = J(r,t)/p a u a (4.4 5 h) 

f(2) _ F^ 2 ^(t,n,0,e) = P(r,t)/p a u 2 (4.45c) 

G (0) = G (0) (t ,n,0,e) = E(?,t)/p o U 2 (U.4 5 d) 

G^ 1 ^ ^ ^(c»n,0,e) = ti(r,t)/p a u3 (4.45e) 


Thus, either one can regard equations (4.45) as boundary conditions on 
the integral functions defined in (4.36), where the macroscopic moments 
(p, J, P, E, Q) must satisfy the conservation equations, (4.30); or the 
macroscopic conservation equations (4.30) can themselves be expressed as 
boundary conditions at 0 = II , <p = 0 on the integral functions 
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(using eqs. (4.45)) in the form: 


9F (0) (t,n,o, £ ) 

3 *4 ' 

3? ( X l(t ^,0 ,£_)_ 

3?, ? 


9G (0) (t,n,o, e ) 

3?,. ? 


F (l) (?,n,0,e) = 0 


• F (2) (t,n,o,e) = 0 

• $ (l) (t,n,o,e) = 0 


(4.46a) 

(4.46b) 

(4.46c) 


Before attempting to solve a given problem using the above formulation, 
one needs yet to: evaluate, or relate in some manner to the other 

variables, the collision frequency 0 (or H) ; obtain a determined 

u<p 

system of equations by some other independent means of relating the 

dependent variables; and prescribe sufficient boundary conditions in the 
->■ 

£ domain. These will be considered in later sections. 


4.2.2 Relatio nsh ips of flow-property distributions and new 
dimens ion less va ri able^ bo_the__more conventional 
macroscopic flow variables 


To facilitate application of boundary conditions for solving a given 
problem by the above formulation, and also to enhance physical understand- 
ing, it is desired to relate the macroscopic flow properties in equations 
(4.29) or (4.45) to the more conventional macroscopic flow variables. 

With the mean mass velocity of the gas (?) defined by equation 
(2.15), the molecular velocity relative to the mean motion (?) defined 
by equation ( 2 . 16 ), and the velocity magnitudes defined by 


5 = ill » (f-f ) 1/2 

V = |V| = (^) 1/2 

c = [ c I = (c-c ) 1 ^ 2 


^ (4.47) 
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one has 


3 = p<J> = pV (4.48a) 

P = p <tf^> = p^f + p<J?]> (4.48b) 

E = p<| |-t> = | p<C 2 > + | P V 2 (4.48c) 

Q = P^T f -tf> = \ P<C 2 C> + V • p<CC> 

+ | <C 2 )> pV + | p V 2 V (4.48d) 


The energy per unit mass relative to the mean motion, or the specific 
internal energy, is 

e = | <C 2 > (4.49) 

Since in equilibrium the internal energy is •— kT per degree of freedom 
of a molecule, where k is Boltzmanns constant, the temperature of a 
monatomic gas of neutral spherically-symmetric molecules in arbitrary 
translational nonequilibrium is defined by 

me = 3 (| kT) 

or 

RT e ! e = | <f 2 > (4.50) 

where R = k/m is the gas constant. 

For a perfect (sufficiently rarefied) gas (implicitly assumed in the 
use of the Boltzmann equation), the "collisional transfer” of momentum or 
energy due to finite size of molecules is neglected in comparison to the 
flux of momentum or energy due simply to the "flow of molecules” (e.g., 
see Hirschf elder , Curtiss, and Bird, 1964), so that the stress tensor 
and heat-flux vector are simply: 
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a = - p<cc> (it. 51 ) 

X = \ p<c 2 ^)> (u. 52 ) 

The hydrostatic pressure is defined by 

p* = - i | : <7 =| p<C 2 > (it. 53) 

so that from equation (U. 50 ) we have 

p* = pRT (4.5*0 


The viscous-stress tensor, T 5 is defined by 


C = - p I + T 


(4.55) 


where p is the thermodynamic pressure. Far from translational equilibrium, 
the "thermodynamic pressure" has no meaning. Therefore, since sufficiently 
close to equilibrium the Chapman-Enskog theory for a perfect monatomic 
gas gives a zero bulk viscosity, k, defined by 



it is assumed that the "pressure" in (4.55) is 
p = p* = pRT = |p<C 2 > 

and therefore that the viscous stress tensor is 
T = - P<c£> + I | p<C 2 > 


( 4 . 56 ) 


(h. 57 ) 


(fc.58) 


T5 



With the definitions in equations ( U • ^4-9 ) to (4.58), equations (4.48) 
become 


5 = P v 

P = p ^v - a 

= p^ + Ip - r 


(it. 59a) 


(b. 59b) 


E = 


p(e + | V 2 ) 


3 ^(e + \ V 2 ) + q - ^ • O 


= p!(e+|? 2 ) +1-1 -T + p! 


(i+. 59c) 


(i+.59d) 


(With equations (4.59), equations (4.30) give more familiar forms of the 
macroscopic conservation equations.) Equations (4.59) will be useful 
in determining boundary conditions on the functions in (4.45). 

The common macroscopic variables of interest can then be calculated 
in terms of p, J, P, E, and Q as: 


* = 

j/p 


(4.60a) 

p = 

(2/3)pe = | (E - 

\ pV 2 ) 

( 4 . 6ob ) 

RT = 

(2/3)e = p/p = 

f (E/p - f V 2 ) 

(4.60c) 

<7 = 

-P + p^ 


(4.60d) 

r = 

-P + p^ + pi 

-P + pl^ + j (E - 

| v 2 )l 

(4.60e) 
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(4.60f) 


q = Q - pV(e + \ V 2 ) + V • O 

“)■ -*>■ Q-> 

= Q - EV - V • P + pV 2 V 


Other macroscopic variables of interest are the specific enthalpy, h, 
and the Mach number, M, defined by 


h 


e + p/p 






(4.60g) 


and 


M 2 


v 2 _ 9 
(5/3)RT 5 



(U. 60 h) 


(where the ratio of specific heats for the perfect monatomic gas is 5/3) . 


k . 3 Indet erminacy of the Eq ua tion s o f Change, and a Proposed Method of 

Closure using Int eg rals^ of_ a_ I^cal-Di rectional Gaussian Distribution 
Function 


It is well known that the conservation equations on the macroscopic 
level are an indeterminate system, containing more unknowns than the number 
of equations. In the directional -mean- free-path method introduced above, 
the equations on the directional level of description contain one more 
unknown quantity than the number of equations if the collision frequency 
is appropriately specified. 

The purposes of this section are to: (a) propose a method of closure, 

on the directional level, that emphasizes the molecular-directional aspects 
of translational non equilibrium and that is compatible with the known 
local-Maxwellian distributions of molecular velocities in limiting regions 
of local translational equilibrium; (b) determine the expressions for the 
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directional property distributions in local translational equilibrium for 
applying local-equilibrium boundary conditions; and (c) consider 
representations of the collision frequencies. 

It is useful to consider what may be called a "local-directional 
Gaussian distribution function." First the property distributions on the 
directional level are to be expressed by integrating such an assumed 
distribution function, for use in an "integral method," to be described. 

As a special case , these property distributions are then found in terms of 
a local Maxwellian distribution of molecular velocity for use in applying 
boundary conditions in regions of local translational equilibrium. 


4.3.1 Assumed local-diifectionaljlaussian distribution function , 
and resulting directional-property distributions 

For making the directional equations of change a determined system, 
one additional equation that relates the variables appearing in those 
equations( (4.31) or (4.43)) may be found by what may be called an "integral 
method." Integral methods entail assuming the form of a certain function, 
and then using integrals of that function to advantage, without necessarily 
implying that the original precise variation, in the assumed form of the 
function integrated, is physically significant. In this same sense, a 
"local-directional Gaussian distribution" is now chosen in the form 


f 


-c 2 E 2 + 2 b c £ 

a e 


where 

a = a(r,e,<p,t) 

b = b(r ,6 ,(p,t ) 
c = c(r ,0 ,<jp,t) 


a e 


-6U-U 


e ‘ /2 V 5 


6 


1/2 


(4.61a) 
( 4 . 6lb ) 


} (4.62) 
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are unknown functions of 0 and q> as well as of r and t. Note that 
tj here is an unknown vector function having the dimensions of velocity. 

The form (4.6lb) is especially convenient. This form for f is specified 
only for the purpose of integrating it to obtain an appropriate relation 
among the directional property distributions. Its use (as in other integral- 
approximation techniques) does not imply that the molecular-velocity 
distribution at any (r,0,*p,t) is thought to be Gaussian. However, the 
assumed directional-Gaussian distribution is especially convenient in that 
it becomes the correct l ocal Maxwellian distribution in any limiting 
region of local translational equilibrium, where a, 5, and U become 
independent of 0 and <p (5 becomes V(r,t), 3 becomes 3^ = m/2kT(r,t), 

and a becomes n(r,t)[m/2II kT( r ,t ) ] 3 ^ ) . A usual procedure in studying 
nonequilibrium phenomena (e.g., see Vincenti and Kruger, 1965 , p. 225) is 
to use to advantage certain relations that are known to be true in 
equilibrium and assume that they hold approximately in arbitrary non- 
equilibrium; one may then ascertain that the use of such assumptions 
does not adversely affect the nonequilibrium results obtained; in other 
words that the relationship between the mathematics and the physics of 
the problem is not overly sensitive to the arbitrary assumptions made. 

If one defines the dimensionless quantities ; 


A = A(r,0,<jp,t) 


2 p (cu ) 3 
a a 


a m u 


(?) 


(It. 63a) 


c = c(r ,6 


cu 


(it. 63b) 


and ( cf . eg. (2.2k)) 


B(b) = 


j n 1 ^ 2 e b2 (l + erf b) 


(4.63c) 


and puts equation ( 4 . 6lb ) into the definitions of directional property 
distributions in equations (4.29), one obtains 
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(1+2 b 2 ) B(b) + b 


(4.64a) 


A(r,0,9,t) = 

A — - — = -=• [b(3+2b 2 ) B(b) + 1 + b 2 ] (4.64b) 

U a c 

A — ^ — = [(3+12b 2 +4b 4 ) B(b) +b(5+2b 2 )] (4.64c) 

u 2 2^ 

a 

v eto ) (^’ e ’ < P’ t ) i 

A — — = [b ( 15+20b 2 +4b 4 ) B(b) + (4+9b 2 +2b 4 )] (4.64d) 

u 3 2c 

a 


One could regard equations (4.64a,b,c) as three equations for the three 
unknowns*: A, b, and c. "In principle", the three equations could be 

solved for A, b, and c as functions of v^iVu and 

0<p a 

These results could then be substituted into equation (4.64d) to give 

v[?J/u 3 in terms of v^Vu and vi^ /u 2 . Then one could find, from 
04> a 0<p a 0<jp a 

equation (4.6 3a), 


r ( *W 

04? a 


tt?) 


( ( 


/p )c 3 


a = 


e<p' 2a 


(4.65) 


In this way, two things would be accomplished in principle: 


(a) A relationship among v 
determined; and 


(i) 


v 


( 2 ) 

Q(p 


and v 


(3) 

0<p 


would be 


(b) The quantities a, b, and c would be determined in terms 


of p e<|> 5 


V 


(i) 

ecp 


and v 


( 2 ) 


A perhaps less obvious, but more tractable, procedure for determining 
equivalent relationships is the following: By considering the four 

equations, (4.64), all at once, one can first determine b and c 
explicitly in terms of v^^ , v^^ , and v^^ ; then find A and B 


(3) 


also explicitly in terms of v^^ * V 0^ 9 and - v 0<f 


Then the expressions 
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,(') v < 2 > 

eq> * etp * 


for both b and B in terms of v^_' , v^_' , and , along with 


/ 3 > 

ecp * 


the definition (4.63c)., provide an implicit equation relating v 


( 1 ) 

0Q> 


, and v^ 3 ^ (cf. the statement (a) above). Again, since A and 

e» — e<j> 

c are determined 9 equations (4.65) could again be used to determine a 

,(1) ... ..(2) ....... ..(3) 

'up * 


in terms of p , v' ^ 7 f and v^ 7 (since v^' is implicitly related 


to and v^* ^ ) * so that a, b, and c are all determined (cf. 

0<P 0<p 5 

the statement (b) above). Following are the details of this more tractable 
procedure just described: 

Combining equations (4.64a 5 b s c), one finds 


Q<? 


( 2 ) 


0<p 


0qp 


V (2) 

e 2 I&£- 


v (l) 

- bTT - 22 - 


u^ 


u 


3 

2 


(4.66 a) 


and combining equations (4.64b 5 c 3 and d) , one finds 


T < 3 > 

c2 - 23 - 


v (2) 

-be - 23 - 


= 2 


U" 


u z 


e<p 

u 


(4.66b) 


From equations (4.64a and b) one finds 


v 


(1) 


(c -23 b ) A - 2b B = 1 

u 

a 


(4.66c) 


and it is convenient to write equation (4.64a) as 


A - (1 + 2b 2 )B = b 


(4.66d) 


If now b is eliminated from equations (4.66a and b) to obtain 


c = u 


m- 


L*> 


(3) 

0<f 


2 v < 2 > 

2 e<p 


- 11/2 


- ( v »Tj 


(4.67a) 
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then (U.66a) gives 


c2(T W >/M a ) - 3/2 


2t (') t < 2 > -1 t ( 3 ) 

6<P 6Q) 2 eq> 


Sr 2 /v (i) ) 2 -i v (2 )l L^V 3 ). / v ( 2 ) f] 

|L V Qt? ) 2 e«p J [_e<?> v etp \ 0 q> / 


Equations (4.66c and d) are easily solved to obtain 


1 - (l+2b 2 )(l + b 2 - Vc v\J / u ) 

0<p a 


b(l + b* - be v^/nj 

1 - (l+2b 2 )(l + b 2 - be - v^/u ) 

0<|> a 


where b and. c are given by (4.67a and b). Now if we define: 


* <3) = ♦ (3) (?.e,4>,t) 3 ^/(v^) 3 
* <2> - * (2) (?, s , i (>, t ) $ v ^/ tv ';') 2 


we can, as described above, use equation (4.63c) with (4.67d) to obtain 
the following convenient form of the desired relationship among the 


directional 


distributions that is used to close the system of 


equations on the directional level: 


i n 1//2 e fe2 (l + erf b) 


b(l + b 2 - bs) 

1 - ( l+2b 2 ) ( 1 + b 2 - bs! 


where (from (4.67a) and (4.67b)): 
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2 <|/ 2) - | / 3) 


2-|*< 2 >][*' 3 >-(* (2 >) 2 ]J‘ 


v (1) 

c -§«- 

u. 

a 


3 ( 2 ) 1 1/2 

2 - f * l ' 

,.( 3 ) /,.( 2)\ 2 


-(t (2 >)' 


Implicitly, equations (4.69) define a direct functional relationship 

( 3 ) ( 2 ) 

between and 9 which may be written in the following form 

that will be found convenient for later use: 


^ (2) = vU) 


where 


4> = ( / 


(2)\V,( 3 ) 


This relationship can be computed from (4.69) and tabulated, once and for 
all, for use in the method described above. One method for computing and 
tabulating the functional relationship (4.70) would be to specify various 
values of b ranging from -00 to and, for each of these values of 

b, calculate s from (4.69a) as 


1 + b 2 + (3b+2b 3 )B(b) 
b + ( l+2b 2 )B(b ) 


(4.71) 


then find both ijj 


(2) ( 3 ) 

and ^ by solving (4.69b) and (4.69c) to obtain 


(fc.72) 


v + - 

; 2 S 
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and evaluate iJj from (4.70b). 

Other relationships that are equivalent, but more convenient for later 
use, are found as follows: From the definition of the function B(b) in 

(4.63c), 

<3B(b)/db = 1 + 2b B(b) (4.73) 


and 


N(b) 

= ^tbE(b)] 

OO ry 

=2 / e _v e 
0 

2bv 2 a 
dv 



= b + (l+2b 2 )B(b) 


(4.74a) 

N' (b) 

s fb»(b) - 

OO o 

r -v 2 2b v 

4 j e e 

o 

v 3 dv 



= 2[l + b 2 + 

(3b+2b 3 )B(b)] 


(4.74b) 

W" (b) 

= ~~2 ®^(b) = 

db> 

4[ |- b + b 3 + 

( 1 + 6b 2 +2b 4 )B(b)] 

(4.74c) 

I (k) (b) 

= ~r" [W(b) ] 
db k 

OO o 

= 2 1+k J e~ v2 

o 

e 2bv v 2+k dv, k=0 ,1,2,. 

(4.74d) 

from which 

can be found the differential 

equations : 



N"(b) = 2b N' (b) + 6 N(b) 


(4.75a) 


and 


N (k) (b) = 2b N (k_l) (b) + 2 ( k+1 ) N (k 2) (b) 


(k=2,3 ,4, . . . , ) (4.75b) 


Thus, it is seen that equation (4.71) is simply 


s 


N'(b)/2 N(b) 


(4.76) 



for use in (4.72), or in fact , from ( 4 . 64 ) , (4.68), and (4.70) * that 


(l/ 2 ^ = ¥(b) ¥"(b)/[¥' (b)] 2 (4.77a) 

/ 3) = [¥(b)] 2 ¥' ”(b)/[¥' (b)] 3 (4.77b) 

tfi = [¥"(b)] 2 /K' (b) ¥' • '(b) (4.77c) 


For large negative b, it is extremely difficult to calculate B(b) 

( 2 ) ( 3 ) 

to sufficient accuracy to obtain ip 9 ip , and ip. In that case, 
the functions can be calculated simultaneously by numerical integration of 
the ordinary differential equations (4.73) and (4.75), starting at large 
negative b with values of the functions found from the asymptotic expan- 
sions: as b -°o : 


B(b) % l 
n=0 


(~l) n+1 ( 2n) ! 


n! (2b) 


2n+l 


Uli ^ y 

db k n=l 


( 2 ) 


( ~l) n+k ( 2n+k ) ! 

2n / -i \ . . 2n+k+l 
(n-1)! b 


(k=0 ,1,2, ... ) 


(4.78a) 


(4.78b) 


The accuracy of this evaluation of the functions B(b) and d k N(b)/db k 
is then easily checked by comparing the values obtained at b=0 with the 
exact values : 


B(0) = /rf/2 , N(0) = /if/2 , W'(0) = 2 


(4.79) 


and with the values of the higher derivatives of N at 
directly from (4.75b) 

The variation of i|/ 2 ^ and ij/ 3 ^ with b can be 
by asymptotic expansions. For example, one finds: 


b=0 found 
examined locally 
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as b -*■ 


/ 2 ) -V 

•| [ 1 - l/2b 2 + 15 Ab 4 + 0(l/b 6 )] 

(4.80a) 

„/ 3 > * 

[1 - 3/2b 2 + 0(l/b 4 )] 

( 4 . 80b ) 

as b 0 : 



A - 

(3n/8)[l - (3n 1/2 - 28/3IT 1/2 ) b + 0(b 2 ) ] 

(4.80c) 

* (3) - 

(n/ 2 )[i - ( 2 in 1/2 /8 - 8/n 1/2 ) b + o(b 2 )] 

( 4 . 80d) 

as b +°° : 




^ 1 + l/2b 2 - 3/2b 4 + 15 /4b 6 - + ••• (4.80e) 

i|/ 3 ^ ^ 1 + 3/2b 2 - 9 /2b 4 + l4/b 6 - + ••• (b.QOf) 

Use of equations (4,68) through (4.8Q) is illustrated later. Figure 4.5 

( 2 ) ( 3 ) 

shows , i/j 3 and \p versus b. 

4.3*2 Special case of local Maxwellian velocity-distribution function 

In order to apply boundary conditions in regions of local translational 
equilibrium, it is convenient to express the directional property distribu- 
tions in terms of the macroscopic variables. If it is assumed that in 
regions of local translational equilibrium the assumed directional Gaussian 
distribution, (4.6l) , becomes the local Maxwellian distribution function, 
(2.23a) , then the functions a, B, and U in equations (4.6l) and (4.62) 
become independent of 0 and <p . Comparison of equation (4.6la) with 
(2.23a) then shows that, as a region of local translational equilibrium 
is approached, we have the asymptotic relations: 
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r “13/2 

a(r ,0,(j>,t) ^ n(r,t) ^~T"~ 1 

L 211 kT (r,t)J 


( 4 . 8la) 


U(r,0,<p,t) -v V(?,t) 


( 4 . 8lb ) 


3(?,6,<p,t) ^ m/2kT(r,t) 


(4.8lc) 


Then, with 3=3 (r,t) defined by equation (2.23b) and with unit vector 
e e 

e = e(r,t) = V/V (eq. (4.13)), equations (4.62) and (4.8l) give, as local 
equilibrium is approached: 


-3 V 2 

a(r,0,<p,t) n(3 /n) 3/2 e 6 


(4.82a) 


b(r,e,<p,t) 3 l J 2 (e 


+ x i /2 

c(r,e,<?,t) 'u 3 e 


\) V 


(4.82b) 

(4.82c) 


and the definitions (4.63a) and (4.63b) give, as equilibrium is approached: 


3 V 2 


/ p V 

A(r,0,<J>,t) 'u 2n 3/2 e 6 (p fl /p) 


(4.82d) 


. 1/2 

c(r,0,<P,t) 'v 3 e u^ 


(4.82e) 


With equations (4.82), equations (4.63c) and (4.64a,b) now give as 
equilibrium is approached : 

— 8 y2 

p Q/ft (?,e,cp,t)i) a V (p/p a )(l/2n 3/2 ) e e (3 3/2 • e V 


0<J> 


+ [1 + 23 e (e 0<p -e) 2 V 2 ] B(B* /2 t Q(f - e V)} 


( 4 . 83a) 
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and 


-6 V 2 


D V 

M 0(P 0<f 


^ ) (r, 0 ,<f,t)/p a u a * (p/p a )(l/2n 3/2 B* /2 u a ) e e a 


+ & e ( ®e<j>‘® )2v2 + B e /2 ®e<p‘* v[3+2 B e ( ®0<f'® )2v2] B(s i /2 *oq>‘* v)} (U * 83b) 

where the function B(h) is defined "by (4.63c). Equations (4.66a) and 
(4.66b) then give also: 


p 0<i> V 0tf )(? ’ e ’ < * > ’ t) 3 (p 0<p /p a ) , *6*'® V 

p u 2 2 6 u 2 u. 

a a e a a 


( 1 ) 


P 0«p V 0<p 

P a U a 


(4.83c) 


p 0» v 0» (? ’ 9a ^ t) ^ _2_ 


p u^ 
a a 


3 u z 
e a 


p 0<j> v 0q> \ e 0tp’ eV 


P u 
a a 


u 


a 


(2) 

P o<p v 0<p 

p 

a a 


(4.83d) 


4.3.3 _E s t imated or ass ume d forms ^of collis ion frequencies 

In § 4.1.3 it was assumed that 0| is independent of velocity 
magnitude so that 0+ = 0 0<p . To calculate both 0^ for use in 

the directional equations of change (4.31) (or in (4.37) and. (4.40) - 
(U.43)) 9 and 0 for the characteristic value 0^ in (4.32) and ( 4 . 37 ) * 
one now needs expressions of assumed or estimated forms of these collision 
frequencies per molecule of the respective classes (cf. § 3-9)- 

For Maxwell molecules (with the intermole cular force assumed to vary 
inversely as the fifth power of distance between molecules) one finds 

0. = 0 = 0 A = 0 = k ' p(r,t ) (4.84) 

If one takes k ! to be a constant, then (cf. eq. (4.37)) 

H = - fillet). ( 4 . 85 a) 

0 p 
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This relation could be assumed as qualitatively appropriate for use in 
the directional' -me an- free -path method. It may be more realistic, however, 
to use a relation equivalent to what has been used in the BGK model (see 
Liepmann, Narasimha, and Chahine, 1962, pp. 1319 and 1321, where their A 
is proportional to our k ! ; see also discussion by Vincenti and Kruger, 
1965, p. 384): 



(4.85b) 


where n is a coefficient of viscosity to be evaluated from some appro- 
priate temperature-dependent viscosity law. For Maxwell molecules, the 

bracketed factor in ( 4 . 85 b ) is unity, to give (4.85a). For rigid spheres, 

1 /2. 

the factor is (T/T ) 7 , which is obtainable also by use of equation 

(2.20b) with (2.23d). Liepmann, Karas imha, and Chahine ( 1962 ) used the 
Sutherland viscosity law for n . In ( 1966 ) they used n « ip 0 * 8 l 6 . The 
latter is appropriate for argon gas up to 4000°K (see Camac, 1965* p. 248) 
Any other appropriate calculation of 0 from equation (3.33b), 
where 0 -> is given by ( 2 . 12 c) , could also be used. 


4.4 A Determined_System of Equations in a Simplified Form of the 

Direct ion al-Mean-Fre e-Path Method for One-Dimensional Flow 

Flow in one dimension is simplified by the fact that the velocity- 
distribution function depends on only one configuration-space variable, 
x^ , and time t, as well as on £ and <p but is independent of 0 . 

The resulting directional-level and full-macroscopic-level equations can 
then be written simply in scalar form. (There is only one non-zero velocity 
component and one non-zero component of heat flux on the macroscopic 
level; and only one component of the viscous-stress tensor is of particular 
interest.) In that case it is convenient to reformulate part of the 
method described above in somewhat simpler terms, as described in the 
following. 
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h.k,l ^Mmensi onle s s variabl es fo r one- di mensional flow 
Let 


a) = cos q> ( U . 86 ) 

Then, in one-dimensional flow, where the distribution function f is 
independent of 0 (and where the velocity vector is V = e 1 u) ^ we have 


->■ -> -* 
p • p =: p • p 

0<P C 0<p 1 


cos <p = 0 ) 


(^. 87) 


If the one-dimensional flow is allowed to be time-dependent , the components 
of X that are of concern are 


^ = x^/L and 


? . = U t/L 

4 a 


(see § 4.2.1). We now define the new dimensionless dependent variables 
(with n now as a dummy integration variable for id): 


p(c»w,e) 


J( C till) E ) 


P( t, tw j e ) 


E( ^ ,u»e ) 




2n / I — ) dn 
-1 \ p a / 


2n / 


w / p,p v <p 


(i) 


_i \ p~ u 


a a 

( 2 ) 


2n / 


<P <P 


-1 \ p a u a 


n dn 


n 2 dn 


0) 

2n / 

-1 


2 V 

p <p> \ 

) dn 

P XL / 

a a / 


0) 

2n / 

-1 


1 / p «p> 


( 3 ) 


p u 3 *" 
a a 


n dn 


( 4 . 88a) 

(4.88b) 

(4.88c) 

(4.88d) 

(4.88e) 
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which are related to the functionsdefined in equations (4.36), evaluated 


at 6 = II , by: 

F^(c»n 5 <p,£) = (4.89a) 

1 ^ (t,n,<p,e) = J(t»ai,e) (4.89b) 

F^ 2 4t,n,<p,e) = e l e 1 P(t:*o 5 e) 

+ (® 2 e 2 + [ E (taW,e)- P(?»a>,e)] (4.89c) 

G^(t,n,<p,e) = E(t,ai,e) (4.89d) 

(^(t.It^e) = e l Q (t.u.e) (4.89e) 

We then have: 

o_/p = (1/2JI) 3p73w (4.90a) 

<P a 

p v^Vp u = (l/2n oj) 3J/3u) ( 4 . 90 b) 

<f K ct a 

P vj, 2 V p u 2 = (l/2n oo 2 ) 3P/3to (4.90c) 

H (p <f K a a 

= (l/2n)(2 9E/9oo) (4.90d) 

P vi 3 Vp u 3 = (l/2n co)(2 30/3(1)) (4.90e) 

a a 

We note from equations (4.88) that 

—y 

for co = -1, all X, : 

-=J = P = E = Q = 0 (4.91) 

and from equations (4.89), (4.45), and (4.59), 
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for a) = +1» all j : 


p* = p"U,l,e) 
J = J(t,l,e) 
P = P(c,l,e) 

E = E(cA,e) 
Q = Q(?,L,e) 


p(x,t)/p 

Utr 

(l/p a u a ) pu(x,t) 

(l/p^U 2 ) [pU 2 + p(x,t) - TjjCx.t)] 
(l/p a u 2 )[pu 2 - ff u (x,t)] 

(l/p a u 2 )[pe(x,t) + j pu 2 ] 

(l/p a u3)[ P u(e + ^-u 2 ) + q.(x,t) - u(x n -p)] 


From equations (4.92) one finds (ef. eqs . (U.6o)): 

= ptt. i, E ) 


u(x,t) 

u 

a 

= ( — 
\ p > 

) J(c.l,e) 

p(x,t) 

p a u a 

2 

3 1 

if-J(k) - 

RT ( x , t ) 

u 2 

a 

2 

3 

a L 

f n (x,t) 

p a u a 

- Cr)Ct) - FiU ’ e) 

x a 7 ' a ' 

: n (x,t) 

2 

(r)(5r) 2 + §^-‘> 

x a 7 x a 7 

p u 2 
a a 

3 1 

q(x,t) 
p u 3 

= Q(?,l,e) -(^~) [E(t,l,e) 

V a ' 


or ' a 1 


a a 


(4.92a) 

(4.92b) 

(4.92c) 

( 4 . 92d) 
(4.92e) 

(4.93a) 

(4.93b) 

(4.93c) 

(4.93d) 

(4.93e) 

(4.93f) 

(4.93g) 
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(U.93h) 


h(x»t) _ 5. RT 

u 2 2 U.2 

a a 



(4.931) 


which can he used to evaluate the conventional macroscopic variables after 
p% J, P, E, and Q are determined. 

For use in the directional equations of change for one-dimensional 
flow, in which for any function F(e,a>,e), ve denote (cf. eq. ( U . 39 ) ) 

i = F(z 9 u,e) (k.9k) 


where 


Z = C+ty (4.95) 

we have now (cf. eqs . (4.40)): 

= ®l (h If/If) + + + *4 | > ^’96) 

(The components u 2 and y 3 are not significant, since there are no 
variations in the x 2 and x 3 directions.) 

b.k,2 Equations of change for one-d-imensional flow 

The directional equations of change (U.31a,b,c), for one-dimensional 
flow in terms of the functions in equations (4.90) (with use of eqs. (4.38)), 
become (cf. eqs. (4.43)) 


9U 



(4.97a) 


£ 


£ 


£ 



P 

3u) 


J 

3a> 


E 
3 co 


+ 


+ 


+ 


a 2 j \ 


3 2 P \ 
3?1 9 “/ 


9 2 Q \ 


(0> a' 


^ 1 ^ ZH ' 


^zh' 




3J 
3 a) 



where 


9E = 1 9P 

3 co 2u)^ 3 oj 


(4.97b) 


(4.97c) 


(4.97d) 


where ( )' indicates a function evaluated at co' = cos <p' according 

to the directional-average collision models discussed in § 4.1.2; where 

( ) is defined by equations (4.94) - (4.96); and where 

— — / . , -¥• . 

Z = Z (C,u),e) = Z(r,0,<j>,t) for one-dimensional flow. 

The macroscopic conservation equations (4.30) become, in terms of 
the functions in equations (4.88) evaluated at oj = 1 (cf. eqs . (4.46) 
and (4.89)): 


3 

N 

F(t,l>£) 

+ a, 

J(t,l»e) 

= 0 

(4.98a) 

9 

^4 

J(t,l,e) 


P(C,l,e) 

= 0 

(4.98b) 

d 

^4 

E(c.l,e) 

+ 1 

Q( C »1 » E ) 

= 0 

(4.98c) 


4.4.3 Directional-Gaussian closure equation 

The equation used to close the system, derived in § 4.3 by relating 
various integrals of a local direct ional-Gaussian distribution function, 
is given by equations (4.69) and expressed implicitly by equations (4.70). 
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The quantities and ijj defined by equations (4.68) and 

(4.70b) are written in terms of the one-dimensional-flow variables as 


( 2 ) 

ip = v. 


( 2 )., (lh2 _ (9P/9 o))(9p/9u>) 

* 3 “ (9J/9to) 2 


(4.99a) 


,(3) _ (3), r ( 1 ) \ 3 2w 2 ( 9Q/3to) ( 9p'/9aj) 2 

ip = T <p /'■ v a> ' — 

V V ( 9 j/9co) 3 


(4.99b) 


* = (+< 2 >)W 3) - — 


2u) 2 (9Q/3w) ( 9J/9co) 


(4.99c) 


Note in particular that ip does not contain 9p/9w , so that if the 
quantities on the right side of (4.99c) are known, one can: calculate 

ifj from (4.99c); obtain i|r from (4.70a); then obtain 3p/9o) from 
(4.99a). This procedure will be used later. 


4.4.4 Local-equilibrium property distributions 

In regions of local translational equilibrium, where the molecular- 
velocity distribution function, f, is Maxwellian, the directional property 
distributions (eqs. (4.83)) for one-dimensional flow are simplified. 

for a monatomic gas is /$/ 6 M, and 


1 / 2 

The quantity g ' V 


e *e = cos <D = a). Then (cf. eqs. (4.62) and (4.82)), 


b = b = Moj = 


0) 


(4.100a) 


where 


M = ^76 M (4.100b) 

and equations (4.64), or (4.83), (with the definitions (4.63c) and (4.47)) 
become 
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(4.101a) 



(4.101b) 


(4.101c) 


(4.101d) 


(4 .101e) 


where subscript e denotes a local equilibrium value. Equations (4.101) 
can be integrated with respect to oT (with use of eqs . (4.75) where 
needed) to obtain, in regions of translational equilibrium in one-dimen- 
sional flow: 


P e (? jU.e) 


r . -m 2 i 

[W / n l/2 M_ 


[u> B(m) + M B ( -M ) ] 


(4.102a) 


J e (c sm»e) 


n r -M 2 

1 [n' / 2 M 3 . 


[u N(co) + M N(-M) - ^B(rn) - M B(-M)] (4.102b) 


P ( C * e ) 0 


7p v -M 2 

(**) n l/ 2 M 5 „ 


[u 3 N(m) + M 3 N ( -M ) ] 


E (5,a),e) = 


Q e ( ? , e ) 


-l T KT I ( _ 1\T ' < _M "I 


e ~, [N'(to) - N'(-M)] 

/2 M 3 _ 


#);? 


(4.102c) 


(4.102d) 


N' ’ ' (-M)-IO N*(w)+10 N'(-M)] 

(4.102e) 
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Equations (4.102) have been made to satisfy the following conditions 
(which are consistent with conditions (4.91) and (4.92) and the definitions 

in (4.93)): 


at co = -1, all g : 


p = J = P = E = Q = 0 

e e e e e 


(4.103) 


and 

at qj = +1, all Z : 


p (i^l.e) = 


J e (?»l,e) 


p(t»l,e) 


= p/p. 


j(C,l,e) = pu/p^ 


P e (c*l,e) 


= P + P u 


P a 2 
a a 


p a u a \ M / 

(■•f) 

^4 (l + \ 

P a\ \ M 2 ) 


2E U,l,e) = 2E(c,l,e) = 


2Q e (? 5 l 5 e) 


2 pu(h + - n 2 ] 

p u 3 
a a 


pu‘ 

P a 2 
a a 


( 4 . 104a) 
( 4 . 104h ) 

(4.104c) 

(4.104d) 

(4.104e) 


Therefore the quantities 
(eqs. (4.100) to (4.102)) 


p/p and M in the local equilibrium expressions 
a 

are, from (4.104a), (4.104b) s and (4.104d): 


p/p 


a 


M = 


ST/Z M = 


3/2 


r2E(t,l,e)" 

P ( £ »1 5 £ ) 

_ P ( C 5 £ ) - 

_j(?,l,e)_ 


ZT 


-1 



(4.105a) 

(4.105b) 


which may he e-valuated in a nonequilihrium solution of the equations. 
These expressions are useful later. 
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4 • 4 ^R eduction of direct ional-avera ge collision models 
for one-d imensional flow 

The quantities in equations (4.97) that are denoted by ( ) f are 
functions of co evaluated at co* , where, for one-dimensional flow, with 
f/V = e = 2 it 


CO = COS (p = 

e • e 

0<P 

<0* E COS q) ! = 

-> 

e • e . 
e f <P 


(4.106a) 


and where, in general. 


(o - co 


1 ( £, }(i)> e) 


(4.106b) 


( £)— 

However, for collision model I, oj ! = 1 and Z = 1, so equations 

(4.97) are simply written without the Z factors and without the "primes” 
on the gain terms. 

For collision model II , with equations (4.12), equations (4.97) 
become: 


e ( 3 2 


8 2 


CD 

W 

II 


H 1 


T 

8(0 


3 to / 

\ 3(0 / 

- H \ 

£ /_al 

JL + 

3 2 

JL\ 

- TT 1 

f^e) 

TT / 

V 3 «4 

8(o 

3 h 

3 to/ 

— IX ' 

e 

\8co / 

“M 


/ 9 2 e + _a 
£ \3?„ 3to 5? 


2 Q \ = g _ H /3e\ 

3to ) e \3co / \3o>/ 


(4.107a) 


(4.10Tb) 


(4.107c) 


3E 

3to 



(4.107d) 


where 3p /3to , 3J / 3 03 , and 3E /3w are given by equations (4.100) and 
e e 0 

(4.101), with P/P a an d M given by (4.105). 
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For collision model Ilia ( from eqs . (4 .l 8 b) , (4.22), (4.23), (4.24) , 
and definitions in (4.90) and (4.92)): 


0) = 0) 


'U,u),e) = 


“ (u-u) 


[(u-u) 2 + - 2co v^tu-u) ] 1//2 


“ [l-A(t,co,e)] [l-A(?,co,e)] 2 '+ 1_ “ 2 ' 


W 


(0 


and 


(4.108) 


U) 


Znia ( ^a),e) 





-1 







1-co" 


-1 


(4.109) 


where 


A( ^ 5 C 1 ) ) 


U-U 



3p(z,o),e) 

3o) 


J(z,l,e) JU,l,e)l 

[dJ( z a a),e) 

L 303 J 


_p"(z>l>e) (T(c,1,£)J 


(4.110) 


The use of the equations given in this section (§ 4.4) is illustrated 
below in Chapter VII in steady-flow shock-wave structure. The equations 
can first be appreciably simplified by a suitable approximation scheme, 
developed in Chapter V, with its application in the present method indicated 
in Chapter VI. 
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CHAPTER V 


VECTOR GENERALIZATIONS OF LAGRANGE’S EXPANSION, 
AND A PERTURBATION-EXPANSION SCHEME 


$.1 Introductory Remarks 

Lagrange* s expansion is a generalization (first derived By Lagrange 
in 1770) of Taylor *s expansion in which the independent variable is defined 
by an implicit equation. Recently Sack ( 1965 a) has reviewed a number of 
approaches to, and applications of, Lagrange’s expansion and has discussed 
generalizations of certain forms of the expansion to several variables 
(see also Sack, 1965b, 1966; Sturrock, i960; and Good, i960). 

In addition to the forms that have been given, an especially simple, 
direct, and useful generalization to N variables of the standard form of 
Lagrange’s expansion can also be found, and is derived here. It can be 
written completely in vector form, which is useful in application to' vector 
formulations of physical problems, especially where transformations to 
different coordinate systems are to be considered. 

Two simple derivations of Lagrange’s expansion with one independent 
variable are first outlined (§ 5-2), and the development of a perturbation- 
expansion scheme shown (§ 5*3). It is then shown (§ 5-^-l) how each step 
in the first derivation of the standard form of Lagrange’s expansion in 
§ 5-2.1 can be simply generalized to two independent variables; and the 
analogous steps in the generalization to arbitrary N dimensions then 
become clear, as shown in § 5-^-2. An alternate derivation which is a 
generalization, to N independent variables, of the alternate derivation 
outlined in § 5-2.2, is then given in § 5-^*3. It is advantageous to 
consider both these derivations because: (a) the first derivation results 

in the more concise final form for one independent variable (although the 
expanded forms are termwise identical) , (b) the alternate derivation 

results in the more concise final form in arbitrary N dimensions (although 
the expanded forms are termwise identical, and (c) .the two derivations are 
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independent and are believed to be both new. The development of the 
perturbation-expansion scheme for arbitrary - N independent variables 
follows in § 5*5 "by analogy with the one- dimensional development in § 5-3. 


5.2 Lagrange's Expansion in One Variable 


If an independent variable z is defined implicitly by the equation 


Z = zU,e) = c + Cy(z) 


(5.1a) 


where ey(z) is sufficiently small, the standard form of Lagrange's 
expansion (cf .Whittaker and Watson, 1927, or Sack, 1965 a) for any function 
f(z) is 


f ( z ) = fU) + l 

n=l 


n _n-l 
c d 

n! 


£—[ y “( C ) 

a c n -iL dc J 


( 5 . lh ) 


5.2.1 First derivation 

A simple derivation of equation (5*rb) is given here (for convenience 
of comparing the steps in the later generalizations). 

Differentiation of (5- la) gives 

dz = d£ + ey 1 (z) dz + y(z) de 

from which 

dz = 1 , 7 - v dc + r — 7 - y de (5.2a) 

l-ep'(z) s l-ey'(z) 

Also 


dz 


z dC 


+ z de 
e 


(5.2b) 
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(where z = 3z/8£ and z = 9z/3e) so that, by comparison of (5-2a) 

E 

with (5-2h) ; or, equivalently, by combining (5.2a) with (5.2b) to obtain 

{[l - ey*(z)]z - l}d? + {[l - ey'(z)]z - y(z)}de = 0 (5.2c) 

and by equating the coefficients of d£ and of de in (5.2c) to zero 
(because the differentials d£ and de are arbitrary , so that equation 
(5-2c) must apply for any dc and any de independently, the coefficients 
of d£ and of de must vanish); one obtains 



For any function of z: 


z = u(z) 

e 1 - ey ' ( z) 

F(z) = F[z(c,e)], one may write 


(5.3) 


3F[ z( C «e) ] = dF 

9e e dz 


(5.4a) 


and 


9F(z) _ dF 

9C dz 

so that, with use of (5*3) » equations (5.4) give 

9F(z) _ Z e 9F(z) _ , . 9F(z) 

9e z^ 9? 9c 

From (5-5). 


9 2 F 

9e 2 


9y(z) 9F(z) 



+ y(z) 


9F(z) 

3? 


9F(z) ~ 
3 ? . 



(5.4b) 


(5-5) 


(5.6) 
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Further, by induction, one finds 


3F(z) 3y 2 (z) j 

3C 3e ) 

, 3F(z) * 3y 2 (z) ) 

+ y(z) ” ~\ 

(5.7) 


3 n F(z) 


9e 


n 



(5-8) 


A direct derivation of equation (5* It) that uses (5-5) through (5*8) 
is now the following: From Taylor* s expansion about e = 0, one writes: 


f(z) = f[z(c, e )l = f[z(c, 0 )l + l 

n=1 ‘ L ae 4=0 


(5-9) 


Since at e=0, z=c , use of (5-8) in (5-9) gives directly 


f (z) 


f(?) 


+ 


00 


l 

n=l 



(5.10) 


which is equivalent to the result sought, (5 * lb). 

The steps in this derivation from equations (5.1a) to (5-5) and from 
(5.8) to (5.10) are equivalent to the corresponding steps in Laplace f s 
derivation in 1780 (e.g., as outlined by Sack, 1965a). However, the steps 
from (5*5) to (5.8) do not employ one additional identity that was used by 
Laplace, but are carried out quite naturally without the additional identity, 
and in a manner that will be seen to provide a simplification in an analogous 
procedure in the generalization to N dimensions, considered in later 
sections . 
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5.2 .2 An alternate derivation for one independent variable 


An independent, alternate derivation of Lagrange’s expansion in a 
form that is termwise identical to (5-lt>) is the following: With (5 -la)* 

one may write Taylor’s expansion as 


f(s) = ,<c) ♦ I 

_ n! , r n 

n=l dC 


- *<>♦ I £ [»(.)]“ ^ 


Therefore also 


00 m 

m(z) = y(c) + I Jr d y(Al 


m=l 


ml 


dC 


m 


(5.11a) 


(5.111) 


Although not as concisely expressed as equation (5.11), equations (5>ll) 
give the same result as equation (5.11) to any desired order. For example, 
to order £ 2 , (5 -Ha) gives 


f(z) 


f(C) + ey(z) 


df(Q 

d? 


+ 



[y(z)] z 


d 2 f ( ? ) 


0(e 3 ) 


(5.12a) 


where, from (5. Hi) 

y ( z ) = y (^) + ey ( z ) d j ^^ + [ y ( z )] 2 - + 0 ( e 3 ) 

= y (?) + e [y ( C ) + sy ( z ) d ^^ — + 0 ( e 2 )] d ^L 

+ £■ {[yU)P + 0(e)} + 0( e 3) 

= y( C ) + ey( C ) + a 2 jy(0 

+ \ [ y (?)] 2 d d ^^ - | + 0 ( e 3 ) ( 5 . 121 ) 
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thus ( 5 - 12 a) becomes 


f(z) = f (?) + evU) + e 2 jvi(?) 

+ |[y (?)] 2 } + 0 (e 3 ) ( 5 . 12 c) 

which is equivalent to ( 5 - lb) to order e 2 . 

Obviously, the form of equations (5.1l) is not as convenient as the 
concise form ( 5 -lb), but this derivation has been given because its 
generalization to N independent variables (§ 5-4-3 below) provides a 
more concise final form than does the generalization of the first deriva- 
tion. 


5,3 Lagrange f s___Expans ion in One Variable with Additional Parameters , 
and Development of Perturb at ion- Exp ansi on Scheme 

In a development exactly analogous to that in § 5-2, it is easily 
shown that equations ( 5 -la) and ( 5 . 1 b) may be generalized to include an 
arbitrary number of parameters; thus, if 


z =5 z 


ep ( z,a ,6 9 y 5 • • • ) 


(5- 13a) 


then 


(5- 13b) 


f(z,a,e*Y» • • • ) = - - - ) + £ TT “ n T [ C , B ,y 5 - - • ) ~ " " " ^" l 

n=l dr L ^ J 


Equation (5.13a) has the expansion 


«) n _ n -1 

z = ? + £ nT — i>T ^ C ,01,6 ,Y . * - - ) J 

n=l ’ 3 5 


(5.13c) 


Consider now the special case of equations (5-13) in which only one 
parameter, a, is included, and suppose f(z,a) and y(z,a) can be 
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expanded in power series in a : 


f(z,a) = l a k 1 f (z) 
k=l 

y(z,a) = l a k_1 y,(z) 
k=l 

Then equations (5*13b) and (5.13c) "become 


f(z,a) = 


+ 


where 


z = C + 


l 


n=l 


l 


n=l 


l 


n=l 


n- 

a 


n 

e 

3 n_1 

n! 

9C n_1 

n 

e 

3 n_1 

n! 




Lk=l 


a k 1 y k (c) 


-in 


These equations may be further specialized to the case where a =e, 
which the terms of equations (5*15) to order e 3 are: 


f(z,e) = + e[f 2 (?)+y 1 (?)f|(?)] + E 2 {f 3 (^)+y 1 (^)f^(?) 

+ y 2 (‘c)f{(c) + (l/ 2 )(< 3 -/ <i ^)[yf (?)f{(?)]} + e 3 {f 4 (?) 

+ y 1 (?)f^(?)+y 2 (?)f^(?) + y 3 (?)f|(s) 

+ (i/2)(d/a c )[ y 2( c )f'( c ) + 2 Vi U)y 2 U)fjU)] 

+ (l/6)(d2/d C 2)[y3( c ) f .( c )]} + 0 ( e 4) 


z = ? + ey^?) + e 2 ty 2 (c)+y 1 (dy;(0] + e 3 {y 3 (?) 

+ y 1 (c)y 2 (c) + y 2 (c)yj(?) + ( 1 / 2 )( d / d ?)[y^(c)y|(c)]} + o( E 4 ) 


(5-l4a) 

(5.14k) 

(5.15a) 

(5.15b) 

for 


( 5 . 16 a) 

( 5 . 16 b) 
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(Equations ( 5 . 16 ) have been found to be useful in applying Lighthill's 
uniformization technique; see Martin, 1967 a.) 


5.j 4 Generalizations of Lagrange T s Expans ion_to __N Variables 


The standard form of Lagrange T s expansion (eq.(5*lb) with (5.1a)) 
can be generalized to N dimensions in vector form. The independent 
variable z (cf. eq. ( 5 . la)) is defined implicitly by 

z = z (t,e) = C+ey(z) (5.17) 


where z, £ , and y are vectors each having N independent components, 

and thus are considered to be vectors in the N -dimensional space having 

orthogonal unit base vectors e^ (k=l, 2 , . . . ,N) , in the directions of the 

respective z coordinates (e.g., refer to Karamcheti, 1967 ) . The vectors 
K 

may therefore be written as 


-> 

z 


N 

I \ z k 

k=l 


C 


N 





k 


N 

y - ,7 e k y k (z l ’ Z 2 ’ ' ' ' ’ Z N ) 


( 5 . 18 a) 


( 5 . 18 b) 


( 5 . 18 c) 


where the unit vectors satisfy the orthogonality condition (which also 
defines the "dot product"): 


e. 

1 



6^. = 1 for i = j 


= 0 for i ^ j 


(5.19) 


in which i and j may have values from 1 to N. 
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Then, for any function f(z) = f (Zj ,z 2 , . . . ,z ) , (which itself may 
he a scalar, vector, or higher order tensor), the result (to he derived 
helow) for the extension of the standard form of Lagrange's expansion to 
N dimensions, in vector form, is 

f(z) = fit) + e{y(t) * fit)) 

+ ( e 2 / 2 ! ) (y • V (y- V f) + (y-V y) • V f} 

Q> Q> 

+ U 3 / 3 !)(y -V [y-V (y-V f) + (y-V y) • V f] 

Q Q *=> 

+ 2 (y • V y ) • V (y-V f) + [y-V (y-V y) 

+ (y-V y) • V y] • V f} + 0 (e 4 ) ( 5 . 20 ) 

hs hs 


where z is defined implicitly by equation (5-15); where the argument of 
each function of the right side of (5*20 ) is % \ and where 


V 

C 




( 5 . 21 ) 


Formulas for the higher order terms will be given below. 

It is to be noted that if y is a constant vector, then equations 
( 5 . 17 ) and ( 5 . 20 ) reduce to the vector form of Taylor’s expansion in N 
dimensions (cf. Korn and Korn, 1961). If, on the other hand, N=l, so 

that f and y are functions of only one variable, z, equation (5-20) 

reduces directly to the standard form of Lagrange’s expansion in one 
variable, equation (5.1b). 

To derive equation ( 5 * 20 ) (to any order in e) , one can proceed 

in a manner directly analogous to that in either § 5.2.1 or § 5*2.2, as 

demonstrated in the following. 
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$.4.1 The two-dimensional case, N = 2; first derivation 


The generalization for N=2 is carried out most simply, which will he 
outlined here. The extension to arbitrary N then follows directly. 

All major steps in this section are analogous to those in § 5-2.1. 

From equation ($.17) 9 for N= 2, 


z i = z 


1 (? 1 ,C 2 ,e) = Cj + ey 1 (z 1 ,z 2 ) 


z 2 Z 2 ~ ^ Z 1 9 Z 2 ^ 


Then 


( 5 . 22 ) 


dz. 


9y l 

d h + e ST 


dz = d t; + e 


3y 2 

7 dz, 

SZj 1 


ay x 

e 3z 2 



iZ 2 * *1 

de 

3y ? 



+ £ 

dz 2 

dz 2 + y ? 

de 


from which 


(’ 


jr) dz i - e jr dz 2 = d ^i + de 


-e 


9y 




dz, = dr + y de 


} (5.23) 


Also: 



3z 


3 z i 


+ 3e 

dz 

1 


d ^i 

+ i ^7 

d£ 


9Z 2 


3z 2 


3z 

dz 

2 

= ^ 

d ^i 

9? 2 

d£ 

+ >r 


de 


de 


i (5.24) 


which may he substituted into equations (5*23) to obtain: 
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I 



(5.25a) 


(5.25b) 


The factors multiplying d^, d^ 2 , and de in equations (5.25a) and 

( 5 . 25 b) must vanish (as in eq. ( 5 . 2 c); see discussion of eqs. ( 5 - 2 ) 
above), so the six equations for the six unknowns, Sz^/ 3 ^ , , 

SZj/Se > 8 z 2 /3 ? 1 , 3z 2 /3c 2 , 3z 2 /3e , yield, with 


3 ^ 


9z l 

£ 3z 2 



3y 2 

9y 2 



3 z ^ 

1 E 3z 2 



9y l \ 

(l 

p 2 1 

9y 2 

3zJ 

V 9z 2 / 

3z 2 

3z 1 


( 5 . 26 ) 


the following: 
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II I I 


I llllll I III 1 1 I I II I II I | II ii min mi ii ■ 


3 z 1 

1 

(*-•5) 

ii 

I 1 — * 

1 ^ 

| ro 

D 

9z i 

1 

[(‘-•St 

9 s 

D 

3 z 2 

1 

(i 

9? 2 " 

D 

V 


l!i _ 1 E 
3? 2 D 


+ e *r 2 U2 


9z 2 = 1 ^2 

D £ 3z J 




y 2 + £ 9^7 *lj 


For any function of z only , 
F(z) = F[z(t,e)] = F[: 


F[z l (? 1 ,5 2 ,e) , z 2 (5 1 »C 2 »e)] 


3F( z) 

3z-^ 

9F 


9Z 2 

3F 

9 e 

3e~ 

9z i 


9e 

3z 2 

3F(z) 

3z 1 

9F 

+ 

9z 2 

3F 

3?1 

= ^7 

9z i 


^7 

3 z 

2 

9F{ z) 

9Zi 

9F 

+ 

9z 2 

3F 

9 ^2 

9 ^2 

9Z 1 


^7 

3z 2 


from which, with use of equations (5.27) , it is readily found that 


Mil = 3F(Z l- Z 2 ) 

3 e 3 e 


" VVh’ — h 


3F(z. ,z„) 3F( z, ,z„) 


Mfli 


y(z) • V ? F(z) 
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Equation (5-29) may be differentiated as follows: 


iMH . + [»£L] 


*3* ■ ™ \ [W ♦ ' '< [« 

[^1 ■ \ 


and by induction one finds 


i -1 

9 F(z) _ v 

Z 1 - 1 ) s •);<?) 

• V7 

’ a 1 J 1 f(z) 

i ^ 

9 e j =0 

\ i 8 , J 

• V 

c 

L Se 1 -')- 1 J 


where the Binomial coefficient is defined by 


«■ 


a! 


( a-b )! b ! 

Equation (5-30) is then used as follows: 

Taylor f s expansion about e=0 for f[z(t,c)] is 


f(z) = f[z(t,e)] = f[z(t 5 0)] + l 


r 8 n f[z(t, £ )]- 


n=l 


n! 


3e 


n 


“ J £ = 0 


Since at e=0 , z=C > this is also 


f(h = f(D + 1 4 r 

n=l n! L 3e n J 


£=0 


where the terms in the summation are found (using 5-30)) from: 


(5.30) 


(5.31) 


(5.32) 


(5.33a) 
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9 1 f[z(c,e)] 
i 


i r - *’ / 

<3 ~ jr .z(- 
9 e 

and also 


i-1 

I 

j=0 


'i-1' 

J / 


9 J y[z(c,e)] 
9 s' 3 


9 1 ~ J ~ 1 f[z(t,e)] ~| 

9e i ”* 3_1 J 


( 5 - 33 b) 


9 y[zU,e) 
9s 1 


i-1 / i-l> 

l f 

j=0 V J ) 


9 j y[z(t,e)] 


Se 1 - 




(5.33c) 


and where 


z = z(c,e) 


t + sp(z) 


(5.34) 


Each of the terms on the right side of equations (5.33b) and (5.33c) can 
be found using equations (5- 33b) and (5- 33c) for a smaller order i. Thus, 
to any order, the terms of equation (5.33a) must be found successively from 
equations (5.33b) and (5.33c). The procedure is illustrated as follows: 

(a) First take i=l to find, from equations (5.33b) and (5.33c): 


9f(z) 

9e 


3y(z) 


9e 


= y(z) • V ? f(z) 


= y(z) • V y(z) 


(5.35a) 

(5.35b) 


Since z=i" at e=0, we have from (5.35a) for use in (5.33a) 




tit) • fit) 


(5.35c) 


(b) Next take i=2 in equations (5.33b) and (5«33c) to obtain 


9 2 f 


9e" 


= y (z) 


raf(z) 1 

. f 3 y(z) 1 

L 3 e J 

L 3 e J 


and, with use of equations (5* 35a) and (5- 35b), 
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(5.35d) 


9 2 f(z) 

3e 2 


y(z) 


V r [y(I)*7 f(I)] + [y ( z) • V y ( z) ] • V t(z) 


Similarly 3 


9 2 y(z) 

3e 2 


= y(z) • V [y(z)-vl(z)] + [^)*V r y(z)] * V r y(z) 


(5.35e) 


Since at e=0, we have from (5*35d) for use in equation (5*33a): 


L J e=0 


y(t) ‘ V Jy(t) ' V f(?)] 


+ [y(?) • V ? £(£)] • til) 


(5.35f) 


(c) Higher order terms follow similarly. Thus, for N=2, the 
"standard form" of Lagrange’s expansion in vector form has heen derived, 
the first few terms of which are given hy equation (5*20). 


5.4.2 Arbitrary-N dimensions;, first derivation 


For arbitrary N the major steps in the derivation are directly analo- 
gous to those above (§ 5.4.1) for N=2. Differentiation of equation (5 . IT): 


= Z(t,e) = t + ey(S) 


(5.17) 


(with arbitrary N in equations (5-18)) gives 


-y -y ~y -y 

dz = d£ + e dz • V y + y de 

z 


(5.36a) 


from which 


dz • [ | 


~y /~y \ -1 
e V z y (z) ] 


= d£ + y de 


(5.36b) 
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where 


N - a 

V = y e — 
z . L - k 3z. 
k=l 1 


(5.37a) 


and where I is the idemtensor (unit tensor) 


N N 


I = I I e. e. 6. . 
i=l j=l 1 J 1J 


i=l 


e. e. 

i l 


(5.37b) 


(cf. eg. (5.19)). Also, since t = z(t»e), one may write 


-> -* 3 2 

dz = d ? ' V z + — de 

^ o £ 


(5.38) 


Substitution of eguation (5-38) into (5-36) gives 


at • {(V z) • [ I - e V y(z)] - | } 

+ de{ • [ I - e V y(z)] - y} = 0 (5.39a) 

d e z 

Equation (5- 39a) is a vector equation with N components. In each of the 
N equations for the components, the coefficients of de and of each d£^ 
must vanish (as in eqs . (5-25) an<i (5*2c)). One therefore obtains 
N(N+l) equations to evaluate the N(N+l) unknowns, dz^/d^^ and 
dz^/de (i,j = 1 , 2, ... , N), as follows: 

From ( 5 • 39a) : 

d£ • (V z) ‘ [ I - e V jj ] = d£ (5.39b) 

Q z 

and 

• [ I - e V y] = y (5.39c) 

□ e z 


u6 



From ‘(5.3913) 



9z 

1 3 Cl 


+ 


9z 

3 5, 


+ 


+ 


-> —» 

3z | 

9 Sd 


[ I 



By expanding this equation into its components and setting the coefficients 
of each d^ k equal to zero, one obtains finally a set of N 2 equations 
for the N 2 quantities 3z./3£., (i,j = 1, 2, ... , N) 9 that may be 
written in the form: 



(i,j = 1, 2, ... , N ) 


( 5 . ^Oa) 


Similarly, from (5* 39c), one obtains N equations for the N quantities 
9z^./3e (j = 1, 2, ... , N) that may be written in the form 



(i = 1, 2, ... , N) 


(5A0b) 


Equations (5-^0a) and (5.U0b) now represent N(N+l) equations, consisting 

of N+l sets of N equations each. The N(N+l) unknowns 9z /3^. and 

^ J 

dz^/de can be found by Cramer's Rule (see, e.g., Hildebrand, 1952) from 
these N+l sets of N equations each. Define 


a ik E 6 ik - e 9 V 9z k 


(5.^1) 


and define the minor of a. . as M. . and the cofactor of a. . by 

ij ij ij 


A. . = (-1) 1+J M. . 

ij lj 


(5.^2) 


Denote the determinant of the matrix a^ by |aj . Then 
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a 


and 


IT 

£ a ik A ik 


k=l 


N 


> (5.43) 


-I , L ®kj \j 


k=l 


and also 


N 


7 a . A._ = 0 for r ^ i 

L rk lk 


and 


k=l 


IT 


> (5.W 


l “ksVl = ° f ” 8 * 3 


k=l 


Then Cramer T s Rule gives the solutions of equations (5.40b) and (5.40a) 
for each j as 


ac j 


i 

fal 


N 


y a. 6.. = 

i=l lr 1J 


[aT A jr ^ r,J 1,2,m 


( 5.45a) 


and 


dz 

2 

d£ 


1_ 

I a I 


N 


T A. y . 
1=1 




( 5.45b) 


Now, since for any function 


F(z) 


— F[z(^,e)] — F [ Z ^ ( L, 1 5 • • • 5 £-^5 G ) 5 ^ 9 • • • 5 £ ) 5 • • • 5 Z -j^l 


one may write 
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(J = 1 , 2 , ... , N) 



N 

l 


r=l 


8z 


r 



8 F 

8z 


r 




and 


3F 

8e 


N 

l 


r=l 



8 F 

3z 


> (5.46) 




substitution of equations (5*45) into (5*46) gives 


N 

8 F 1 y 8 F 

8 C . [ a | n jr 3 z 

j 1 r=l ° r 


(j = 1, 2, ... , N) 


(5.47a) 


8F 

8e 



N N 

I l 

r=l i=l 


lr 


8 F 

8z 


r 


(5.47b) 


Upon multiplying equation (5.47a) by y. 

J 

N, one finds 


and summing over j from 1 to 


3F(z) 

dc 


N 

I M,(z) 

j=l J 


3F( z ) 


(5.U8a) 


This result may also be written in the form 

9 ^ 2) = y(z) • V F(z) (5. Mb) 

de £ 

The remainder of the development for arbitrary N is identical to that 
from equations (5-29) to (5*35). Thus, an extension of the "standard 
form" of Lagrange’s expansion has been derived for arbitrary N dimensions, 
with the first few terms given by (5-20). 

It should be noted that in equation (5.17) and the result (5.20) the 
appearance of e is extraneous (since it appears only in the product ey) , 
and so e could be taken to be unity. The nature of the convergence of 
(5*20) then depends on each of the components y^ being sufficiently 
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small (rather than each evu). Each y ^ could then he replaced in (5.17) 
and ( 5 - 20 ) hy e^y^ > and since the y^ T s are independent, the e ^ ’s 
could he considered as representing different orders of magnitude (such 
as = e , = e 2 log e , = e 2 , = e , etc.) (Actually the 

same results had originally heen found hy carrying out the entire derivation 
starting with arbitrary different ! s, rather than the simpler deriva- 

tion in terms of one e given above. Since the results are equivalent , 
the simpler derivation that uses only one e has heen given here.) 


5.^.3 An alternate derivation for* N independent variables 


A generalization, to N independent variables, of the derivation 
given above in § 5-2.2 is described here. The final result is somewhat 
more concise than that given in §§ 5*^.1 and 5.^.2 above. The expanded 
terms are identical. 

It is convenient to first introduce some additional vector and 

tensor notation not previously used. Consider the vectors z, y , and 

the Cartesian base vectors e^ defined as in equations (5.17) to (5.19 ) 9 

and the vector operators V and V defined as in (5.21) and (5* 37a). 

£ z 

For arbitrary N -dimensional vectors 


1 


N 


l 

k=l 



B 


N 


I 

k=l 





B,_ 


(5.49) 


define the nth order tensors 


i (n) s m...t , s (n) = m...s (5.50) 

n times n times 

(We might call and "polyadics, tf since the special cases for 

n = 2 , 3 , and 4 are known respectively as dyadics, triadics, and 
tetradics; cf. Morse and Feshbach, 1953.) Further, define the scalar 
products (which follow from ( 5 * 19 ))- 
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(5.51a) 


A*B 


N 


I A. B. 
l i 


i=l 


AA:BB = 


+ + +* N N 

A* (A*BB) = y y A. A. B. B. 

x ,1 J i 


i=l j=l 


(5.51b) 


AAA:BBB S 


+ + + ^ N N N 

A* [A* ( A'BBB) ] = l l l A A ABB B 

i=l j=l k=l 1 J k k J 1 


(5.51c) 


and, in general, define the following notation for the nth scalar product : 

A (n) (?) B (n) = y y ... y a a a b b b (5.51a) 

* -1 • _ -i • -t -L -L ., • • • 1 _ 

1=1 i=l 1 =1 12 n n n-1 1 

1 2 n 


Then, in particular. 


N N N 


[y(t)] (n) = £(?)£(?)...£(?)= y y ... y e. ...t y . (t) 

n times 1 1 =1 X 2 =1 1 n =1 1 n 1 


(5- 52a) 


•Mj U) 
n 


and 


,(n) _ 


V ?V " v c 

*- 

n times 


N 

I 


1 1 1 X 2 1 


N N 

l I 

i =1 


v n 


e . 


n 


i 3c. 3C. 

n i, i rt ... l 

12 n 


(5.52b) 


and with equation (5-17) the N -dimensional Taylor's expansion is 


r? ^ ^*( C 1 > • • • 9 C-jvr) 

tlz ,- z 2 V " f( 5j.5 2 . •••■%)+ X <V 5 1> sT 


i=l 


N N 


2 T I I (z -c..)(z -£ ) 

i=l j=l 11 J J 




35. 35. 
i J 


(5.53) 
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which, in the above-defined notation, is_ equivalent to : 

td) = tG) * l £ v< n) fG) 

I ni Cf 

(compare with eq. (5.11a) )• Then also 

h?> = ZG) + i £ [Jh)i <n) (h v< n) set) 

In • c, 

(cf. eq. ( 5 . lib ) ) , and these latter two equations determine the N 
sional Lagrange's expansion to any desired order. For example, to 

f(z) = fit) + ey(z) * fit) 

+ \ £ 2 [y(z)y(z)] : 7 V ? f(c) + 0 (e 3 ) 

where 

y(z) = y(^) + ey(z) • V ? y(c) + 0(e 2 ) 

= iiit) + £y(?) * y(c) + 0(e 2 ) 

Thus (5* 55a) becomes 

f(z) = f(t) + £y(t) * fit) + £ 2 ([y(t) * tit)] * V ? fit) 

+ \ [p(t)y(t)]: \\ tit)} + 0(e 3 ) 

Since it is easily verified that 

y(?)y(t) : V V f it) = tit) * V [y(t) * V f(t)] 

- [y(t) • V ? y(t)] • V ? f(t) 


(5.5^a) 


(5.5^) 

-dimen- 

order 


(5.55a) 


(5.55b) 


(5.55c) 


(5.56) 
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equation (5.55c) is seen to be equivalent to equation (5-20) to order e 2 . 


5»5 Inclusion of Additional Parameters in the N-Dimensional Lagrange 
Expansion, and Development of Perturbation-Expansion Scheme 

As in i 5.3 one may also include additional parameters as arguments 

~y _Jy 

of f , y , and z in equations (5-17) through (5-20) or in equations 
(5.33) and (5.34). 

Thus, equations (5-33) and (5-34) may he generalized to: 

If 


z = z(? ,e,a,3, . . . ) 


C + ey(z,a,3, . . . ) 


(5-57) 


then 


(5.58a) 


f ( z ,a, 3 , . . . ) = f(S,a,3,. 


) + l 

n=l 


£ 

n! 


3 n f[z(g»e 9 a,g a . . . ) , 


8s 


n 


Je=0 


where 


(5.58b) 


9 t [ z ( ^ <»£ s0c>3>»..’)a 


9e 


•i(T) 


3 a y( z,a,3, , 


_ti| v fs 1 J 1 f(z,q,B,..) 

- 1 . 


and 9 1 y(z 5 a 5 (3 5 ...)/9£ 1 


£(?)[**'* 


9 u(z a cx,3»...) 


3 e 


9* ^ y(z,06,3a*..) 

3 e 1 - J - 1 J 


(5.58C) 


The expansion of (5*57) may also he found from ( 5 . 58 a) hy taking 
f = p( z,a, 3 , . . . ) . In more explicit form for the first few terms, the 
corresponding generalization of equations (5.17) and (5*20) is: 
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If 


z = z(t,e . . . ) = t + ey(z,ot,3, . . • ) (5. 59a) 

then 

f(z 5 ot 3 3»***) — f(c»ot 9 3»...) £{y(Cj^535**-)*V^ f(£ 9 ot 9 35...)} 

+ [terms same as in eq. (5-20) 9 hut with f 

and y having arguments (c" 9 a,3 9 • . - ) ] (5.59b) 

and 

y(z,a,3,...) = y(t,a,g,. .. ) + e{y (c,a,3 , • • • ) • V^y (t>a,B , . . . ) } +••• (5.59c) 

Consider the special case of only on£ parameter a (as in eqs . 

( 5-1*0 ) ; that is 9 


f(z,a) 

= l « k 1 f k (z) 

k=l * 

(5.60a) 

P(z,a) 

oo 

r k-1 /-K 

= 1 a y k (z) 

( 5 . 60 b) 


k=l 

Then equations (5-59) become 

f(z,a) = l a k_1 f At) 
k=l K 

+ e{[ l a k_1 t -Vl[ l a k_1 + 0( e 2 ) (5.6la) 

k=l k C k=l k 

where 


z - £ + e{ J a k " L p,(?)} + 0(e 2 ) 

k=l 


( 5 • 6lb ) 
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In particular, for a = e 


f(z,a) = f(z,e) = f ( z ) + e f (z) + e 2 f (z) + ••• 

1 2 3 

y(z,a) = p(z,e) = ^(z) + e p 2 (z) + e 2 p 3 (z) + ••• 


(5.62) 


and equations (5*59) become (cf. eq. (5-20)): 

f(z,e) = f (?) + e f (?) + e 2 f (?) + 0(e 3 ) 

1 2 3 

+ e{[p (?)*V + e y C?)*V + •••][f ("?) + e f C?) + ***] 

1 ? 2 ? 1 2 

+ (1/2!) e 2 { [p i ( ? ) • + 0(e) ] [^ ("? ) *V C? ) + 0(e)] 

+ [^(C)-V ? ^ Cc) 3 -V ? f x (?)} + 0(e 3 ) 

where 

z = ? + efpj?) + e Z 2 (D + 0( e 2 ) } 

+ e 2 (u 1 (?).V c ^(?) + 0(e)} + 0( e 3 ) 
or 

f(z, e ) = ^(?) + e{f 2 (?) + ^(c).V f^?)} 

+ * 2 {f 3 (?) + V* ),v ? f 2 (c) + y 2 (?)-v ? f x (t) 

+ (1/2) p (?)-V [p (?)-V f (?)] + (l/2)[y (?)*V p (?)]*V f (?)} 
1 1 1 1 1 Q»1 

+ 0(e 3 ) (5.63a) 
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where 


/ -K p f -> 

z = ^ + e y n U) + e (y 0 U) + y, 


y x (t)} + 


This result (eqs. (5*63)) has been derived for use in the method of 
Chapter IV, to be shown in Chapter VI. 
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CHAPTER VI 


APPLICATION OF LAGRANGE-EXPANSI ON PERTURBATION SCHEME 
IN THE DIRECTIONAL-MEAN- FREE- PATH METHOD 


6 . 1 Introductory Remarks 

The Lagrange- expans ion perturbation scheme can be applied, purely 
formally, directly to the gain term in each of the directional equations 
of change (i.e., to the first term on the right side of each of eqs . 
(4.43a,b,c), or of eqs. ( 4. 973 t,b ,c ) in the special case of one-dimensional 
flow) . 

For validity of the scheme (sufficient convergence or asymptotic 

convergence), it is not necessary that e be small. In fact, e can be 

taken to be unity (so that L = u /© ), since the M arbitrary length, " 

a a 

L, has not previously been defined. Lagrange* s expansion is often written 
without a small parameter ( cf. Whittaker and Watson, 1927), which is equiva- 
lent to taking e = 1 in the forms given in Chapter V. It is only necessary 
that the product ey be "sufficiently small" in equation (4.4l) or in 
equation (4.95), that is, that ?-£ be sufficiently small. However, the 
parameter e can still be left in the equations for convenience in 
"ordering" of the terms, discussed further in § 6.3 below. (See also 
related discussion of " artificial parameters * 1 by Chang ( 1961 ), pp. 8 l 6 , 
820-825.) 

Each of the vector and scalar functions of z (denoted by ( )) 

in the right side of equations (4.43a,b,c) for general three-dimensional 
flow can be expanded by use of equations ( 5 . 63 ) , and subsequently treated 
by the scheme discussed below in § 6 . 3 . 

It may be instructive at this point to recall that the relation 
(cf. eq. (4.4l)) 

z = ? + ey(z,6 ,<jp, e) (6.1) 
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is a dimensionless vector equation with four components (cf. eqs . (k.2), 
(t.3), (ii.28), and (l;.4o)) representing 


Xi = Xj - (c os <p) 

x 2 = x 2 - (sin (p cos e)x^ ) (x 1 ,x 2 ,x 3 ,e , 4 >,t) 

x 3 = x 3 - (sin(P sin e) A* (xj ,x 2 ,x 3 ,9 ,<p,t ) 

T . ^ 
t “ t - / x 




) ( 6 . 2 ) 




Since the procedure in three configuration-space dimensions (with time 
as the fourth dimension) is directly analogous to that in the simplified 
one-dimensional- flow case, for simplicity and economy of space only the 
latter is treated in detail in §§ 6.2 and 6.3 below. In one- dimensional 
flow, only the components ^ and 5 of X are of significance, and 
y is therefore a function only of £ , £ , w, and e (but is independent 
of c, in steady flow) . 


6.2 Expansions for One-Dimensional Flow 

In the equations for one-dimensional flow (§ U.4), assume that each 
function of X and £> or of X> w, and e, can be expanded in a power 
series* in e: 


FU ,£) = Fjft) + e F 2 (t) + e 2 Fg(t) + 0(e 3 ) 
F(t,co,e) = F 1 (t,t») + e F 2 (e,u>) + e 2 F 3 (c,w) + 0(e 3 ) 


(6.3) 


* These expansions may be only asymptotically valid as e 0, and 
not "analytic” in e , but this does not preclude use of the symbol 

because the order symbol (e.g., 0(e 3 ) as e 0) is used (see, e.g., 

Erdelyi, 1956, or Martin, 196 Tb). 
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and each function of z, go, and e has the expansion 
F(z,oj,e) = F 1 (z,oj) + e F 2 (z,<o) + e 2 F 3 (z,id) + 0(e 3 ) 
Thus: 


p"(c jO),e) 
J(Z sO)»e ) 

P(t ,u),s) 
eCc,u),e) 

Q( C £>,£ ) 


+ e P 2 (C»w) + 0 (e 2 ) 

JjCCjto) + e J 2 (? 5 u) + 0(e 2 ) 

P^t.u) + e P 2 (t s u) + 0(e 2 ) 

E (?,to) + £ E ( £ ,u)) + 0(e 2 ) 
1 2 

Q 1 (?,w) + e Q 2 (c»o>) + 0(e 2 ) 


and 


(6.4) 


a 


> (6.5a) 




H(^,o),e) 
Z(^,U)»E ) 
m( c ,u,e) 
M(t.e) = 


= H^c^o) + e H 2 (c,co) + 0(e 2 ) 

= Z 1 (^,o)) + e Z 2 (c,w) + 0 (e 2 ) 

= ^(£,0)) + e y 2 (t»w) + 0 ( e 2 ) 

M x (c) + e M 2 (t) + 0(e 2 ) 


> (6.5b) 


J 


Since for one-dimensional flow the functions vary with only two components 

V 

of £ (or zj , as well as with go, the Lagrange- expansion perturbation 
scheme for 


z = C + £m(z,0),e) 


( 6 . 6 ) 


and for any function such as (6.4), takes the form (cf. eqs . (5*63)): 
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F = F(z,u),e) = + e[F 2 (t,u) + F^t.o))] 

+ 0( e 2) 

= F^ ,oj ) + e[F 2 (|',a)) + (v x ) x (£»<*>) ~~ F(?,o)) 
+ (yj^c.w) 0(e2 ) 

where (see eq* (4.69)): 


N _ — 3 J /9 03 . 

e l.p ( C,co,e) - STpVTaT " 


SJ^/Sco + e 3J 2 /3 w + 0(e 2 ) 


[H^ + e + 0(e 2 ) ] [3p ^ /3a) + e 3p 2 /3(o 


— 3 /3(o 


3J 2 /3u) 


3 /3co 


H, 


P 3p 9 /3cu N 

+ z 


H x Sp^Su) 35^/30) " \^H 1 Sp^StoM ^ 3^/30) 


= (u x ) 2 ( C ) + e(y 2 ) 2 (? ,a>) + 0(e 2 ) 

and 

* -*,± > _ 1 -1 H 2 (e,o>) 

V^\“’ e) --f = Hl (t.«.) C H^,*) * ° (e) 

= (p 4 ) 1 (t»u) + e(p 4 ) 2 (?,oo) + 0(e 2 ) 


Thus, in the expansion form of equation (6.7), for any function 


3J (c,o)) / 3p (^,w) 

^ / H, ( c. CO) — ^ 


3o) 


and 


(y 4 ) 1 (c,w) = - l/Hj ( C ,to) 
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(6.7) 


+ 0(e 2 ) ] 

(6.8a) 

+ 0(e 2 ) 

( 6 . 8t ) 

(6.9a) 
(6.9b) 
F(z,m,e) : 
(6.10a) 

(6.10b) 


I II 


I I II I 


i mil 



For equations that use collision model II , only the expansion forms 
(6.3) and (6.7) vith (6.10) are needed. Thus the directional equations 
of change for that case (eqs. (U.107)) "become: 


(!!!!L + !5_) * + !5_) +0( 

3^3(1)/ \3£ 4 3u) 


e 3 ) 


= H 


3p 


le 


3p 


1 e 3 uj 


- H 


1 


1 du 


+ e 


( H .«[^ + Vi a JrS^r)] 


3p 


le 
3 oo 


[ 


H 


2e 


3H_ 3H. ‘ 

* irf * '"-’i -5^7 


- H 


3p 2 


9 Pl 


1 3o) H 2 3 u 


+ 0(e 2 ) 


(6.11a) 


/9 2 Jl 

9 2 

j. 

If) 

/3 2 

. f 2 

j 2 

3 2 P 2 

\ 3C 4 9u) 

3C 

1 3<*)/ 

+ t 1 

Vh, 

^30) + 

3^ 4 3b) 

= H 

3 J 

le 

- H i 

j 

!„ 

rsj 0 

2e 

le 

3o) 

3o) £ j 

| .e 

3d) 


) 


+ 0(e 3 ) 


‘*.>1 ♦W.i f-(^) 


3J 


le 

3(A) 


t 


3H 


011 ie 

H 2e + ^ y l t + ^4^1 


3H 


1 e 


n L 


3 J 


- H 


1 9o) 


2 h ■— L 1 

H 2 9u) j 


+ 0(e 2 ) 


(6.11b) 


/ 3 2 E 1 3 2 * l \ 9 (* 2 *2 * Z 

\ + ) + e z I : + I + 0(e 3 ) 

\3^ 4 3o) 3C 1 3o)/ \3^ 4 3a) 3^30)/ 


3E 


= H 


le 


le 3o) 


( r 9E o ~ 

/ 3E \ 


/ 9E \" 

| H le 3o) + ^ y l t 3? 1 

( le ) 

+ 3C 4 ^ 

16 ) 

V 30)/ 

V 94) / 


3E 


le 


3(d 


9H le 3H l 

H 2 e + Tc7 + (y 4^1 ~n 


le 

4 


3E, 3E, 

H 1 3u) “ H 2 So) ^ + °( e2 ^ 


(6.11c) 
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and 


3E X 3E 2 

3 co 9oo 


+ 0(e 2 ) 


1 

2w 2 



+ e 


1 

2 u) 2 



+ 0( e 2 ) 


(6. lid) 


* “T 

For other collision models, evaluation of functions of z and to ’ 
in equations (4.97) is required. Since w' is given in general by 
( 4 . 106 b) , it is convenient to write the expansion 


(t,w,e) = ajj(t,t*)) 


e w. 


,(£»“) + 0(e 2 ) 


(6.12) 


and then to write 
F(t»w,e) = 
F'e F(z,o)' ,e) = 


(cf. eqs . (6.3) and (6.7)): 

F^ (^ ,oj) + e F^ (^,o>) + 0(e 2 ) 

F^ (^,o)' ) + e[F 2 (^,o)')+ (u 1 ) 1 (c»a)) 


9 




Fj (£.«•) 


(6.13a) 


+ (y ) (c,») ^-F^c.u,’)] + 0(e 2 ) 


(6.13b) 


and, writing each F.(^,(o') in terms of a Taylor's series about 

J 


)' = a>jU,w) : 

3F,(t # w{) 

f .( tv) = F j(t,a>;) + (w , -o ) ;) J 9u) , ■ 

(oj'-o)! ) 3 2 F ( t, ,ouI ) 

j. : —u 1 


21 


9w{2 


one obtains (6.13b) in the form (with use of (6.12)): 

aF^c.uJ) 


F' = F(z,aj',e) = F^CjCoJ) + e[<u^(c,w) 


(Ul)l(t,o)) 3 ^-F : (t^{) + (y 4 ) 1 (t,«o) 3 ^- Fi(c,03’ )] + 0 ( e 2 ) 


9toJ 

_9_ 


+ F (?,«’) 


(6.13c) 
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As a corollary of (6.13c) one may also write 


(SL\ 3F(z.to' « e ) = 3F i ( ^ >ai P + e f 9F 2 ^ >a) i 

V 3 cj / 3 co f 3 co ^ 3 co ^ 


+ u' 2 (t»u) 




9cu 


» 2 


(6.13d) 


_ /3F (?,a ,{)\ ^ . / 3F 

+ ^ (-‘-nr)* d: (-V)J* 0U2) 

Note that, for any function G(l;,u)’) in these expressions, one may write 


w( *-;> . I L i ) 

/ 3u>J Of ,to) \ 

\ 3^5 / 


9o) ’ 


(6.14) 


For collision model Ilia , the terms in (6.12), from (4.108), are 
ITT [l-A^c.u)] I [1-A (|,m)] 2 + ^r 2 


a)j(^sO)) 


1 


(6.15a) 




l-A^t.u) 




(6.15b) 


where the expansion 


A(^ 5 co 9 e) = A^^joj) + e A^(^ 9 oj) + 0(e 2 ) 


(6.16) 


has been used. The expression for Z _ from (4.109) 9 combined with 

JLX.L a 

(4.108) , is 


Z IIIa^ ,a3,e ^ 


g) 


u>[l-A(t,u),e) ] l-(D 2 A(^,a),e) 


(6.17) 
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With the expansion of (6. IT) given in (6.5^), one finds 


Z 1 (t,u) = 


Z 2 (t,u) 


U)J(C,U)) 

wCl-A^^jO))] l-u) 2 A i (c,w) | 


= Z 


-,*4# 

L »< 


(C,oo) 

, 0 )) 


A 2 (? ,w) 

l-A^t.o)) 


a ) 2 A 2 (c,w) 1 

1-u) 2 A i (^,o)} J 


( 6 . 18 a) 


or, with use of (6.15), 


Z U,w)A (e,w) 

z,(f,») = r 

' l-Ajtc.u) 


m 2 ( 1-Aj ) 


j 2 (1-A 1 ) +1-0) 2 1-o) 2 A 1 


( 6 . 18 b) 


From the definition in (4.110) one finds, upon using the appropriate 
expansions , 


A(? ,o),e) 


+ 


= 0 + £ 

9p 1 ( C,m)/9m 

-j^inr^yto,) 

BJ^C.l) 

_ 3 / 3 co_ 

_p' 1 (f,i)JL jj(t,i) 

3?1 

(y 4 ) 1 (t,w) 

3J^ ( C >l) 

(p ) 1 C »co) 9p 1 ( C ,l) 


Jjtf.i) 


Pjd.i) 


(y lt ) 1 (^,w) 

3p 1 (t,l) "| 

f 0(e 2 ) 


P,(t,D 

3? 4 J 



( 6 . 19 a) 


= A l (t,w) + e A 2 (^,w) + 0(e 2 ) 


( 6 . 19 b) 


The vanishing of the first term in (6.19a), its consequences, and a remedy 
are discussed below in § 6.3. 



0) 


In the following, denote by superscript ( )" a function of £ and 

with (o replaced by ; that is 


F" = F(t,<oJ) 


( 6 . 20 ) 


Then, in terms of the above expansion forms ((6.13) with (6.10)), 
equations (U.97) become: 


/3 2 p 1 92 Ji\ / 92 P 2 9 2 J 2 \ 

t V3C„9« + HjSu/ * E Ui t 3to + a C] a«j/ + o!e * 


<V S. - H 
1 1 




• r ^2 

ttM ^ 

11 [ 9w{ 


+ (oi 


3 V; 


2 amp" 


3 2 pV 


9 2 p" 


+ ^1^1 + ^4^1 9?,. 3m,' 


4 1 J 


t (o)f *L 

Z 1 9wJ 


9H" 

h ; + “^ 


9H" 9H" ' 

^ P l^l 9c, + ^ y 4 ^ 1 3c, 


(0)- „ iJl 

+ Z 2 H 1 3a,} 


- H 


1 3o) 


9p 


- H 


1 


2 9(o 


+ 0( e 2 ) 


(6.21a) 


/ 92J, 3 2 PA 2 / 92 J 2 92 P 2 \ 

'\3C 4 3(o + 3C 1 3(o/ +£ \3C 4 3 (o + 9^ 1 3to/ 


+ 0(e 3 ) 


r i )- 3 j’; 9j ! 

- Vl’ip - H l — +e 

9 2 J" 9 2 J" 

+ (yJ, ~ *TT + (p 4 )i 1 


f(i) SH „ [!fi . , _ei 

Z 1 H 1 [3(o; “2 9(0**- 


11 9? 1 9(oJ 


3 ? l 3 “Ij 


+ d>z 

Z 1 9(0* 


3H" 

H" + to' — r- 
2 2 3o)J 


3H" 9H” 

^ y i^i 9c7 + ^ y 4 ^ 1 FET 


(D- 3J i 

+ ( )z X irr - 


9J„ 


2“1 3(o{ H 19 (o H 2 


9^) 

9(0 J 


+ 0(e 2 ) 


(6.21b) 
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1 1 am iai ii i i i ■ 


/ 3 2 E 3 2 Q x \ /B 2 ®! 3 2 Q x \ , 

Vfijjlio + 3 T^) + £ \3C 4 3o) + 3^W + ° (£ } 

3E" 3Ej ( (2) _ fSEj 

- Vl' 3 ^ - H 1 — + E ( Vl [l^ + 


3 2 E' 1 ' 


2 3to{ 2 


3 2 e’; 


3 ^E ?i 

+ ^l\ SCjSioJ + 3? 4 3 oi* 


+ ^Z. 


3E" f 3H" 

i aSf ^ 


3H" 3H” 

+ (p, ), TT“ + (p ) TT" 


l'l 3? 


Vi h l 


( s 3E" 3E, 

. IV TTtt L _ II i 

Z 2 H 1 3oa» H I 3 oj 


-h 2 ^I +o< £2 ) 


(6.21c) 


and 

3E 


3E 


- — + e t + 0(e 2 ) = 


3oj 


3oj 


i !!l 

2 03 2 3 co 


+ £ 


1 ^2 

2cjO Z 3(a) 


(4.97) for one-dimensional flow, become: 



o, « 

*3p 2 (t.l) 

9J 2 (C,D] 

n ' 3? i J 

"T £ 

9? 4 

J 

(■SJjtCvL) 3P,(tA>" 

-1 . r~- 

"3J 2 (t,l) 

3P 2 (C,1)] 

[ 9 «4 ' 9 S J 

T e 

9C 4 

9 S J 

fsljfc.l) dQjtc.l) 

J r- 

'3E 2 (c,1) 

3Q 2 (t,l)] 

l 9 «3 ‘ 9 ^1 J 

■ b 

^4 

9?1 J 


+ 0(e 2 ) 

(6.21a) 

e form of equations 

+ 0( e 2 ) = 0 

(6.22a) 

+ 0(e 2 ) = 0 

(6.22b) 

o 

II 

CM 

W 

o 

+ 

(6.22c) 


( 2 ) ( 3 ) 

The expansions of the functions ij; » and ip defined by 

equations (4.99) for use in the implicit closure equation (4.70a), 
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representing equations (4.69) » are: 


/ 2) = + e ^ 2) + o(e 2 ) 


(3P 1 /3o))(ap- 1 /9to) 


1‘ + e 


( 3^/30 )Y 


i 3p i 


/dhl 


3p 2 /3o) 




/3oj apj/au 


8J 9 /3w I 

2 I + 0(e 2 ) 

3J 1 / 3u J 


> (3) = ^ 3) + e ^ 3) + 0( e 2 ) 


2w 2 ( BQ 1 /3o)) ( 3p j / So)) 2, 

(a'ji/Su)) 3 


1 + e 


3Q 2 /3(jo 

3^/30) 


3p /3cc 

+ 2 -==-7 — 

3p ^ /3oj 


- 3 


3J 2 /3w ' 
3^/3^ 


+ 0(e 2 ) 


and 


i|»l + e ip + 0(e 2 ) 


(3F 1 /3ca) 


2a) 2 (3Q 1 /3aj)(3J 1 /3a)) 


+ 0(e) 


6.3 A Truncation Scheme by Introduction of an Additional 

Artificial Parameter, Illustrated in One-Dimensional Flow 


As alluded to briefly in § 6.1, an asymptotic expansion in an 
"artificial parameter" e may be used in a scheme of successive 


(6.23a) 


(6.23b) 


(6.23c) 
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approximations in different ways: (a) In a most conventional manner, one 

may substitute the expansions into the problem and let e 0 successively 
to obtain equations and boundary conditions for the respective terms of the 
expansions; or (b) one may let e = 1 (thus defining the "arbitrary 
length" in the definition of e) and suitably truncate the expansions 
to obtain the solution at each step of the approximation scheme. (An 
important illustration of one way of using the latter approach is the 
Chapman-Enskog procedure of kinetic theory; cf. Chapman and Cowling, 

1961, Chap. T; or Grad, 1958, pp. 253, 259-266.) 

To circumvent nonuniformities in the mathematical solution in the 
first approach, it is often necessary to prevent the system of. equations 
from degenerating (e.g., because of elimination of higher -order derivatives) 
in the limit as e 0 by some means such as introducing suitable trans- 
formations of variables (see, e.g.. Van Dyke, 196 ^; or Martin, 196 Tb). 

For similar purposes, the "suitable truncation" in the second approach 
above can be taken to be whatever appropriate means prevents the degeneracy . 

The second approach, including the "suitable truncation", can be 
formalized by: leaving the parameter e in the equations, to aid in 

ordering the terms, and by introducing another artificial parameter e T 
(which in the final analysis will also have the value unity) at appropriate 
places in the equations; then letting e -> -0 to obtain an asymptotic 
solution and at the same time letting e 1 -* 0 such that, for example, 
e , /e = 0(l) as e 0. This procedure and its motivation are discussed 
and illustrated in the following treatment of the equations from § 6.2. 

Consider first equations (6.3) to (6.23) in which e is not yet 
defined. Suppose we used these equations as they stand and let e -* 0 
to obtain an asymptotic solution (i.e., the first approach mentioned above). 
If we used either collision model I or collision model Ilia (with Z 1 
and Z 2 found from eqs . ( 6 . 15 ) to (6.19)) in equations (6.21), we would 
find (for the steady- flow shock-structure problem of Chapter VII ) that 
the first-order term on the right side of each of equations (6.21a,b,c) 
would vanish identically; and with the resulting set of equations, in 
the limit as e -* 0, it would not be possible to obtain variations of 
the flow variables with the independent variable x because there would 
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be no way to determine a length scale in this formulation. (Any scale 
factor on £ would cancel out.) The only possible solution that would 
satisfy both upstream and downstream boundary conditions would be discon- 
tinuous in x, and a continuous shock-structure solution could not 
therefore be obtained. 

On the other hand, the limit e -*■ 0 in equations (6.1l) for collision 
model II would cause the equations to degenerate to lower-order differential 
equations. The first approximation would then correspond to local-transla- 
tional-equilibrium flow (with a local-Maxwellian velocity distribution) 5 
goverened by the Euler equations of inviscid flow, as found in the Chapman- 
Enskog procedure, with only a discontinuous shock solution possible for 
collision model II (as for I and Ilia), 

The nonuniform limits described above can then be prevented for 
collision models II and Ilia as follows: 

For collision model II : Arbitrarily divide the left side of each of 
equations (6.11a,b,c) by and require e ! /e - 1 as e -+ 0. (This 

step is motivated also by the fact that the left-side and right-side terms 
of the Boltzmann equation are of the same order very near and within a 
shock; cf. Grad, i 960 , p. 117*) By this procedure' we obtain the first- 
order equations as e -* 0 with e'/e “ Is 


9 2 p 9 2 J 

L + — 

dZ 3 co 9 c; 9 w 
4 1 


9 2 J 9 2 P 

1 + 1 


3 £ ^ 9 (0 9 £ ^ 3o) 


= H 


= H 


9p 


ie 


1 e 3 oj 


3 J 


ie 


la 3 to 


- H 


- H 


3Pj 

1 3(0 

5 

1 9(0 


9 2 E 


3 Z Q, 


9 C 9(0 9 £ 9(0 

4 1 


= H 


3E 


ie 


ie 3(o 


- H, 


9E 
1 9(0 


3E 
3 to 


1 ffi 

2 uj 2 9w 


(6.24a) 


(6.24b) 


(6.24c) 


(6.24d) 


where, from equations (4.100), (4.101), and the expansions (6.5), 


03 = /5 7 6 Mj 03 


(6.25) 
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and 



3o) 



3(o 



3co 



( 6 . 26 a) 

( 6 . 26 b) 

( 6 . 26 c) 


where = M^(^) is given by (4.105b) with subscript 1 on all functions. 

For collision model Ilia : Divide the left side of each of equations 

(6.21a,b,c) by e' and at the same time divide the right side of equation 
(6.19a) by e' (to prevent degeneracy of the collision term); then 
require e’/e = 1 as e 0. By this procedure the first order equations 


as e -* 0 

for 

collision model Ilia 

are 

obtained: 




32j i 


z i h i 

3p\' 

- H i 

3p x 

(6.27a) 

3^3(0 

+ 

3^30) 


3u)J 

3(1) 

3 2 7 


3 2 P 


Z H" 
1 1 

3J" 

- H 

1 

3?i 

(6.27b) 

3 £^3go 

+ 

3^3(0 


3(o| 

3d) 

3 2^ 




VI 

3E" 

- H 1 

3E 

(6.27c) 

3 £^3u) 

+ 

3 ^ ^ 3(i) 


3io’ 

3d) 



3Ei 
3 co 

= 

1 

2m 2 

3?1 

3o) 



(6.27d) 

vhere co 1 = 


(£,u>) 

is 

given by ( 6 . 15 a) 

, Zj = Z^t.oo) 

is given by 


(6.l8a) , in which , from equations (6.19)» becomes 
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I 


A 1 (t > 4 ) ) 


/3p 1 /3(o\ 

"j x (t.n" 


Vi 

3Jj(c.D (y 4 )! aJ^c.i) 

\3 /3o)/ 

p^t.n 


ji(t.n 



(p x x 3p 1 (t»l) 


(u lt ) 1 3p 1 (c»i) 
3? 4 


( 6 . 28 ) 


and where (y^^ and (y ^ are given by (6.10). 

For either collision model II or collision model Ilia, the systems 
of equations (6.2U) or (6.27) are closed by 


(2 ) 

= vt^) (6.29a) 

where 

= (*[ Z) ) 2 /*[ 3) (6.29b) 


and. where ^ and are the first terms in (6.23). Also 

imposed are the conditions from macroscopic conservation of mass, momentum, 
and energy ( from (6.22)): 


3 

9 ^4 

pjC?* 1 ) + 

3 

Jjtc.l) 

= 0 

(6.30a) 

3 

9 ^4 

+ 

3 

Pjft.l) 

= 0 

(6.30b) 

3 

9 ^4 

1(1,1) + 

3 

9 ^i 

Qj(tA) 

= 0 

(6.30c) 


For a given problem in one-dimensional flow, these equations would of course 
have additional initial and boundary conditions in the X domain, as well 
as the conditions at a) = -1 (cf. § U.4): 






Q^C.-l) = 0 (6.31) 


lUl 



Recall that w = cos <p , a) 1 = cos <p f , and that 

5 x i /L and = u^t/L 

where (with e' = e = l) L = u a /0 a . In a given problem, and 0^ 

would be suitably determined as discussed in §§ 4.1.2 and 4.3.3. For 
example if equation (4.85a) is used, then from (4.92a): = Hj(£) = p"^^,!). 

Recall also now. that the above procedure has simply constituted a formal 
means for truncating the equations to obtain a hondegenerate first approxi- 
mation, and that actually s' = e = 1. Just as the initial introduction 
of the arbitrary length L was a convenient artifice, so also is the 
introduction of e’, to obtain a first-order system of equations that is 
nondegenerate. We take it to be significant that the form of the above 
first-order equations is very similar to the original unexpanded equations 
in § 4.4.2. The expansion process and the truncation scheme have simplified 
the first term on the right side of the directional equations of change in 
the first approximation and have provided a means for obtaining higher 
approximations. The utility of higher order solutions by this scheme may 
be limited in a manner similar to that in which the Burnett equations , and 
higher order equations, are limited in the Chapman-Enskog method (see 
Grad 1963, p. 149.; or Vincenti and Kruger 1965, pp. 4l6-4l8). However, 
in view of the very recent findings by Khosla (1967) that nonuniformities 
in the Chapman-Enskog expansion are artificial, and are removable by use 
of Lighthill's technique (1949, 196l) with the extension by Kuo (magnifica- 
tion of the independent variable; see Tsien 1956) , the present higher- 
order solutions may be expected to be valid (see further discussion below). 

Although this method is certainly not equivalent to the Chapman- 
Enskog procedure, there are similarities that can aid understanding of 
some aspects of this method with sufficient prior understanding of the 
Chapman-Enskog procedure. (Substantial understanding of the Chapman- 
Enskog procedure can be gleaned from the discussions by Chapman and 
Cowling, 1961 ; Grad, 1958, 1963; or Vincenti and Kruger, 1965; on the 
respective pages of those works referred to earlier in this section.) 
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Conversely, because of similarities, understanding of the present method 
may also contribute to understanding some aspects of the Chapman-Enskog 
procedure. The most striking similarity between the present method and 
the Chapman-Enskog procedure is the occurrence of successively higher 
derivatives in the higher approximations, which are seen to occur here 
directly by introduction of the Lagrange- expans ion perturbation scheme, 
developed in Chapter V above, into the formulation of the equations of the 
directional-mean-free-path method. (However, Khosla, 1967, finds that use 
of the PLK method eliminates the higher derivatives in the Chapman-Enskog 
expansion. Similarly, in the higher order terms in the present method, 
the higher derivatives are only those of previously-determined lower-order- 
approximation functions; so there is no increased number of boundary 
conditions required. ) 

The first approximation in the Chapman-Enskog procedure (resulting 
in the Euler equations of inviscid flow) would correspond to the method 
described here in the first approximation ±f_ we did not introduce e T = e 
into the equations, as described above. The second approximation in the 
Chapman-Enskog procedure (first translational- nonequilibrium correction 
to the local Maxwellian distribution function, or to the Euler equations) 
yields the Navier-Stokes equations. The mathematical mechanism by which 
the higher order approximations in the Chapman-Enskog method eliminate the 
degeneracy, and hence the discontinuous solutions of the first (Euler) 
approximation, is described by Grad (i 960 ) , p. 119- 

In the present procedure, the use of the truncation scheme with 
e ? = e gives a nondegenerate solution in the first approximation, and 
so is fundamentally different from the Chapman-Enskog procedure in that 
respect. The first nonequilibrium approximation here is not a correction 
to local translational equilibrium. It is suggested therefore that it 
may be capable of describing the flow approximately, but realistically, 
for arbitrarily large deviations from local translational equilibrium, 
which the Navier-Stokes equations can not. 


CHAPTER VTI 


AN ILLUSTRATION: 

PROBLEM OF SHOCK-WAVE STRUCTURE 


T . 1 Introductory Remarks 


As an illustration of the method developed above, consider the 
equations in § 6.3, for the first approximation in the direct ton al~me an - 
free-path method, applied to one-dimensional steady flow through a normal 
shock wave in a perfect (sufficiently rarefied) gas composed of neutral 
spherically-symmetric monatomic molecules. 

Since this chapter is not intended as a thorough study of the shock- 
structure problem, but only as an illustration of the method, the reader 
is referred to works such as: Lighthill (1956), Hayes (i960), Vincenti 

and Kruger (1965) 5 Sherman and Talbot (i960), and Talbot (1962) for 
comprehensive treatments and further references on shock-wave structure. 

The conditions at upstream infinity (x x = -00) are denoted by sub- 
script a, those at downstream infinity (x 1 = +») by subscript 3. 

Since there is only one component of particular interest in one-dimensional 
steady flow, we drop the subscript 1 corresponding to the direction x 1 
and denote the dimensionless coordinate (with e = l) as 

C = (0 /u )x (7-1) 

a a 


where u is the macroscopic average velocity in the +x direction, and 

0 is the upstream collision frequency, which could be represented, for 
a 

example, by equation (2.20b) with (2.23d). Thus 


u 

a 


m 

(8n) 1/2 a 2 p 

a 



(7.2a) 


lkk 



for use in equation ( 7 *l) if* desired, where m is the mass of a molecule, 
a is an fT effective molecular diameter" (see § 2 . 2.2 above), y( = 5 / 3 ) 
is the ratio of specific heats, and M^ is the upstream Mach number. 

It is also convenient, for comparison with other results, to find the 
results simply in terms of x/A^ (e.g., see Vincenti and Kruger, 1965 , 
sec. X-9)* Thus from (2.21) and (2.23d): 

u /0 = (l/2)n 1/2 A (y/ 2) l/2 M ( 7 - 2 b) 

u u cx ct 

where y = 5 / 3 , so that ( 7 * 1 ) may be written 

f - = “ r - ( /I M <,> 5 ' ( t . 3 ) 

a 


where is the upstream mean free path relative to the mean mass motion 

(cf. Chap. II). 

In the following sections of this chapter, use is made of equation 
(4.85b), with (4.92a) and ( 6 . 5 ), to represent the function as 


Hj = HjU) = H^O/K^T) 


(7. to) 


where 


Kj(T) 




HjC?) = P^C.D 


(7. *rt>) 


and where q is a specified temperature- dependent viscosity coefficient. 
It is also convenient to define 

? s - / [l/K.(T)]dc ( 7 . to) 

C* 

where is some value of £ , to be defined later. Directional- 

average collision models II and Ilia are considered, and a numerical 
solution scheme is outlined for collision model II. 



7.2 Equations and Boundary Conditions 


For simplicity, in the following we denote partial derivatives with 
respect to u) hy subscript a). Since as represented by equations 

( 7.^)3 does not depend on w, equations ( 6 . 2 U) for collision model II 
are written for steady flow as: 


K — (J ) = 

1 9 C 1(J ° 


- - 2 = (j ) 
3C la) 




- 2= (p ) 

3? lui' 


K 1 Tr (2ti J = 


- If 


2E = iy P, 

1 (jl) 0) 1(a) 


h [t - p", 1 

1 K leo) K l(o 

(7-5a) 


(7.5b) 

- 25 lJ 

(7.5c) 


(7.5*) 


where 


and E are given by equations (6.26) with (6.25). 
16(0 


■ lew 5 l6(o 

Note now that equations (7*5a,b,c) are directly integrable to obtain 


K. 


3 

■ 9C 

(J 1 } = 

3 

- 9T 

(J i> 

= H i[p- ie -pi ] 

(7 • 5e) 

3 

(P_) = 

3 

(p. ) 

= H [ J - J ] 

(7-5f) 

- 9? 

1 

* ar 

1 

l ie 1 

3 

H 

( 2 Qj) = 

3 

" ar 

(2Q x ) 

= H 2 [ 2 E ie - 2 E 2 ] 

(7 . 5g) 

Pie 

, J 

9 le 5 

and 

E 

le 

are given by equations (U. 102 ) with 

o/p « 


replaced by and subscript 1 on M and of (cf. eq. (6.26)). 

For collision model III equations (6.27) become, for steady flow, 

] ( 7 . 6 a) 


K f- (J, ) = 

1 d£ 1 (a) 


!r ( 7„> 


H 1 [Z 1 e'iai’ - ^ 0 ," 
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(7.6b) 


K 


, hr (P. ) = - 

1 3 1, lu) 


3 /- \ 
w - (p i > 

3Q 1(0 


= H 


JZ.J" , - J ] 

1 1 l03 1 l0) 


K i k - 


J ‘ h i t z i 2E r«>' - 2 e i«.i 


2E la) = <l/» 2 )P lal 


(T.6c) 


(T.6d) 


where oj’ = co ^ ( 5 ) is given by (6.15a) and Z x = Z 1 (^,o)) is given by 

(6.l8a), where A^^ is given by (6.28) (which is considerably simplified 
by conditions given below). 

With equation (7-4), either equations (7. 5a) to (7-5d) or equations 
(7*6) provide four differential equations for the five unknown functions 
of £ and co : P x * » P 1 > Eq » and Q x * r ^ cie needed fifth equation 

(the "closure equation") is provided by equations ( 6 . 29 ), obtained from 
equations (4.70) (which represent (4.69)), in the form: 


'1 2) 5 P l„ P»u/ (J lJ 2 - T< *1> 


where 


#, = 


i (2)>2 /+! 3) 5 (p j 2 / 2 “ 2 v. j i. 


(7.7a) 


(7-7b) 


A sufficient number of boundary conditions on the dependent variables 
in co and £ must be provided. 

At co = -1 ( all c ) , the functions p’ 1 , «J 1 , P x , E 1 , and Q 1 
vanish identically (eqs. (6.31)). 

At co = +1 ( all £ ) 3 equations (6.30), with use of equations (4.92) 
and (4.60), give, in the first approximation 




( l/p a u a )( p u) 


(-“,!) 


constant 


(7.8a) 


lU7 


1 



I 


Pj(c.l) 


2Q 1 (c.D 


(l/p u Z )(pu 2 + p - t) = constant 
a a 

P.(-oo,l) = 1 + 3/5 M 2 

1 /a 

(2/p u 3 )[pu(h + 7 T u 2 ) + q. - ut] = constant 
a a ^ 


= 2Q r (-oo,l) = 1 + 3/k 2 


(7.8b) 


(7.8c) 


Note that the usual forms of the exact one-dimensional steady-flow conserva- 
tion equations are contained in equations (7*8a,b,c) and that their 
evaluation at x = ±«> , where q and x vanish, gives the well-known 

Rankine-Hugoniot equations : 


P u 
a a 


p u + p 
a a a 


h + \ u 2 
a 2 a 


p e u e 


p e u e + P 3 


h e + 2 u 3 


(7-9) 


For the perfect monatomic gas, with specific heat ratio of 5/3, one can 
obtain from equations (7-9) (see, e.g., Liepmann and Roshko, 1957, pp. 56- 
59): 


3+M 2 

a 

5 M 2 -l 

a 


(7.10a) 


V p 


a 


4 M 2 
a 

M 2 +3 

a 


1 + 


It 

5/2 

( 5/6)M 2 


(7.10b) 


At C = to ° (all a)) , equations ( 4 . 102 ) 
p", J, P, E, and Q in terms of P/P a 9 M, 


give the conditions on 
and Z) = Mw , where. 
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at £ = -°° : M = M 


p/p = 1 

a 


at £ = +°° : M = M 


6 5 


given by (7- 10a) 


p / p a = Pg/P a > given by (7-10b) 


> (7.11) 


J 


With conditions (7.8a) and equation (7.4b) we now note the simplified 
results from (4.105b) and (6.28): 


M 


1 



and 


A 


1 


1 




3/2 

2E 1 (C,1)H 1 (C)-1 


1/2 


(7.12) 


(7.13) 


for use, respectively, in (7-5) and (7*6). Equations (4.102) and (7.11) 
are used in the first approximation simply by putting subscript 1 on 
H, M, and ol (cf. eqs . (6.25) and (6.2 6)). 


7.3 Description of a Solution Procedure for Collision Model II 


Since 


P x ( C»w) 


GQ 



d u 


and 


Ej (?»«*>) 


E !0, d “ 


(T.l^a) 




lh9 



we can regard the set of ten equations including: (7-l4a,b), (7. 5a) to 

(7*5g)» and (7.7a), as a complete set to determine the ten dependent 
variables : 


lco 


1 03 


1(0 


1(0 


* 1(0 


Pi , 


E 




(7.15) 


as functions of £ and to (with c"( C ) » H^e), and Mj(c) given by 
(7.4b) and (7-12)). (To avoid, as much as possible, numerical integrations 
over to except for equations (7-1*0 » all of equations (7-5) can be treated 
as independent.) For each value of to, equations (7-5) can, in fact, 
then be regarded as a set of seven ordinary differential equations in the 
£ direction, with (o as a parameter, and with the needed values of the 
three quantities p^^^to) , p” 1 (C^)) ? and E^^oj) at each step in £ 
calculated respectively from (7-7), (7-l4a), and (7-l4b) for each io. 

Thus, the following numerical solution procedure can be used: 

(a) First, for a large number of values of b, compute and store, 

versus b : and ip from equations (4.6 9 ) through (4.78). 

The computer program can be designed to interpolate between values in this 

( 2 ) 

table to obtain a value of tJj corresponding to any specified value of 

ip in the range calculated. This provides the functional relationship 
(4.70a) . 

(b) Define £ = 0 as the location where p = (l/2)(p + p ). The 

a p 

numerical integration can be started at a location in £ where the state 
of the gas is very nearly that at either upstream infinity (£ = - co ) or 
downstream infinity ( £ = +°°) . The choice appears arbitrary. However, 
it has previously been found in the Navier-Stokes solution of shock 
structure that it is easier to integrate from the downstream state to the 
upstream state rather than vice versa (Gilbarg and Paolucci, 1953) 3 
because of a saddle point in the temperature-velocity phase plane. 

Therefore the integration in the present approach will be started at 
£ = where the state of the gas has a specified arbitrary small devia- 
tion from the downstream state. It is presumed that the values of the ten 
dependent variables in (7-15) can be approximately determined for all co 


150 



from -1 to 1 at E = C* from the known conditions at downstream 
infinity (see § 7*3.2 below). The values of H ^ and M 1 are then also 
known from (7.4b) and (7-12 ), and the local- equilibrium values of all the 
dependent variables can be computed, for subsequent use in equations (7«5). 
Then it is convenient to define the new independent variable 


C = C* - C (7.16) 

so that ( cf. eq. (7.4c)): 

dC/d? = - dC/d? = K X (T) (7.17a) 

K = C = 0 at £ = C* (starting point for (7.17b) 

integration) 

and 


it K 

? = C* at ? = 0 (where P = i (p a +p g ) ) (T.17c) 

and the integration proceeds upstream with increasing positive values of 
X ( as £ decreases from £*) . 

(c) Divide the range -1 < to < 1 into an odd number, £ , of 
intervals ( Ato = 2/£), such that 

a). = -1 + J Au) = -1 + 2 j / £ , J = 0, 1, 2, ... , £ (7.18) 

J 

Calculations are to be made at each ok , and the specification of £ as 
an odd number insures that no calculations are made at exactly oj = 0 
(where there would be numerical difficulties introduced by dividing zero 
by zero) . 

(d) Starting at = 0 with the known values of the variables in 

( 7 . 15 ), increase £ by A£ . At ^ = 0 + Ac 5 and at each of the to. 

_ J 

in (7.l8) 3 calculate J laj , P lw , ^lo) 3 ^1 3 ^1 5 ^1 f rom equations 
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(7*5) and from values at z = 0, using a standard technique for ordinary 

differential equations. With each value of P laj calculated, evaluate 

E 10J from (7*5<l)‘ Also, with known J 1(jJ , P la) , and Q la) at each to , 

( 2 ) 

calculate ijn from (7*Tb), then find the corresponding from the 

1 (2) 1 

previously- computed relationship between and (cf. eq. (7- 7a) 

and step (a) above), and determine iT lw from the definition in (7.7a). 

Finally, determine and E 1 at each to by the integrations indicated 

in (7*1*0 • Then, again, and are known from ( 7 • ^b ) and (7-12), 

and the local-equilibrium values of all dependent variables are computed 

for subsequent use in equations (7-5). 

(e) At each Z , and at to = 1, evaluate all flow quantities of 
interest (to first order) from equations (^.93)* Evaluate Z using 
(7 • 17a, b ) . 

(f) Increase Z again by a small increment, and repeat the procedure. 
Continue the integration until all the quantities no longer vary with 
increasing ~Z > so that the upstream equilibrium state has been reached. 

(g) Where pT (^,l) = ^ (l + p /p ), evaluate (eq. (7.17c)), 

1 d- p Ot 

so that Z is known versus Z from (7.16) and x/X is known from 

a 

(7.3). 

The values of the dependent variables at each u> at z ~ 0, for 
starting the numerical integrations, can be evaluated from integrations of 
the Chapman-Enskog distribution function for a very small deviation from 
the downstream state. This procedure is simple and straight forward, 
either in terms of the BGK model or for Maxwell molecules, in combination 
with the method of Gilbarg and Paolucci (1953) for evaluating flow 
variables near the downstream state. It is hoped in the future to obtain 
numerical results for the shock-structure problem by the procedure outlined 
above. 
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CHAPTER VIII 


CONCLUDING REMARKS 


The most important results of this study are considered to he: 

(a) the development of concepts and equations for the directional level 
of description; (b) the introduction of the directional-mean-free-path 
approximation for the collision integrals in the directional equations of 
change; (c) the development of new forms of vector (higher dimensional) 
generalizations of Lagrange f s expansion, and a perturb at ion- exp ansi on 
scheme based on those generalizations; and (d) the outlining of the 
calculation procedure for shock-wave structure according to the direc- 
tional -mean- free-path method. 

After successful performance of the numerical calculations outlined 
in Chapter VII for shock-wave structure, it may be desirable in the future 
to extend this study in the following several directions: 

(a) inclusion of the effects of boundaries, in some approximate 
manner, in the formulation of the general method; 

(b) inclusion of body forces in the equations; 

(c) calculation of higher order solutions in the shock-structure 
problem; 

(d) attempt to solve the shock-structure problem with collision- 
model III; 

(e) possible extension of the concepts in an approximate manner 
to polyatomic gases; and 

(f) possible extension of the concepts to gas mixtures. 

The first extension suggested, (a), probably will be the most important 
next step in further development of the directional-mean-free-path method. 
Future application of the method to any problem with boundaries, however 
simple, will depend on inclusion of the effects of boundaries in the 
formulation of the directional-mean-free-path approximation. The main 
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value of the directional-mean-free~path method would therefore he realized 
only after such an extension of the theory has been made so that it can 
be applied to problems that may not be tractable by kinetic theory methods. 
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